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PREFACE 


This  report  describes  work  performed  by  ITT  Gilfillan  for  Rome 
Air  Development  Center  (RADC)  under  Contract  So  FJOf-02-  ~S-C-OI  1 9 
The  Project  Engineers  for  RADC  were  Mr.  Vincent  Vannio-la  and 
Russell  Brown.  The  objective  of  this  study  i«i s  to  optimize  radar 
systems  by  using  the  properties  of  the  scattering  matrix,  including 
(Ktiarizalion  and  statistical,  to  design  the  waveform  and  receiver.  The 
report  contains  residts  of  computer  simulations  which  compare  target 
detection  performance  for  derived  dual  chaiuiel  systems  which  incor¬ 
porate  processing  techniques  utilizing  the  scattering  matrix,  against  that 
of  conventional  single  channel  systems. 

Analytical  and  systematic  approaches  are  described  for  five 
techniques  which  lead  to  dual  channel  optimum  processing.  Finally . 
the  statistical  target  and  chaff  models  are  formulated  and  described  as 
they  are  applied  to  computer  simulation  for  system  evaluation. 

The  re/>nrl  was  edited  by  Mr.  Albert  Klein  who  also  contributed  to 
Sections  S.  5.  and  6.  Dr.  David  Hammers  supervised  the  signal  pou  estine 
and  analysis  tasks,  which  were  performed  by  Mr.  Masaaki  fT.-mt  Fuiita. 
Dr.  Gmrge  loanniais.  and  Dr.  Shan  Levan.  Dr.  Jeffery  Bell  wrote 
Section  4  on  modeling  and  simulation  and  contributed  iu  .jCC tl'22lS  I 
and  2.  Mr.  Charles  Lucas  was  responsible  for  the  signal  processor  design. 
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Air  Force  radar  systems  must  operate  against  natural  and  man-made 
interference,  particularly  clutter  and  chaff.  This  program  has  investi¬ 
gated  new  methods  of  target  detection  *n  dual  polarization  radar  systems. 
Contrasting  scatter  in  polarization  between  targets  and  chaff  is 
exploited  to  enhance  target  detection. 

Of  significance  in  this  study  are  design  methods  for  polarization 
coded  transmit  waveforms  and  matched  receivers.  A  simulation  was  per¬ 
formed  for  analysis  of  potential  improvement  factGr.  This  contract 
supports  TPG  P.48,  Surveillance  ECCH. 


SUSSELL  D.  BROWN 
Project  Engineer 
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Section  I 

INTRODUCTORY  SUMMARY 

1.1  BACKGROUND  AND  PROBLEM 

Conventional  radar  signal  processing  techniques  use  Doppler  and  polarization 
information  in  discriminating  target  returns  from  chaff  and  ether  forms  of  interference. 
These  two  discriminants  are  used  independently;  most  often  in  single  channel  system:.  For 
example,  the  polarizations  of  the  transmit  and  receive  antennas  are  fixed  and  determined  b\ 
the  combiner.  On  the  other  bend,  spectral  processing  is  performed  by  the  receiver  processor 
on  signals  that  have  already  been  combined  by  the  antenna.  As  a  result,  such  systems  do 
not  exploit  the  statistical  properties  of  the  target  and  chaff  scattering  matrices,  which  include 
both  spectral  and  polarization  information. 

The  utilization  of  this  information  which  exists  in  the  scattering  matrix  should  be 
more  fully  exploitable  in  2  dual-channel  system,  relative  to  the  more  conventional  single 
channel  system.  Thus  the  basic  tbjective  of  this  effort  was  to  explore  the  merits  of  per¬ 
forming  the  signal  processing  in  2  dual-channei  system,  to  develop  the  optimum  combined 
spectral  ana  polarization  processing  approach.  This  involves  determining  the  optimum  trans¬ 
mit  waveform  and  coherent  receive  polarization  and  time  weights  for  each  pulse  in  a  burst. 
The  criterion  of  optimality  is  to  maximize  the  probability  of  target  detection  for  a  given 
false  2lanr.  rate  attributed  to  the  interference. 

1.2  APPROACH 

The  basis  c-f  our  approach  was  to  formulate  the  system  ir  an  operator  theoretic 
manner.  si-'F.  tliat  the  resulting  approach  would  be  readily  amenable  to  be  validated 
by  computer  simulation.  We  star'ed  by  developing  a  range.' Doppler  spread  dual  channel 
scattering  formulation,  since  this  affords  complete  information  on  the  target  operate..  How¬ 
ever,  \vc  reduced  the  problem  to  a  Doppler  spread  point  target  in  2  dual  channel  formulation 
for  computer  simulation,  and  concentrated  on  this  during  the  bulk  of  the  study. 

Several  system  models  were  developed  and  used  in  the  analysis.  They  include  a 
matrix  receiver  weighting  method,  a  vector  receiver  weighting  method  (essentially  a  subset  of 
the  matrix  method)  a  state  space  approach,  a  Fredholm  integral  approach,  and  a  Stokes 
vector  approach.  (The  Stoke*  approach  is  a  single  pulse  polarization  processing  scheme, 
which  has  no  spectral  processing  capability.  It  was  investigated  in  a  prior  study  and  reported 
hare  only  for  the  sake  of  completeness.)  The  matrix  and  vector  methods,  as  well  as  the 
Stokes  approach  were  simulated  on  a  digital  computer  to  evaluate  their  performance.  The 
state  space  approach  was  carried  to  the  point  of  an  optimum  receiver  design,  given  the 
transmit  waveform  and  target/clutter  scattering  properties.  The  Fredholm  integral  equations 
were  formulated  to  the  extent  of  developing  the  system  of  equations  with  target  clutter 
scattering  matrices. 
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In  order  to  perform  the  simulations,  it  was  necessary  to  have  adequate  models  of 
the  target  and  chaff.  We  chose  the  BQM-34A  as  our  target  model  since  we  have  validated 
simulated  res  plots  with  actual  measurement  data.  This  targe!  model  is  deterministic,  and 
consists  of  mathematical  expressions  for  the  elements  of  the  scattering  matrix  which  are 
functions  of  the  target  orientation  with  respect  to  the  radar,  and  to  the  carrier  frequency  of 
the  emitted  signal. 

Chaff  was  modeled  as  a  collection  of  dipoies  with  selectable  preferred  orientations. 
The  dipole  cloud  was  further  described  by  Doppler  and  dipole  rotation  statistics.  From  the 
distribution  of  the  dipole  orientations  and  the  spectra!  statistics,  a  theoretical  covariance 
matrix  was  derived.  This  matrix  was  then  used  in  conjunction  w;th  the  chaff  model.  Sim¬ 
ulated  signatures  were  derived  from  the  target  model  by  sampling  its  scattering  matrix  at 
various  aspect  angles  corresponding  to  a  trajectory  .  The  simulated  chaff  signatures  were 
derived  from  the  same  sample  rates.  Once  these  simulated  signatures  were  generated, 
correlation  functions  (covariance  matrices)  were  created.  Tne  optimum  dual-channel 
waveforms  and  receivers  were  then  computed  through  application  of  the  various  techniques 
developed  (matrix,  vector,  and  Stokes  approaches). 

The  figures  of  merit  used  in  the  evaluation  were  probabilitv  of  detection  and 
target-to-chaff  ratio  improvement.  The  various  optimum  dual  channel  systems  were  simulated 
and  compared  with  single  channel  simulation  results,  and  with  the  single  pulse  Stokes 
formulation  results. 

1.3  RESULTS  AND  CONCLUSIONS 

The  most  significant  results  are  plots  of  probability  of  detection  vs  signal- to-duff 
power  ratio  generated  from  simulation  outputs  which  compares  dual  channel  performance 
against  single  channel.  An  example  is  shown  in  Figure  1-1.  Here  we  see  she  probabiiitv 
of  detection  vs  input  signa'-to-chaff  ratio  for  a  dual  channel  system  using  matrix  receiver 
weighting.  Three  curves  are  also  shown  for  conventional  single  channel  systems.  The 
conventional  systems  are  assumed  to  have  the  same  polarization  during  both  transmit  and 
receive.  Vertical,  horizontal,  and  right  circular  polarizations  are  shown.  Each  of  the  con¬ 
ventional  systems  is  assumed  to  have  a  filter  matched  to  the  transmit  waveform.  The  dipole 
distribution  of  ±45  degrees  means  that  the  dipoles  in  the  chaff  cloud  are  uniformly  distributed 
between  =45  degrees  from  the  horizontal;  yielding  a  part  ally  preferred  orientation  in  that 
direction.  The  pulse  sample  rate  was  3  msec,  such  that  the  target  model  scattering  matrix 
is  reasonably  iecorrelated  pulse-to-pulse  (approaching  a  Sverling  11  model).  Vhe  curves  show 
that  the  probability  of  detection  is  far  greater  in  the  dual  channel  system  than  in  an>  of  the 
three  conventional  single  channel  systems. 

AH  of  the  simulation  results  reported  upon  yield  to  the  same  basic  conclusion ; 
tha>  is.  the  exploitation  of  the  polarization  and  spectral  spread  scattering  properties 
of  the  target  and  cnaff  in  a  dual  channel  matrix  con  figuration,  provides  superior 
target  detection  performance  compared  to  conventional  single  channel  matched  filter  systems. 
Tire  degree  of  target  detection  enhancement  is  a  strong  function  of  the  pulse-lo-pulse 
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Figme  1-1.  Dual  vs  Singh  Channel  RT3  Target 


correlation  property  of  the  target  scattering  covariance  matrix.  A  Swelling  II  type  target 
which  is  uncorrelated  pulse-to-pulse.  will  perform  better  than  a  Swelling  I.  which  is 
completely  correlated  pulse-to-pulse.  In  the  latter  case  however,  dual  channel  performance 
is  still  better  than  the  conventional  single  channel. 

In  view  of  the  results,  it  is  recommended  that  further  study  be  performed  to 
more  fully  define  the  practical  advantages  of  the  dual  channel  configuration.  This  would 
involve  application  of  a  typical  surveillance  radar  implementation  incorporating  MTl  type 
processing  in  conjunction  with  the  polarization  discriminant.  A  preliminary  design  of  the 
dual  channel  system  would  be  undertaken  to  provide  qualitative  cost  assessment  against 
single  channel  systems. 

Further  analytical  work  is  also  recommended  to  develop  the  state  space  approach 
for  waveform  design  and  computer  simulation.  Finally,  the  range  spread  pjoperties  of  the 
target  should  be  considered  as  an  added  discriminant  toward  optimizing  targe*  detection 
in  clutter. 
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Section  2 


PROBLEM  FORMULATION 

The  purpose  of  this  section  is  to  mathematically  formulate  the  problem  with  the  aid 
of  scattering  operator  theory.  Using  this  approach,  we  treat  the  target  as  a  system  which  is 
characterized  by  a  scattering  matrix  (operator).  The  function  of  this  scattering  matrix  is  to 
operate  on  incoming  or  incident  radar  signals  (waves)  to  produce  backscattered  or  output 
signals.  This  represents  a  linear  coherent  dual  or  two  port  system  which  can  account  for 
the  horizontal  and  vertical  polarizations.  If  the  transmitted  radar  signal  is  a  two-clement 
polarization  vector  which  :$  operated  on  by  the  2x2  target  scattering  matrix,  ih-n  the  cor¬ 
responding  two-element  vector  scattered  back  to  the  radar  is  considered  to  be  the  radar 
receiver  input.  This  receiver  input  is  then  operated  on  by  the  radar  in  processing  the  signal. 

A  block  diagram  cf  tiiis  total  process  is  depicted  in  Figure  2-1.  Here  we  see  dual 
channel  signal  modulation  during  the  transmitting  process.  Likewise,  the  dual  channel  receiver 
process  is  shown  operating  on  the  backscattered  signal  from  both  the  desired  target  and  the 
•mdesired  clutter.  The  basic  problem  treated  in  this  study  is  to  determine  the  optimum 
combination  of  dual  channel  signal  waveform  (including  the  transmit  antenna  polarization)  and 
the  dual  channel  signal  processing  (including  the  receive  antenna  polarization).  By  optimum 
we  mean  that  which  yields  the  maximum  probability  of  detection  for  a  given  false  alarm  rate. 

The  approach  is  to  solve  for  the  dual  channel  transmit  wavetorm  loptimum)  which 
provides  the  maximum  separation  of  received  target  and  clutter  covariance  functions.  This 
requires  a  vector  fonn  of  the  scattering  operator  to  effectively  manipulate  the  mathematical 
functions  inherent  in  the  covariance  matrices,  so  as  to  pi>perly  isolate  the  clutter  return 
from  the  target  plus  clutter  return.  Once  this  is  obtained,  the  receiver  function  is  derived  to 
maximize  the  probability  of  detection. 

Since  a  key  element  in  designing  the  optimum  dual  channel  system  is  knowledge  of 
the  target  and  chaff  covariance  matrices,  we  will  also  develop  a  methodology  in  this  section 
for  computing  these  matrices  for  use  in  optimizing  the  system  design. 

2.!  POLARIZATION  PROPERTIES  OF  TARGETS  AND  CLUTTER 

It  should  be  fairly  obvious  that  more  information  about  a  target  and  background 
exists  in  a  dual  channel  polarization  sensitive  system  than  in  the  conventional  single  channel 
system.  Since  our  goal  here  is  to  include  dual  channel  scattering  characteristics  into  a  system 
design,  let  us  briefly  review  a  few  important  known  polarization  sensitive  characteristics  3bout 
targets  and  clutter. 

In  the  area  of  dutter  suppression  we  know  that  polarization  sensitivity  of  rain  ha^ 
been  well  established,  and  that  drcuhif  polarization  can  be  used  to  cancel  isotropic  rain.  More 
recent  approaches  (Regain  (1975).  Nathanson  (1975)j  have  improved  cm  this  by  estimating 
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U(t)  ;  BASIC  UNMODULATED  CARRIER  SIGNAL 

M(t)  ^  DUAL  CHANNEL  TRANSMITTED  SIGNAL  MODULATION  OPERATOR 
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i  TIME-RANGE  TARGET  SCATTERING  OPERATOR 
Hc(t.r)  =  TIME-RANGE  CLUTTER  SCATTERING  OPERATOR 
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G(t}  DUAL  CHANNEL  RECEIVER  SIGNAL  PROCESSOR  OPERATOR 

2(t)  DETECTED  SIGNAL 

Figure  2-i.  Dual  Channel  Scattering  System 


the  ellipticity  of  nonisotropic  rain  and  then  adaptively  varying  the  horizontal  and  vertical 
polarized  receive  channel  to  enhance  suppression.  For  target  enhancement,  it  has  been  found 
that  pulse-to-pulse  variation  of  polarization  of  a  transmitted  sisnaJ  will  tend  to  decorrelate 
target  returns  in  a  manner  similar  to  that  achieved  by  frequency  variation.  Poelman  (1978) 
analyzes  this  problem  relative  to  simple  targets  and  shows  that  considerable  improvement  can 
be  achieved  in  target  detectability  through  polarization  agility.  The  S3me  concept  has  also 
been  applied  to  improve  angle  estimation  in  a  multipath  environment  {Ewell,  et  al.  (1971)j. 

The  pulse-to-pulse  correlation  phenomena  for  targets  and  rain  may  also  be  extended 
to  other  types  of  gro  .id  ana  sea  clutter.  For  example,  backscattcr  from  manmade  cluster, 
such  as  buildings  (cities),  tends  to  exhibit  much  variability  (10  to  20  dB)  between  H  and 
V  responses,  whereas  mountainous  areas  tend  not  to  exhibit  such  deviations.  Angies  of 
incidence  and  wavelength  influence  this  heavily.  More  variability  exists  at  smaller  aspect 
angles  and  lower  frequencies.  For  example,  at  angles  of  incidence  less  than  10  degrees,  the 
variation  in  the  VV/HH  ratio  has  been  measured  at  C-band  to  be  about  14  cB  {Daley  (1968), 
city  of  Phoenix  j. 

Chaff  drops  (of  interest  in  this  study)  are  characterized  by  metallic  type  dipoles 
which  are  cut  to  resonate  at  the  radiated  frequencies  of  the  radars  they  are  designed  to 
confuse.  As  with  target  scattering  surfaces,  maximum  backscattcring  occurs  when  the  majority' 
of  the  dropped  dipoles  are  oriented  parallel  to  the  F-ficld  of  the  transmitted  signal.  After 
the  chaff  has  been  dropped  dispersion  occurs  (due  either  to  air  conditions  or  by  design)  so 
that  the  principle  plane  of  the  dipoles  in  the  chaff  cloud  changes  with  respect  to  space  3nd 
time.  The  changes  that  occur  can  be  correlated  with  wind  characteristics.  Relative  to 
sampling  time,  it  is  possible  for  a  polarization  diverse  radar  to  make  estimates  cf  the 
polarization  properties  ?x.d  thus  adapt  in  such  a  way  as  to  maximize  the  target-to-clutter 
ratio.  The  bases  for  utilization  of  these  polarization  proj*crtics  are  formulated  in  the 
following  sections. 

2  2  SCATTERING  SYSTEM  FORMULATION 

The  relationship  between  the  polarization  sensitivity  of  targets  and  background  as 
applied  to  radar  systems,  car.  he  described  by  a  dual  channel  (or  2-peT-D  opeiator  as  shown 
in  Figure  2-1.  In  this  formulation  we  consider  changes  in  the  wideband  range  dependent 
scattering  properties  as  well  as  longer  time  ( pube-to-pube)  scattering  properties  due  to  motion. 
As  a  result,  both  interpube  and  intrapulse  signal  modulation  and  the  corresponding  optimum 
receiver  can  be  derived  to  maximize  the  signal-1  o-background  ratio. 

Referring  to  Figure  2-1.  we  find  that  the  transmitted  signal  is  the  2  x  1  complex 
vectc-rX(t)  where  elements  xj(t)  -  the  vertically  polarized  component  -  and  \-><t)  -  the 
horizontally  polarized  component  -  are  functions  of  the  time  variable  t.  <  t  <  ». 
Similarly,  the  received  signal  will  be  a  2  x  1  vector  denoted  by  Yrt). 
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The  scattering  matrix  -  which  relates  the  relative  amplitude  and  phase  of  the 
incident  and  backscattered  electric  field  vectors  -  of  a  complex  radar  target  is  given  by  a 
2x2  matrix  of  the  form: 


H(t.r)  = 


[" hj  j(t,r) 


hvjtU 


h.  i(t,r)l 

1  —  I 


hntt,r) 


Here,  the  entries  hjj(t,r),  i,j  =  1,2,  are  taken  to  be  complex  Gaussian  processes.  •  is  the  time 
variable,  and  r  is  the  space  variable  —  measured  in  units  of  two-way  travel  time  to  the 
reflection  point  on  the  target. 

From  the  scattering  approach,  a  complex  radar  target  then  admits  the  input-output 
dual  channel  or  two-port  description: 


Y(t)  =  /  Hft,r)  X  ft  -  rkir, 


=  f  H(t,t  -  r)  X  (t)ut. 


These  equations  form  the  familiar  linear  relationship  between  the  input  signal  (X>  and  output 
(Y)  when  operated  on  by  a  system.  At  this  point  it  is  convenient  to  introduce  a  scattering 
rector  which  has  tire  same  elements  as  the  matrix  H(t.r)  and  is  defined  by: 


S(t.r)  -  J h j  jtt.r),  hj-»(t,r),  h->j(tjr),  ti22fM’M 


where  ~  denotes  matrix  transposition. 


It  then  follows  that  <  2-2)  can  be  written  as 


f  ' r 

'J 

•/-oo  Q 


jxjft  -  r)  x->(t  -  r)  0  0  ~|  |  hjj(t.r)  | 


1 

Xj(t  -  r)  X2<t  -  r) j 


j  ^1  ^  j 

j  i 

j  h2l(t.r)  j 


h->->{t.ri  • 
_  ”  _l 
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or. 


where 


Y(t)  =  /  Xf t  -  r)S(t.r)dT  (say). 


X(t-T)  = 


r  x^t 
L  o 


-  r)  x~>(t  -7)  0  On 

0  XjU  -  -)  x^U  -  7) 


(2-5) 


(2-6) 


Smflariy,  (2-3)  can  be  written  as 

Y(t)  =  r  X(r)S(U  -  rkir.  (2-7) 

-00  2 

Equations  (2-5)  end  (2-7)  now  define  the  received  signal  vector  Y  in  terms  of  the 
transmit  matrix  X  and  th-_  sca'tering  vector  S. 


2.3  COMPUTATION  OF  RETURN  SIGNAL  COVARIANCE  MATRIX 


2.3.1  Definitions  of  Important  Matrices 

Since  knowledge  of  the  covariance  matrix  of  the  target  and  clutter  returns  are  of 
major  importance  in  designing  the  optimum  system,  wc  will,  in  this  subsection,  define  the 
covariance  matrix  and  devise  two  approaches  for  computing  it.  Here  we  define  a  general 
expression  for  the  covariance  matrix  Ky(t,w)  of  the  received  signal  Y(t)  -  which  is  a  complex 
Gaussian  process  with  zero  mean.  Thus,  let  E(*)  denote  the  expectation  operation  and 
*  stand  for  the  complex  conjugate  transpose,  then 


KyXw)  =  E[Y(t)Y(w)*l 

E{yj(t)yj  <»>*}  -lyj(t)y2V'v)*! 

Efysftjyjfw)*)  Elyiftfeifw)*) 


(2-8) 


Before  proceeding  to  derive  the  two  methods  of  computing  the  covariance  matrix,  it 
is  instructive  to  define  two  new  matrices  which  have  the  .'ame  information  as  the  scattering 
matrix  H(t,r),  Equation  (2-1).  They  are 


?nd 


MU.-) 


N(t.r> 


hj  j(f.T)l7 
h->j(t,T)I-> 

H(t,r) 

°2 


n12^’r^2 

h-»-i(t.r)I-> 

°2 

H(t.r) 


(2-9) 


(2-10) 


where  f?  is  the  2  X  2  identity  matrix,  while  0?  is  the  2X2  zero  matrix. 
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tV 


2.3.2 


Fits!  Method 


if 


1 

5 

S 

$ 

IS 


We  now  present  the  first  method  of  computing  Kv;t.w)  using  the  target  scattering 
vector  S(t,r).  For  this  wv  define  the  scattering  vector  covariance  matrix  K$(t.r,w,o)  as 


K^O.t.w.o)  =  EJS(t.r)  •  S(w .0)*! 


Next,  we  have 


Mtt.rJ 


pi 

[ 

• 

0 

0 

Uj 

«  hi:  °]  P  I  Pn’« 


h2i  0  h22  0  jO  h2] 

I 

0  h-»»  0  h-*->  i  i  h-»->  ! 

-•  — J  ! _ <_  — _i 


!  =  Slt.r) 


where  the  hjj’s  are  functions  of  i  and  r.  Similarly,  it  is  easy  to  verify  that 


(2-11) 


{:-l2l 


E* 


I 

I 

£ 

1 

-s 

1 

5 

1 

I 


Therefore 


0  i 

— —  | 

N(t.r)  =  S(t,7) 

0 


S(i.r)  *  = 


Mn.rl  i  |1  00  1]  N(w.o» 
0 


(M3) 


S(i.r)  *  S(w.o)*  =  M(t.7>  * J_*  Nfw.o>  =  Nlt.r)  •  Mtw.o)* 
where  —  indicates  complex  conjugate  and 


(2-14) 


where 


I  = 


i  0  0  1 
0  0  0  0 
0  0  0  0 
10  0  1 


Consequently. 


(2-1 5 > 


Kja.r.w.o)  =  EI(M(!,t)  - 1  *  N(w,o)J  =  EfNu.rl  •  X*  M(v.o)*}  (2-16) 


It  then  follows  from  Equations  (2-5)  and  (2 S)  that  the  covariance  matrix  Kv(t.w)  of  the 
received  signal  Y(t?  is  given  by 

K.,{t,w)  =  r  r  XO  “  ">  *  Kc(t,T,w.o)  •  X(w  -  a)*  dodr  .  (2-17) 

'  — CO  — DO  *  - 


where  K<j(->  is  the  target  scattering  vector  covariance  matrix  and  X(-)  is  the  transmitted  signal 
matrix  defined  ir.  Equation  (2-5). 

23.3  Second  ,\iethori 

The  second  method  does  not  lead  to  a  direct  computation  of  the  covariance  matrix, 
instead,  a  -t-ctor  is  computed,  which  has  the  same  elements  as  those  found  in  the  covariance 
matrix.  This  covariance  vector,  as  we  will  call  it.  is  4x1  and  defined  by 

KyJt.W)  =  {Kyjj(t,w),  Kyp(t,w).  Ky,j(t.W).  Ky-,,(*.W>I  f2-18) 

where  Ky~(t.w)  are  the  elements  of  the  covariance  matrix  in  Equation  (2-S). 

If  we  further  define  the  vector 

X(t-r,w-o)  =  Ixj(t-r)xj(w-o).  xj(t-T)x2(w-o), 

xifi-rxitw-cr),  xift-rlx^fw-a)!  (2-19) 

the  covariance  vector  is  given  by 

rc..(t,w)  -  /  /  ElM(t.T)  •  Nfw,o)j  X( t-r.w-o )  dodz 

-  — oo  —oo 


(2-20) 


This  may  be  written  as 


Ky(t,W) 


Kjj(t,r.w,o)  X(t-r,w-c)dedr, 


(2-21) 


where 


Kj|(t,r,w,'7)  =  t{M(t,r)  •  S(w.c)J  (2-22) 

is  defined  as  the  scattering  matrix  covariance  matrix. 

If  the  entries  of  H(t,r)  and  H(w,o»  are  uncorrelated  for  7  ~  a  that  is.  when  the 
return  signals  from  different  points  on  the  target  are  uncorrelated  then 


Ky(t,W) 


OO 

r 

j 

-00 


fi[  M(t,r)  -  N(w.r)j  X{t-r,w-r)dr. 


(2-23) 


It  b  important  to  note  that 


M(t,r)  -  N(w,g) 


h 


h 

h 


il^ll 

hl,h,2 

hl 2hl i 

hj  -»h  j  -> 

1 1^21 

hj 

h12h21 

h 1 2h22 

2l*Ml 

*121*112 

h22*Ml 

h-»->h|  -> 

2I**2l 

h-»  j  h  -»-> 

h->->h->j 

h  -*  -»h  ->  -* 

N'(w.o)  *  M(t.r). 


(2-24) 


where  the  h::  arc  functions  of  t  and  7,  while  the  h:;  arc  functions  of  w  and  n. 

h 


23.4  Further  Notes 

It  is  important  to  note  th2t  in  the  dual  channel  description  above  (Equations  (2-2) 
and  (2-3))  a  radar  target  is  characterized  by  a  linear  system  whose  impulse  response  matrix  is 
the  target  scattering  matrix  (H(t.r)  (Equation  (2-5 ».  However  in  this  case,  the  entries  of 
H(t.r)  arc  random  processes  and.  in  general,  arc  unknown.  Thus  from  a  system  theoretic 
viewpoint,  we  have  a  linear  random  system  and  oar  problem  is  basiraHv  that  of  "identifying" 
the  random  entires  hjjft.r)  of  its  impulse  response  matrix  H(t.r).  To  thb  end  we  have 
transformed  Equation  (2-2)  into  Equation  (2-5)  by  introducing  the  scattering  vector  Sft.r) 
from  Equation  (2-4)  and  the  transmit  matrix  _X(t-r)  from  Equation  (2-6).  This  is  justified  by 
the  fact  that  Equation  (2-2)  when  written  out  in  full  becomes 


yj(t) 

90 

—  n 

hjj(t,T)  sj(t-r)  +  hj2(t,r) 

y?(t) 

_ 

=  / 

-OO 

h->  jft.r)  xj(t-r)  +  h->->(t,r)  x-*(t-r) 
L~  _J 
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Then  each  of  tne  scalar  terms 


J°°  Mt-')  Xw  (t-r)  dr  .  i.  j.  k.  =  1.  2.  ( 2-26) 

-ao  v  K 

is  actually  the  output  of  a  linear  scalar  (i.e ,  single  input-output)  system  whose  impulse 
response  function  is  the  scalar  random  function  h;j(t,r)  and  whuse  input  is  the  scalar  function 
\j.(t-r).  Since  the  system  is  iinear.  Equation  (2-26)  can  also  be  regarded  as  the  output  corre¬ 
sponding  to  the  input  hjj(t,r)  (in  time  and  space)  of  the  system  whose  impulse  response  is 
xj.(t-r).  He  net  Equation  (2-2)  becomes  Equation  (2-5)  in  which  the  four  scalar  inputs 
hjj(t.r)  ,  i,  J  =  1.  2,  are  arranged  into  the  4  X  1  scattering  vector  Sit.r).  Hcr.ce,  it  is 

easy  to  see  that  Equation  (2-5)  is  a  linear  transformation  from  a  4-dimens:ona!  space  into  a 
2-dimensionai  space.  Since  the  input  vector  X<t-r)  is  a  2  X  1  vector.  Equation  (2-2?  is  just 
a  linear  transformation  from  a  2-dimensional  space  to  a  2-dimensional  space. 

The  idvantages  of  Equation  (2-5)  are  many  fold  first,  the  problem  of  identifying 
the  matrix  H(t,r)  now  becomes  that  of  identifying  the  vector  S(!.r)  from  the  output  data 
Y(t).  This  allows  one  to  use  techniques  of  signal  design,  processing  and  simulation  in  radar 
problems.  Secondly,  in  computing  the  output  covariance  matrix  Kv(t.w)  in  Equation  (1-17). 
we  have  relied  heavily  or  the  covariance  matrix  K^t.r^v.o)  in  Equation  (2-11)  for  the  vector 
process  S(t,r).  This  in  turn  is  expressible  in  terms  of  the  two  matrices  M(t.r)  and  N(w.o)  of 
Equations  (2-9)  and  (2-10)  respectively.  These  matrices,  as  we  have  shown,  aie  on  the  one 
hand  directly  related  to  the  vector  S(t,r)  and  therefore  allow  one  to  compute  Kv(t.w)  from 
simulated  or  experimental  data,  and  also  on  the  other  hand  lead  to  the  type  of  "‘spectral 
polarization."  of  Kjt't.r.w.o)  in  Equation  (2-16)  and  KH(t.T.w.o)  in  Equation  (2-22). 

finally,  it  is  noted  that  for  the  case  of  an  isotropic  target,  hp  =  h->  j .  the  matrices 
M(t.r)  and  N(t.r)  arc  both  symmetric  and  so  is  the  matrix  Kjj(  ).  whereas  K^i  1  aid 
are  always  hermitian  matrices. 

The  two  methods  discussed  earlier  arc  used  in  Sections  5.13  and  3.1.2  respectively. 
It  is  appropriate  to  comment  on  the  fact  that  due  to  the  i.uture  of  the  Generalized  Matrix 
Approach,  the  second  method  cannot  be  used,  whereas  ;he  Generalized  Vector  Approach  may 
use  either  methods.  However,  from  the  computational  aspects,  the  second  method  was 
preferred. 


m 


m 
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2.4  COVARIANCE  MATRICES  OF  POINT  TARGET  IN  CHAFF 

\W  now  apply  the  result?  found  a  bow  to  the  case  of  a  point  target.  We  have  in 


tnis  case 


hjj(t)  hp(t) 

=  *6(t  -  Tj 

hijft)  h  t ) 


q)  -  H(t)  •  6(7  -  Tq) 


where  tq  is  the  Two-wav  travel  time  to  the  point  tarset. 

U  is  seen  from  Equations  (2-9)  and  (2-10)  that 


fclUON 

t.r>  = 


En(t)l 


h-»-)(tjl  -> 


*S(t  -  Tq)  =  M(t)  *  5(7  -  Tq) 


(2  28) 


1 


H(i)  02 

^‘-7)  =  *  5(7  -  Tq)  =  \'( 

0,  H(t) 

L- 


t)  *  5(7  -  Tq) 


Therefore,  from  Equation  ? 2-23) 


(2-29) 


Kh<!.t.w,o)  =  Ef\i(t)  -  NVw))  •  5(t  -  tq)  •  5(o  -  Tn». 


nd  from  Equation  (2-16) 


(2-30) 


KS<t,7,W,0)  =  E|M(i)  •  1*  N(W)1  -  6(7  -  Tq)  -  5(0  -  Tq). 


Then,  using  Equation  (2-23)  we  find 


*y(;.w)  -  E{\!a)  ♦  N(w) j  -  X(t  -  tq,  w-  -q). 


Similarly,  from  Equation  (2-17) 


KyO.W)  -  X(t  -  Tq)  *  Ej\!(t)  •  j  •  N'(\v)j  •  X(W  -  Tq)*. 


(2-31) 


(2-32) 


0-33) 


Let  HjO.t)  Hj(t)  ♦  6(r  -  tq)  be  the  scattering  matrix  of  a  point  tareet  surrounded 
by  a  chaff  cloud  -  whose  scattering  matrix  is  denoted  by  H^t.r).  As  in  the  above,  let  X't) 
and  Y(t)  be  the  transmitted  and  received  signals,  respectively,  then 


I 

1 

3 


Y{t)  =  HT(t)  *X(t  -tq )+_T  -X(t  -  t)  •  dr 

=  Vj<t)  +■  Y^t) 

Let  Ky(t,\v)  be  the  covariance  matrix  of  Y(t),  we  have 

Ky(t,w)  =  Ei(Yj<t)  +  Yc(t))  (Yj(w)*  +  Y^w)*)}. 

Therefore,  since  Yy  and  Y^-  are  uncorrelate'? 

Kyft-w)  =  E[YT<t)Y7<w)*l  +  E[Y0J)Yc(w)*j. 

=  Ky^t.w)  +  Ky  (t,w),  (2-34) 

where  the  point  target  covariance  matrix  KV  (-(t,w)  is  computed  from  Equation  (2-33),  while 
the  chaff  cloud  covariance  matrix  KyrXt,w)  can  be  derived  from  Equation  (2-17). 

It  is  also  evident  from  Equation  (2-34)  that 

Ky(t.w)  =  KyT(t.w)  +  «yc(t,w)  (2-35) 

where  Ky  (t,w)  is  g»ven  by  Equation  (2-32)  and  kv  (t,w)  can  be  computed  from 
Equation  (2-21).  *'■ 


As  noted  in  the  following  sections,  the  received  signal  covariance  functions  become 
key  dements  in  the  design  of  the  transmit  waveform  and  associated  receiver  weighting  for 
maximizing  the  output  s/c  ratio. 


Section  3 


OF1IMAL  TRANSMIT  WAVEFORM  AND  RECEIVER  DEVELOPMENT 

The  previous  section  described  the  manner  whereby  radar  targets  and  clutter  can  be 
characterized  by  scattering  operators.  Associated  with  these  scattering  operators  are  scattering 
vectors,  which  can  be  regarded  as  stochastic  vector  processes.  With  this  description,  the 
signals  received  from  the  target  and  clutter  can  be  formulated  as  the  outputs  of  a  determinis¬ 
tic  system  whose  inputs  are  the  scattering  vectors,  see  Figure  5-1.  Since  a  major  study  task 
was  to  maximize  the  probability  of  target  detection  for  a  given  false  alarm  rate,  the  more 
immediate  problem  becomes  one  of  transmit  waveform  and  associated  receiver  design  to  meet 
this  objective.  Given  then  the  scattering  vector  descriptions  of  both  target  and  clutter,  it  was 
necessary  to  rind  the  optimal  deterministic  transmit  waveform  and  optimum  receiver  weighting 
for  maximizing  the  target  return  signal-to-chaff  plus  receiver  noise  ratio. 

Five  v  uiants  of  the  scattering  theory  approach  were  investigated.  Two  of  these 
methods  fall  into  the  domain  of  finite  dimensional  discrete  time  and  frequency  space.  These 
were  more  fully  devef  ■«>ed  than  the  other  methods  in  terms  of  obtaining  both  theoretical  and 
simulation  resuits  for  waveform/receiver  design  and  associated  detection  performance.  The 
state  space  and  matched  filter  Fredholm  Equation  approaches  were  both  mathematically 
formulated,  but  they  were  not  simulated  on  a  computer. 
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Figure  3-1.  Simplified  System  MoJ*I 
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finite  dimensional  space 

The  two  methods  related  to  finite  dimensional  discrete-time  and  frequency  space 
were  developed.  They  are  herein  designated  as  the  vector  and  matrix  approaches  The 
matrix  approach  has  the  potential  for  utilizing  all  of  the  information  a-.ailabie  'within  the 
dimensional  space  constrained  by  both  the  transmit  signal  and  target/chaff  scattering 
properties.  Assume,  for  example,  that  the  target  and  chaff  sc-ttering  sectors  each  consist  of 
four  variables  (hj  j,  bp,  h->j,  Iit-i),  and  that  the  transmit  signal  is  ~omp.*seii  of  N  poises, 
each  containing  separable  V  and  H  polarized  components.  The  four  vLment  scattering 
rector,  in  combination  ith  tile  N’  dual  polarized  transmit  pulses.  •.-•.iaMs'i.ei  a  4\  \  4N 
receiver  weighting  matrix  for  optimum  processing  of  the  return  signal.  As  ibsequently 
shown,  the  weighting  matrix  results  from  eigenvector  solutions  ’’'voicing  th.  target  and  .nail 
scattering  vector  covariance  matrices.  However,  since  the  observable  data  consist  or  only 
2N  separable  signal  returns,  the  receiver  weighting  dimension  must  be  reduced  to  a  3  x  4N 
matrix.  The  matrix  elements  then  consist  of  the  2N  eigenvectors  as%oci:ited  with  the  2N 
largest  eigenvalues  resulting  from  the  scattering  vector  covariance  matrix  solutions.  This 
resultant  weighting  matrix  maximizes  the  signal-to-clutter  ratio  (s  c)  under  conditions  shown 
•:n  Section  3.1.2. 

By  contrast,  the  vector  approach  utilizes  only  the  one  eigenvector  which  is  associated 
with  the  largest  eigenvalue.  Consequently,  the  receiver  weight  is  then  a  4\-ck-mvni  vector. 
Although  the  vector  approach  can  be  considered  as  a  subset  cf  the  matrix  implementation,  it 
was  developed  from  a  somewhat  different  mathematical  approach  and  provided  valuable  insight 
into  subsequent  development  of  the  matrix  method.  Additionally,  the  sector  approach  is 
simpler  to  implement  both  from  a  software  and  hardware  standpoint,  which  can  prove  to  be 
o!  further  benefit  in  the  realization  of  actual  systems,  even  though  performance  is  not  as  eood. 

Associated  with  the  vector  and  matrix  approaches  is  the  concept  of  the  maximum 
likelihood  receiver,  from  which  the  false  alarm  and  target  detection  probabilities  were  derived. 
The  likelihood  implementation  is.  therefore,  described  in  the  following  sections  for  complete¬ 
ness.  It  should  be  realized  however,  that  the  design  of  the  transmit  waveform  and  receiver 
weighting  function  n-.cd  not  depend  upon  the  maximum  likelihood  concept. 

3.1.1  Generalized  Vector  Approach 

The  vector  approach  is  described  below  relative  to  the  interaction  of  a  coherent 
pulse  train  with  a  target  characterized  by  a  polarization  scattering  matrix  possessing  both 
temporal  3.nd  range  variations.  The  analysis  presented  here  is  general  and  can  be  applied 
to  both  the  coherent  pulse  train  anJ  the  single  coded  pulse  wave-form.  In  wha:  follows,  we 
assume  that  the  radar  system  under  consideration  employs  two  orthogonal  polarizations  on 
transmit  and  receive. 
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The  envelope  of  transmitted  vector  is 


X(t)  = 


xj(t) 


X'.(t) 


(3-1) 


and  the  target  scattering  matrix  is  represented  by 


11*1^1.7)  — 


h;2(t.r) 


hj  i(t.ri 

h->j(t,7)  hn(t,T) 

The  taraet  covariance  matrix  is  a  4  x  4  matrix  riven  bv 


^HT<,_1,t’  =  *^ij  T,>  '  =  1.  -  -  .  4.  j  =  1.  ...  4. 
where,  referenced  to  equation  (2-24); 

kj,(t-u;T)  =  Hfhj  j  (t."%  h  j  j  (u.  t)J 

k  j  t  (t-u:  t)  =  HJh j  j  (t,r)  h  j  2  (u,r)} 

kj3(t-u:r)  =  H{hJ2(t.r)  hj  j  (u,t)1,  etc. 

and  where  H(  *  J  =  expectation  operator.  (Note  that  when  the  scattercr  is  isotropic 
hj2  =  h->j)  then  Kjj  is  symmetric.) 

For  3  receiver  weighting  function 

A 

Wi;)  = 

the  video  amplitude  at  the  output  of  the  receiver  is  given  by 

Y(t)  =  £  dt,  W*  (t;-t)  £  H-jft,.  r)  X  (tj  -  r)  dr 
From  above  we  find  that  the  average  output  power  is  proportional  to 
EJV(i)  V*(t)j  =/°°dtj  /“dl2/°d71  /°°dr2  E|W«(trt)HT(i1.7,)X(ir7,)X*<t2-72)HT^t2.r2)W<t_ 


(3-2) 


(3-3) 


(i.e.. 


(3-4  > 

(3  3) 

-0! 

(3-6) 
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If  we  assume  that  the  elements  of  ll-Ht..  —  j »  and  Hytts.  r~)  are  uncorrelated  for  rj  =*  r-»,  then 
we  can  write  (3-6)  as 


o©  oo 


Pm  =  F.[V(t)  V*(t)l  =  J  dt,  j  dts  \v*{t,-t)  W(!,-ti  (3-7) 

— OO  1  -OO  —  1  >11-  - 


where  Kv(tj-t->)  is  a  2  x  2  cocariar.ce  matrix  whose  ekments  are  riven  by  (3-8)  below  as 
previously  derived  in  liquation  (2-20). 

j  x,(tl-r)  x  jits-T)] 


KvT<llt2‘  = 


KyTl  5ftl12i 
KvT12<Il12) 


Kvyj  j  ^  XjUpr)  x 

=  L, dr  KH7l'rl2- 

-j — » j ( t j ,t-»)  x^ttj-r)  x  jlti-') 


Kvi2^(tj  -t  -> ) _ j 


X">(ti-r)  \  ->(ti— r)j 


Introducing  the  discrete  vector  W.  where  N  is  the  .unber  of  pukes 
W*  =  J\V*  « 1 ».  W*  f2»  .  .  .  W*  *Nm 

I  the  astensk  indicates  complex  conjugate  transpose),  we  note  that  in  discrete  form  PiO» 
obtained  from  liquation  (3-7)  for  t  =  G  is  rivcn  by 


1X0)  =  W*K>tW  ,3-10) 

where  Kvyis  an  N  x  N  covariance  matrix  whose  elements  are  the  2x2  matrices 
Kv  (lit.  jit),  where  it  is  the  interpulse  period  Thus,  using  liquations  1 3-<> )  and  i3-~i.  we 
can  write  the  sjgr.al-to-chaff  ratio  as 

i\i”  HtXi  »X*  Ht*  \V) 


is"  K.  .  W 

\  s 


where  K>c  is  the  N  x  N  clutter  covariance  matrix  defined  ir.  a  similar  manner  as  Kvy  above. 


vector  is  a  function  of  the  target  and  chaff  scattering  covariance  matrices,  and  also  of  the 
transmit  signal  vector.  Although  the  scattering  properties  are  known  stochastic  functions,  the 
weighting  vector  cannot  be  determined  until  the  transmit  waveform  is  specirxd.  The  vector 
approach  as  derived  abor-e,  therefore,  decs  no*  yield  a  direct  solution  for  the  optimal  transmit 
waveform.  An  iterative  procedure  is  possible,  however,  which  does  nrs-aii  hi  ^  cuasi-opnnwm 
transmitter  waveform  design  under  certain  conditions,  as  described  below. 

3.1.2  Vector  Approach  Waveform  Design 

The  approach  of  Rummler  (1966.  ?  967)  is  adapted  herein  to  the  formulation  of  an 
optimum  set  of  transmitter/receiver  operators  for  maximization  of  the  s/c  ratio  for  dual 
channel  systems.  The  generalized  solution  for  the  receiver  weighting  sector  \V  as  given  in 
Equation  (3-14),  will  be  applied  to  the  case  of  a  Doppler  spread  point  target  in  Doppler 
spread  clutter.  For  convenience.,  the  notation  denoting  the  time  varying  properties  cf  the 
transmit  vector  and  scattering  functions  will  be  omitted.  The  signahto-chaff  ratio  of 
Equation  (3-10)  can  be  written  for  the  above  point  target  case  as: 

w*kvT  w 

s/c  =  - 

W*KVCW  (-'-!:>> 

where  and  Kyc  are  the  target  and  clutter  covariance  matrices,  and  where 

KyJ=  E[(HX)(HX)*J  =  X  E  I KH*1  X*  =KHT^  (3-16) 

H  is  the  previously  defined  2x2  target  scattering  matrix.  The  parameter  Ky.  is  defined 
similarly  for  clutter.  It  is  again  emphasized  that  both  Kyy  and  Kv..  are  functions  of  the 
transmit  signal.  X.  Given  K..y  and  Kyc.  the  problem  becomes  one  of  finding  the  value  of  W 
which  maximizes  the  s/c  ratio. 

It  is  clear  that  s/c  (kW)  =  s/c  (W)  where  k  is  any  arbitrary  scalar  constant.  So  we 
may  assume,  without  less  of  generality',  that  \V*  K^.  W  =  1.  Hence,  the  original  u-Kcrsitained 
optimization  is  equivalent  to  the  following  constrained  optimizaicn  problem: 

max  W*  Kvy  W 

subject  to  W*  Kyc  W  =  1  (3-17) 

Applying  Lagrange’s  multiplier  rule  (see  Appendix  C).  we  sec  that  the  optimum  W  and  s  c  must 
^>tisfy  the  following  conditions: 

‘S'T-^max  Sc’  W  =  0 

^c^max  ~  \nax 
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(3-18) 


where 


=  maximum  eigenvalue  of  KyC  *  Kyj 

The  optimum  signal-to-chsf!  ratio  is  then  the  largest  eigenvalue  of  Kvc~*  Kyj.  and  the 
optimum  weighting  sector  i»  the  eigen.ee tor  corresponding  to  this  largest  eigenvalue. 

Note  th2t  the  s/c  ratio  was  maximized  for  a  given  transmit  vector,  which  had  to  be 
chosee  a  priori.  The  choice  would  generally  be  made  based  upon  experience  dealing  with 
specific  target /clutter  situations.  According  to  Rummler,  this  initial  choice  could  be  the 
unweighted  uncoded  vector.  In  any  cast,  it  may  be  possible  to  further  improve  the  s/c 
through  an  :terathe  procedure  devised  by  Rummler.  This  is  based  upon  the  fact  that  *ne 
s/c  ratio  is  unchanged  if  the  transmitted  waveform  and  the  receiver  response  are  conjugated 
and  interchanged.  This  is  dons  by  using  the  receiver  weights  from  one  iteration  as  the 
transmitted  signal  vector  for  the  next  iteration,  for  which  the  procedure  will  converge  to  a 
maximum  s'c. 

Although  the  solution  converges  to  an  optimum  s'c  ratio  there  is  no  guarantee  that 
a  better  solution  could  not  have  been  reached  if  the  procedure  had  been  ••tsrted  with  a 
different  initial  waveform.  Rummler  has  noted,  however,  that  the  s  c  ratio  convergence  does 
not  appear  to  be  strongly  dependent  upon  the  choice  of  initiai  vector.  Nevertheless,  the 
vector  to  which  the  procedure  converses  is  strongly  dependent  upon  the  initial  choice. 


This  iterative  procedure  was  developed  for  computer  solution  of  transmit  signal  ano 
receiver  weighting  configurations  for  optimization  of  the  s/c  ra*  o  for  both  dual  and  single 
channel  configurations.  Results  of  this  Technique  are  report-  a  von  in  Section  5.  A  summary 
of  the  iterative  procedure  is  given  in  Table  3-1. 

TABLE  3-1.  VECTOR  METHOD  WA  VEFO/  M  DESIGS 

I)  Maximize  s/c  ratio  given  scattering  properties 
W  Ky-j-'X)  W 

s/c  = 

W*  X  IX)  W 

VC 

2}  For  given  transmit  vector  X.  find  reat-vs»  weight  W  to  i.vaximize  s/c 
3  Constrained  optimization  prcSem,  eigenvalue  solution  of: 

,KvT  -  -'max  V  Wopt  =  0 

where 

-  maximum  eigenvalue  of  Kyy 

.  A)  W^-^eigenvector  of 

r* 

ir  51  Iterate  by  interchanging  X  and  Wr  obtain  convergence  to  optimam  s/c 
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One  additional  fact  should  be  noted  in  regard  to  receiver  noise.  The  effects  of  noise 
arc  taken  into  account  by  adding  an  identity  1  matrix  to  Kvc.  where  KyC  is  appropriately  scaled 
to  reflect  the  received  clutter-to-noise  ratio  (c/n). 


3.1.3  Generalized  Matrix  \pproach 

In  the  analysis  below.  we  consider  the  defection  of  a  point  Doppler-spread  target  in 
Doppler  spread  clutter,  when  the  radar  transmits  a  coherent  pulse  train  composed  of  N  equally 
spaced  rectangular  subpulses.  It  is  also  assumed  ‘hat  the  radar  is  capable  of  modulating  the 
transmitted  polarization  fro'  i  pulse- to-pulse.  A  simplified  block  diagram  of  me  system  is 
shown  in  Figure  3-3.  For  convenience.  the  time  notation  will  be  omitted. 

Given  3  transmitted  sequence  of  X  pulses,  the  transmil  vector  X  will  be  a  column 
vector  shown  as  a  (2S  X  i)  dimensional  matrix.  The  transpose  of  X  is  represented  by: 


X  -  «XV|.  X|||.  .  .  XVN.  XHN) 


<3-19) 


where  \yi  is  the  complex  variable  rcpresentating  the  amplitude  and  -'base  of  the  vertical 
component  of  the  pulse-!  transmitted  electric  field,  X^j  is  the  horizontal  component,  etc. 

The  scattering  properties  of  the  target  and  chaff  are  represented  by  their  respective 
1  x  2  scattering  matrices.  My  and  fly.  Given  that 


*11 

h12 

*21 

h:: 

—  J  T.C 


<3-:o) 


and  that  the  target  (T)  and  clutter  (O  scatterers  are  isotropic  (hj->  =  h->j).  then  the  elements 
of  the  scattering  matrix  can  be  rearranged  for  mathematical  convenience  as  a  scattering  vector 
h.  of  dimension  <3.\  X  !). 


—  (  « 1 1  < i >  S<Nl  ;.(N>  |X}\ 

h  =  y  h  j  j  ,  h]2  -  h->2 . ‘  |  j  *  *]->.}•■»->  J 


<3-:i> 


where  hjj  represents  the  vcrtical-to-vcrtical  scattering  proper’y  relative  to  pulse-1,  etc.  Note 
that  the  values  of  the  elements  varv  with  time  for  Doppler  spread  tarsets.  su'h  that  in  senera! 
n“>  *  h0>. 


The  signal  received  at  the  receiver  antenna  is  denoted  by  the  vector  V,  as  noted  iii 
figure  3-3.  For  mathematical  convenience,  the  transmit  vector  X  is  rearranged  into  a 
(2N  X  3N)  transmit  matrix  X,  such  that  Y  can  be  given  by 

Y  =  p  (3-22) 

where  Y  is  a  (2N  XI)  vector  with 

Y  =  <>VJ.  >'H1 . >VN*  >'HN*  (3-13) 

The  dements  of  Y  represent  the  complex  vertical  and  horizontal  components  of  the 
received  electric  field  vector.  For  Equation  (3-22)  to  be  dimensionally  vaiid. 


X  = 


XH1 

0  \ 

) 

X\T 

XH1  / 

(XV2 

XH2  0  \ 

\  0 

XV2  xHll 

(3-24) 


XYN  xilN  0  V 
0  XYN  XHN'/| 


The  above  rearrangement  of  terms  is  a  special  case  of  Equation  (2-6).  where 
hp  =  h-.j.  resulting  in  a  dimension  reduction  of  h. 

TARGET 
A -NO  CLUTTER 
MOOEL 

d<Dd06 

l,T  «  II  .  (2N  i  1> 


TRAJiSVlT 
VECTOR 
<2*  *  1) 


OBSERVATION 

NOISE 

!2N  t  1J 


Figi:re  3-3.  Simplified  Block  Diagram  Matrix  Method 
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lixpression  \3-24)  shows  the  matrix  ^  as  a  (2N  x  3N)  arrangement  of  the  elements 
of  the  transmit  vector  X.  Observation  noise  Nq  is  adoe*!  to  V  to  produce  the  received 
<2N  x  i)  vector  Y.  Tltis  realization  of  the  vector  Y  is  simply  the  result  of  the  transmit 
signal  X  Ming  reflected  back  to  the  radar  by  the  target  and  chaff  plus  the  additive 
receiver  noise. 

Ihe  preliminaries,  to  this  point,  of  rearranging  the  mathematical  form  of  the 
transmit  signal  and  scattering  properties  into  X  and  h.  was  necessary  to  properly  -Jimensionalize 
the  matrix  multiplication  for 

X  ~  X  h  +  Nq  (5-25) 

The  effect  of  a  linear  receiver  upon  the  signal  vector  Y.  can  then  be  expressed  in 
terms  of  the  matrix  G,  which  operates  on  Y  to  give 

Z  =  GY  =  G^ht  GN0  (2-26 } 

The  receiver  output  Z  would  then  be  subject  to  hypothesis  testing  for  one  c f  (wo  conditions: 
the  presence  of  target  plus  clutter  or  of  clatter  only.  The  implementation  of  this  likelihood 
process  will  be  deferred  to  Section  5.2.  A*  this  juncture,  the  itaior  consideration  is  the 
des«gfi  of  *.»  X  suc^  that  the  s.c  ratio  a!  Z  is  maximized. 

3-1-3. 1  Matrix  Transmtt-Rccehe  Design  -  Define  a  weighting  matrix  \\\  where 


W  =  G  X  (5-27) 

Ihe  object  then  is  to  find  the  optimum  W  for  maximizing  the  s. .  ratio,  pten  the  statistics  of 
the  stochastic  scattering  sectors  for  both  the  target  and  clutter. 

The  target  and  clutter  covariance  matrices  are  given  by 

*Sy  =  i:  [hy  hy*I  «5-2Sat 

and 


Kv.  =  i*»c  hc*i 


<3-2Sh> 


where 


*  represents  the  conjugate  transpose. 


Since  Kyj  and  Kyc  are  positive  definite  hermit  ian  matrices,  they  can  be  simultaneously 
diagonalized  by  a  matrix  p  (see  Appendix  D).  where 

O*  h.,T  O  —  A 


<>‘Kvcp  =1 

The  columns  of  O  are  the  generalized  eigenvectors  pj  satisfying 
Kyj  oj  =  X,  Kyc  p; 

where  the  Xj  are  the  eigenvalues  of  the  matrix  equation 


KvT  ~  X  Kvc5  =  0 


(3-29) 


( 3-30 ) 


(3-31 ) 


The  joint  diagonal  zati  on  of  Kvj  ~nd  Kv,c  by  this  process,  orthogonaiiz.es  the  target 
and  chaff  covariance  matrices  through  the  matrtx  "transformation  p.  This  in  effect  provides 
a  maximum  separation  of  the  target  and  clutter  vectors  at  the  Z  output,  which  is  the 
desirable  condition  for  hypothesis  testing.  Consequently,  the  matrix  \V  can  be  chosen  as  the 
transformation  matrix  P.  There  is.  however,  a  dimensionality'  problem  in  that  p  is  a 
(3N  x  3N)  matrix,  the  result  of  the  arrangement  of  the  h.  On  the  other  hand,  there  are 
only  2N  observations  at  Z  as  iimi‘ed  by  the  transmission  of  N  dual  channel  pulses. 
Consequently,  the  dimension  of  the  W  matrix  must  be  reduced  (k  x  3\T.  where  k  2N. 

The  problem  then  arises  as  to  the  method  of  selecting  k  out  of  3N  eigenvectors  for  the 
matrix  \V,  where 


W  =  (W.  W, 


(3-32) 


and  the  W-  are  <3N  x  1)  column  vectors,  representing  the  k  eigenvectors  of  the  transformation 
matrix  p.  that  correspond  to  the  k  larges!  eigenvalues  of  matrix  equation  (3-31). 

The  criterion  for  selecting  the  x -optimum  eigenvectors  is  based  upon  the  Chernon 
bound,  whereby  the  probability  of  error  is  ..linimized  with  respect  to  the  hypothesis  test  on 
Z.  Given  Hypothesis  Hq  (clutter  only)  and  Hj  (target  plus  clutter),  i:  car.  be  shown  th-at  the 
probability  of  error  Pc  is  bounded  between  JVan  Trees.  19'*IJ. 


lisnin  *p.a)  s~  <  P„  <  li  J 


(3-33) 


where  p  is  the  a  priori  probability  that  Hj  is  true,  q  is  the  a  priori  probability  th2t  Hq  is 
true  and  I  is  the  integral 


!P(Z.Hj)  pfZ/Ho)]'1  dz 


(3-34) 


For  zero  mean  Gaussian  densities  (see  Appendix  E), 


fc/i  *  Kzo:  I 

B  =  -?n(J)  =  Cn  <  .  w, 

-  lkZll  '  >kZ0* 


(3-35) 


where  B  =  -Cn(J)  is  known  as  the  Bhattacharyva  distance.  N<  te  that  since  0  <  J  <  1.  the 
minimization  of  J  is  equivalent  to  the  maximization  of  B.  Thenelore.  we  select  the 
transformation  matrix  W  so  as  to  maximize  the  Bhattacharyva  distance. 

in  Equation  (3-35).  KZ1  and  K z0  are  the  covariance  matrices  of  Z  under 
hvpotheses  Hj  and  I Iq  respect ivelv .  These  quantities  are. 


KZ1  =  W  KvT  W*  +  w  Kvc 


KZ0  =  W  Kvc  W 


(3-36) 


Substituting  ot  the  W  matrix  as  resulting  from  the  eigenvector  solution  of 
Equation  (3-32)  into  the  above  for  KZ1  and  Kz0.  and  then  substituting  of  KZ1  and  Kzo 
into  the  Bhattacharyva  distance  B  as  given  in  Equation  (3-35)  yields 


il£n  (Xj  +  1)  +  r—TT  +  -  I  +  k'2  kn  a 

i=l  |  ‘  J 


(.so  ) 


Note  that  due  to  the  restriction  or.  k.  k  must  satisfy  <2N.  Then,  maximization 
B  requires  finding  the  appropriate  k  eigenvalues  a- . Xj.  out  of  3N  eigemalues  that 


satisties 


-  1>  *  <X,  *  If1;  >  {«> -  -  ! i  -  <X,  *  lf’*i  >...>!•«♦  i*  *  *  if 

*  I  i  *  *■  *“ 


>  ;>  *  <>5X  *  if’ 


implies 


tX  -  ! i  -  <X  -  if1  >«.\  -  ii  *  ♦  if1 


.\  »  for  ail  X.  X  >  o 

i  j  «  j 


mi  .-:  eigenvalues  that  maximize  B  an'  then  the  k  largest  eigenvalues  of  the  ordered  set 


,Xma\  =  V  >  X2  -  -  >  H  >  -  -  >  X3N 


1 3-3f>  t 


It  then  follows  that  tin  k  eigenvectors  of  W  should  be  selected  from  tho^c  which  correspond 
to  he  k  largest  eigenvalues  of  -  Xkv-.  =  0 


’  •v,£^ si**. 


3. 1.3.2  Transmit  Waveform  Determination  -  Once  the  aoove  optimum  matrix  W  has  been 
determined  the  receiver  transformation  matrix  G  and  transmitter  waveform  matrix  X  arc- 
selected  bv  minimizing  the  Frobenius  norm 


8  G  X  -  W  5  F 


(3-39) 


which  is  defined  by 


“ai!f-  15  'u 


<V40: 
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A 
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The  choice  cf  the  Frobenius  norm  here  is  for  mathematical  convenience,  since  in 
finite  dimensional  spaces  all  norms  are  equivalent  -  as  far  as  topological  properties  are 
concerned.  For  a  given  waveform  matrix  X.  the  minimizing  G  is  given  tv 


G  =  W  X*  (X  X*)_! 


(3-41 ) 


Substituting  Equation  (3-41)  into  the  expression  for  the  norm,  yields 
1  W  (XMXX*)"1^-  !*  1  }• 


(3-42) 


J? 
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Optimum  performance  is  obtained  if  liquation  (3-42)  is  identically  zero.  The  above 
discussion  indicates  that  for  a  given  matrix  W  the  transmit  matrix  X  can  in  theory  be  computed 
by  minimizing  the  Frobenius,  norm.  In  reality,  this  direct  method  for  solution  of  X 
extremely  difficult.  An  alternative  approach  is  to  start  with  a  large  library  of  reasonable 
waveform  structures,  and  search  for  the  waveform  which  minimizes  the  norm.  Thi;.  was  done 
during  the  study  and  will  be  reported  in  Sections  4  and  5. 

Once  X  is  selected  through  this  procedure,  the  receiver  matrix  G  is  computed  from 
Equation  (3-41).  A  summary  of  the  matrix  approach  is  given  in  Tabic  3-H. 

TABLE  3-1L  SUMMARY  MATRIX  WAVEFORM  APPROACH 

1)  Solve  (or  eigenvalues  of  IKyj  -  *  =  0 

2)  Select  the  Sc<2N  largest  eigenvalues 

3  rovn  'i¥  matrix  from  the  '«  e  gemrectors  of  largest  eigenvalues 

4)  Solse  for  transmit  matrix  by  minimizing  norm  of 

riv  (X*  (£5p_1  X-IJ12 

5)  Step  4  a^ove  can  be  satisfied  by  comparing  a  library  of  stored  transmit  vectors  against  W,  and 
selecting  the  one  resulting  in  the  minimum  error 

«?  65  Solve  for  receiver  weight  G  using 

S  .  -  W./IXXT1 


3-  *  -4  Maximum  Likelihood  Detector 

The  matrix  approach  yields  a  transmit  waveform  with  an  associated  receiver  weight 
which  maximizes  the  s/c  ratio  based  upon  the  Bhattacharyya  distance  criterion  for  optimalitv. 
Up  to  this  point  of  the  design  (output  Z  of  Figure  3-3),  the  hypothesis  test  for  the  target  or 
chaff  decision  process  has  not  been  a  factor.  For  a  surveillance  radar  application,  one  is 
interested  in  optimizing  the  probability  of  detecting  a  target  ir.  chaff  for  a  given  false  alarm 
probability  created  by  chaff  plus  receiver  noise.  Referring  to  Figure  3-3.  the  two  conditions 
for  hypothesis  testing  are 

H0:  Y  =  Yc  +  N  o 

Hj :  Y  =  Yj  -r  Y;-  -r  N0  1 3-43 > 

Here  Hq  represents  the  false  alarm  situation  of  chaff  return  plus  noise,  while  Hj  represents 
the  target  return  combined  with  the  chaff  and  noise. 

The  above  receiver  inputs  are  lineariy  operated  upon  by  the  weighting  matrix  («  to 
produce  the  observed  receiver  output  Z.  A  decision  must  then  be  made  as  to  whether  Z 
satisfies  hypothesis  Hq  or  Hj.  Tne  transmit  vector  X  and  receiver  matrix  G  are  now  known 

design  parameters.  The  statistics  of  the  target  and  chaff  scattering  properties  as  well  as  the 

receiver  noise  are  known  stochastic  processes;  consequently,  the  signal  distribution  at  Z  is 
also  known  for  each  hypothesis.  This  then  lends  itself  to  a  conventional  maximum  likelihood 
decision  process  for  Hq  and  Hj.  The  maximum  likelihood  receiver  configuration  can  *cund 
in  standard  texts  (Van  Trees,  1968J,  and  will  not  be  elaborated  upon.  There  are.  however, 
some  design  features  peculiar  to  the  matrix  approach.  These  arc  discussed  below  with  the 
aid  of  Figure  3-4,  which  is  an  extension  of  Figure  3-3. 

At  the  threshold  T,  a  binary  decision  is  made  for  condition  Hq  and  Hj.  Since  it  is 
desired  to  utilize  a  conventional  receiver  implementation  such  as.  for  instance,  a  square  law 
detector,  it  is  necessary  to  transform  the  output  at  Z  to  fulfill  certain  conditions.  This  is 
indicated  by  Figure  3-4  by  the  matrices  0  and  P  which  operate  upon  Z  to  produce  the 
output  signal  r.  This  signal  is  then  squared  and  summed  to  provide  a  signal  u,  upon  which 
the  threshold  test  is  performed.  The  mathematics  associated  with  the  decision  process  involves 
knowledge  of  the  probability  density  function  of  u,  and  this  computation  can  be  facilitaicd  if 
the  elements  of  the  signal  matrix  at  r  are  uncorrelatcd.  This  is  accomplished  by  finding  a 
transformation  matrix  Q  which  diagonalizes  the  covariance  of  signal  Z;  that  is. 

Z  =  QZ  (3-44) 


where 


E|Z  *  i  =  diagonal  matrix. 
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!l  is  noted  that  the  covariance  of  Z  may  itself  satisfy  the  orthogonality  criterion. 
This  would  be  true  if  in  the  matrix  waveform  design,  the  quantity  GX  were  identically 
equal  to  the  matrix  \V.  In  general,  however,  the  Frobenius  non  i  ?  GX--V  h  ^’ssfied  by 
a  least  squares  fit  procedure.  Mich  that  GX  =4  W.  In  »his  cast.  Z  must  be  orthogonalized 
by  tiic  transformation  matrix  Q.  Additive  receiver  noise  also  requires  imposition  of  the 
Q  matrix. 

The  matrix  P  which  operates  on  Z  to  produce  r.  is  simply  a  >cale  factor  to 
properly  weight  the  elements  of  r  prior  to  the  square  law  operation.  The  final  observed 
receiver  output  u.  is  then 


r*  r  =£  1  fj  j- 


(3-45) 


The  likelihood  ratio  hypothesis  test  is  dependent  upon  the  probability  density  of  u 
under  the  criteria  Hq  and  Hj.  The  probability  of  false  alarm  for  a  given  threshold  T,  is  then 


Pfa  "  L  P(ui,J0)du 


(3-46) 


For  a  desired  Pp^,  the  threshold  T  can  be  computed  using  the  Newton-Raphson 
tcchtiique.  Given  T,  the  probability  of  detection  then  becomes 


J 

1 

I 


w 


PD  =  fj  PfolHjJdu 


(3-47) 


A  target  is  declared  if  u  >  T.  Note  that  the  optimum  receiver  overall  weighting 
function  is  sow 


Gopt  -  P-O-G 


(-MS) 


Appendix  8  provides  a  detailed  derivation  of  the  foregoing  procedure,  assuming 
that  Z  is  jointly  Gaussian  under  both  hypotheses  Hq  and  Hj.  It  should  be  noted  that  the 
maximum  likelihood  detector  decision  process  can  be  applied  to  any  of  the  waveform  design 
techniques  considered  during  this  study. 


3.1.5  Preliminary  Design  of  the  Processor 


The  processing  algorithm  is  presented  in  a  functional  flow  form  in  figure  3-5.  A 
sample  vector  (k  x  1)  is  multiplied  by  a  we  shting  matrix  (k  x  k).  The  elements  of  both  the 
sample  vector  and  the  weighting  matrix  are  <.  implex.  This  produces  a  (k  x  1)  result  vector  r 
with  complex  elements.  The  elements  of  the  result  vector  r  are  squared  and  the  squares  are 
summed  to  form  a  scalar  u  which  is  then  compared  with  a  threshold  T. 


For  the  preliminary  design,  k  is  assumed  equal  to  6,  Under  this  assumption,  the 
sample  vector  will  contain  six  complex  (I.Q)  sample-  one  sample  from  a  channel  with 
horizontal  (H>  polarization  and  one  sa^-Ic  -fern  a  channel  with  vertical  (V)  polarization  fer 
each  of  three  pulses.  A  pr be  repetition  lime  of  three  milliseconds  is  assumed.  Further,  the 
quantizer  is  turned  to  quantize  the  1,Q  quadrature  components  to  10-bits. 


Under  these  assumptions,  the  processing  load  can  be  estimated  as  follows: 
Processing  Per  Sample  Vector 


Matrix  Multiplication: 

Vector  Squaring  and  Summing: 
Tnreshold  Detection 


144  real  multipiys 
•32  adds  subtract^ 


12  real  multipiys 
1  1  adds 


«  sub?  rail 


These  data  may  be  summarized  as  follows:  each  three  millisecond  period,  the  processor  is 
required  to  perform  156  multiplications  and  144  add /subtracts  for  each  sample  vector 
processed.  Input-output,  bookkeeping.  3rd  ether  overhead  operations  are  in  addition  to  these 
requirements. 


If  these  operations  were  tc  be  pcTormed  in  a  general  purpose  16-bit  microprocessor 
or  computer  without  any  specialized  multiplication  hardware,  a  reasonable  assumption  would 
be  one  16-bit  by  16-bit  multiplication  each  16-clocks,  ar.d  one  add  or  subtract  operation  each 
clock.  Using  these  assumptions  the  minimum  processing  load  u-oula  be  (156  x  16)  + 

1**4  —  2640  clocks.  If  In  addition,  we  estimate  four  clocks  of  overhead  for  each  element  of 
the  sample  vector  and  of  the  weighting  matrix,  we  require  an  additional  4  x  <12*  +  72)  «=  336 
docks  for  a  total  of  2976  clocks.  This  can  be  rounded  up  to  3,000  clocks. 


figure  3-5  contains  a  processing  block  diagram  which  contains  a  1 6-bit  microprocessor 
identified  as  the  ITT  Giifiikm  Radar  Processing  Module  (RPM).  The  RP.M  is  capable  of 
operating  at  a  clock  speed  of  250  nanoseconds.  A  first  estimate  will  be  made  assuming  that 
the  entire  -ask  were  performed  in  the  RPM.  This  requires  3000  x  0.250  =  750  microseconds 
during  each  three  millisecond  period  for  processing  of  the  target  in  a  single  range  bin. 
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figure  3-5  shows  the  use  of  3  specialized  programmable  Outboard  Processor  designed 
to  unload  the  RPM  tG  improve  processor  efficiency.  The  detailed  block  diagram  of  the 
Outboard  Processor  is  given  in  figure  3-6  and  represents  an  alternative  to  RPM  only 
processing,  noting  that  the  RPM  is  still  required  for  data  management  and  other  ranging 
functions.  This  processor  is  built  around  an  available  12-fcit  by  !  2-bit  LSI  multiplier 
accumulator  chip.  This  processor  is  designed  to  receive  quadrature  samples  and  matrix 
weights  from  the  RPM  in  one  set  of  memories  while  it  is  processing  data  from  the  other 
set  of  memories.  This  processor  is  designed  to  perform  the  matrix  multiplication  and  vector 
squaring  operations.  Addition. subtractions  into  the  accumulator  are  fully  overlapped  with 
multiplications.  A  total  of  26  multiple  and  accumulates  are  required  to  produce  the 
squares  of  the  real  and  the  complex  parts  of  the  elements  of  the  result  vector  r.  >  5  for 
the  real  part  and  13  for  the  imaginary  part.  This  requires  a  total  of  6  x  26  =  156  docks 
for  each  sample  vector  processed.  The  elements  of  the  matrix  weights  are  changed  only  as 
required.  In  this  configuration,  assuming  a  250  namsecond  clock,  a  matrix  mu’i.ply  and 
vector  square-  operation  can  be  done  in  156  x  0.25G  -  3‘»  microseconds  (muen  more 
effident  than  the  750  microseconds  associated  with  the  RPM  for  this  function). 

With  reference  to  Figure  3-5.  the  RPM  designates  the  range  at  which  samples  are  to 
be  collected  to  the  Range  Gate  Logic  which  strobes  the  outputs  of  the  H  channel  and  the 
V  channel  A/D  converters  and  passes  the  samples  to  the  RPM.  The  RPM  passes  the 
samples  and  the  matrix  weights  to  the  outboard  pro^essc.  which  computes  the  squares  of 
the  elements  of  the  result  vector  rand  passes  these  squares  back  to  the  RPM  where  the 
squaies  are  summed  and  threshclded. 

The  configuration  of  Figure  3-5  would  be  effective  for  implementing  a  track  function 
in  the  radar,  resident  in  the  RPM.  To  apply  the  same  processing  to  an  acquisition  function 
in  the  radar,  additiora!  processing  capability  would  be  required.  Adaptation  of  the 
ITT  Gilfillan  developed  Optimum  Filter  Transform  (OFT)  modular  approach  would  provide 
the  required  additional  processing  capability  (the  OIT  module  implements  a  vector 
times  a  matrix  multiplication  and  was  developed  for  high  performance  Doppler  processing 
in  radar  signal  processing  applications). 


3.2  IKE  STATE  SPACE  APPROACH 

Tiiis  section  is  also  devoted  to  modeling  of  radar  targets  and  clutter  and  to  the 
design  of  an  optimum  receiver  for  detecting  a  target  in  clutter.  Here,  as  in  previous  sections, 
both  target  and  chaff  are  taken  to  have  a  dual  channel  scattering  operator  description,  but 
our  approach  will  be  that  of  the  modem  state  space  theory  of  sysiems. 

From  an  external  viewpoint,  a  system  can  be  described  by  the  family  of  ail  possible 
input-output  pairs.  This  may  be  the  most  general  description  of  a  system  but.  unfortunately, 
it  does  not  give  a  useful  relationship  between  a  given  input  of  a  system  and  its  corresponding 
output.  Because  of  this,  we  introduce  the  state  space  theory  of  systems. 

The  state  of  a  system  is  the  least  amount  of  information  which,  together  with  the 
input,  will  uniquely  determine  the  output.  In  what  follows  we  shall  be  mostly  concerned 
with  discrete-time  systems  nith  a  finite  dimensional  state  space.  Thus,  let  t  be  the  time 

variable  which  takes  the  values  0.  1.2 . and  let  A.  B,  C  be  matrices  of  appropriate 

dimensions.  Then  a  linear  discrete-time  system  is  described  by  the  pair  of  equations 


sf!  +  1)  =  A  sft)  +  B  u(t ) 


►41 


qui  =  C  s(t) 


where  s(t)  is  a  vector  called  the  state  st  time  t.  while  u(t)  and  qtn  are  the  input  vector  and 
the  output  vector,  respectively. 

3.2.1  State  Space  Modeling  of  Radar  Targets  and  Qmter 

3.2.1. 1  Point  Target  Model  -  Tire  state  space  theory  is  not  onlv  developed  for  studying  sys¬ 
tems,  but  it  is  aiso  intended  for  ‘‘generating”  stochastic  processes  -  by  means  of  systems 

for  applications,  for  instance,  in  communication  and  control  systems.  It  is  evident  from 
Equation  (3-49)  that  if  the  input  u(t)  is  a  stochastic  process  then  so  is  q{t).  Therefore  if  a 
stochastic  process  can  be  generated  from  a  dynamic  system  with  a  state  space  description, 
then  we  have  3  state  space  model  for  the  process 

We  now  derive  a  stale  space  model  for  a  point  target.  Let  Ky(t.  r)  be  the 
scattering  matrix  of  the  target,  then 

Hj(t,  7)  =  H(t)  •  5 ij  -  r0).  H(ti  =  fhgft)  J  .  i.  j  =  1.2.  (3-50) 

where  hy{t)  are  Gaussian  processes.  If  x(t)  is  a  transmitted  signal  -  which  is  a  2  bv  ? 
deterministic  vector  -  and 


S(t)  -  l  hj  j(t».  hp(!).  hij{0,  hii(t)  J 


(3-51) 
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ts  the  transpose  of  the  target  scattering  vector,  then  the  corresponding  received  signal  is 


Y(t)  =  2^<t  -  fQ)  •  S(t) 


(3-52) 


where  the  2  by  4  matrix  X(*)  is  given  by 


-  Tq)  - 


\j(t  —  Tq;  Xstt  ~  "q) 
0  0 


\j(t  -  Tq)  >.««  -  Tq)  J 


I  :>->:>) 


which  is  deterministic.  It  fellows  from  Equation  (3-52)  that  the  receiv  -J  signal  vector  Ym 
can  be  regarded  as  the  output  of  a  linear  time-varying  dynamic  sv  stern  whose  state  vectar  is 
the  scattering  vector  Sri).  Tnus  the  problem  of  modeling  the  point  target  becomes  that  of 
modeling  the  vector  process  S(l). 


To  model  the  process  S(t)  we  use  an  Gutore%re$she  model;  that  is.  Sri)  is  taken  to 
be  of  the  form: 


N 

Sri)  =  A,,  Sti  -  n)  +  u(t). 

n=l 


( _>-54) 


where  the  4x4  constant  matrices  An  aic  to  be  computed  and  the  4  \  1  vector  r-u>  «s  z 
white  noise  process  with  zero  mean.  The  N  step  autoregressive  model  Equation  «3-54)  is 
actually  a  form  of  a  state  space  model  and  is  equivalent  to  the  pair  of  equations: 


Sri)  =  A  *  Sri  -  I)  +  U(t). 
Sit)  =  C  -  Sri). 


(3-55) 


where  St  -)  and  U(*>  are  4  \  i  vectors,  while  A  and  C  are  constant  tr-iirices  of  dimension- 
4N  x  4N.  and  4  \  4N  respectively. 


The  basic  problem  in  representing  S(t)  by  an  autoregressive  model  is  to  estimate 

the  matrix  coefficients  Aj.  A-s . A^-.  from  the  observations  S(l).  St2- . ?»N).  I.ct 

us  now  illustrate  the  estimation  procedure  for  3  one  step  autoregressive  model.  Thus  let 
Sri)  be  generated  by  the  equation 


Sit)  =  Aj  •  SO  -  1)  +  u(t). 


(3-56) 


where  Aj  =  (aj-j.  i.  j  =  1.  2,  3.  4.  is  tc  be  estimated.  For  this  we  define 


a  -  1  3j».  aj2«  aj3-  3j4-  a2i-  -  -  -  a24»  a3l*  -  -  a34-  a4j  -  -  -  a44  i-  ^ 3-571 


Then  Equation  (5-56/  can  be  written  as 


S(t)  =  Q  (t  -  1)  -  a  u(t). 


(3-58) 


where  the  4x1'  matrix  p  is  made  up  of  entries  of  the  4  x  !  vector  S(t  -  1).  Now  if  the 
values  S(l)  .  .  S.L)  where  (L>5)  are  known,  then  we  form  the  4L  x  1  vector 


SC.  .  .  .  L)  =  JSC).  .  .  St  L>1 

and  the  4L  x  *6  matrix 

p(l . L  -1  ‘  =  |p(  ] ) 

1  ! 

|p<L  -1) 

It  then  follows  from  Equations  (3-58).  (3-59)  and  *3-60)  that 

SC.  .  .  .  L)  =  o  (1.  .  .  .  L  -i)*  a  +  UC.  .  .  .  L» 


(3-59  > 


(3-60) 


(3-61) 


where 


UC.  .  .  .  L)  =  JUC).  .  .  .  L<L)1 


To  estimate  the  vector  a  -  which  in  tum  gives  the  estimate  of  the  matrix  Aj  -  we  use  the 
!east  square  criterion,  i.e..  w*»  minimize  the  functional 

1  SC.  ...  L *  -  £  tl ....  L-l )  -  a  1-.  (3_5: 

A 

Then  the  least  square  estimate  a  of  a  is  given  by 


A  -1 

-  =  !  1 .  ...  L-l)  -p  tl.  ...  L-l  i]  *  pM.  .  .  .  L-l)  -  SC.  .  .  .  L) 


(3-63) 


3.2.1. 1.1  Computation  of  Matrix  Coefficients  -  In  the  case  of  our  problem,  the  observations 

are  actually  the  transmitted  signal  ri3ta  set  [X(!) . X(N)j  and  the  received  signal  data 

set  [Y(l> . YtN)  j.  I*  is  evident  from  Equation(3-52)  that 


where 


S(t)  =  |  #  (t  -  r0)  J*1  •  X(t  -  7q)  -  Y(t) 


= 


XjXj  X,X2 

X2X1  X2X2 


0 

0 


0  1 
0 


XjX,  xjx2 
X2Xj  X->X-I 


=  X*  -  X 


(3-64)- 


(3-65) 


Therefore  the  data  set  (S(l),  S(2) . S(N)J  can  be  computed  from  Equation  (3-64),  using 

the  transmitted  and  received  data  sets. 
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The  matrix  wiildonts  Aj.  A->.  .  ..  Aj^.  can  also  be  computed  from  the 

Yuie-Waiker  equations 

AjKcim  -  t  AsXg£;*>  -  2t  ‘  .  .  +  A^Kg  (m  -  N)  =  Kg  (m). 

‘or  n:  =  !.  2.  3 . N. 

Here  Kg  is  the  covariance  iirasfj*  of  the  process  Si*.:  ‘‘■hid:  can  K  approximated  by 


(3-66) 


N-* 


Kgi'K)  -  >  Stn  +  k)  *  Sm**.  fo 


-  0.  i.  2.  . 


(3-67) 


Hquivalcntiy.  one  can  also  compute  K^t*)  from  the  :vrdve-y  cenal  covariance  matrix  Kyl*l. 
as  discussed  in  Section  2.  : 

Kgtt  -  w',  =  (X  <t  -  r0)J-'  -  XC  -  rQf  KyM  -  w>  -  Xiw  -  rQt  * 

I3thv.  -  -0»)“5  (3-68) 

3.2.1. 2  Ciatf  Model  -  In  ths.  raj*-  **c  ha.c  :he  t s cnsrr.i  1 1  cd -'e re  i « c-d  Signal  relationship 


Yin  =  j_oo  ~  4)  *  S::.  r;  a*. 


(3-6Q) 


where,  as  in  the  abo.e  X(- 1  is  the  transmitted  signal  matrix  and  St*.  *s  the  range-spread 
chaff  scattering  vector.  Tfn.  problem  row  hcct-P'es  that  of  modeling  the  input  Sit.  r> 
which  is  a  stochastic  process  m  t  >'T  cadi  fixed  r  -  of  a  tinea:  s\  stem,  given  the  sy  stem 
weighting  matrix  \<  •)  and  the  outpu»  y  We  chan  now  gjv,.  an  approximation  procedure 
for  this  case. 

Let  I6n  (t)  1  =°n=  |  be  a  sequence  of  crthcnoimal  4  x  3  sectors  over  the  interval 
(- °o.  =»).  i.e.. 


|  i  for  n 

/J^  l‘4n  <r).  ym  (r)|  di  -  2>n  m  -  j  0  for  n 


m 

m  • 


l3-~0J 


where  J...J  denotes  the  inner  product  in  C*.  i.xpanding  the  4  x  1  vector  S(l.  r)  in  terms  ol 
the  sequence  i  y»n  (r)l  wc  find 


S(t.  r)  -  ^  ^  sn  (l)  y’n 


(7) 


(3-~l » 
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where  the  scalar  functions  s^I)  are  random  processes.  Let  Kgt.t,  r;  w.  o)  be  the  covariance 
matrix  of  the  process  S(!  r). 

Ks(t.  r:  w.  a)  =  ^  £  Anm  «-  w>  '  <r>  *  v'm  <o)*  (3'72) 

n  m 


where 

Anm  (t-  w>  =  E  lsn  (I>  *  sm 


(3-73) 


Therefore 


oc  oo 

Aj..  (t.  w)  =  /  /  (T>*  -  Kg  (t.  t:  w.  0}  -  oj.  fci  -  dodr. 

j.  k  =  1.  2.  3 .  (3-74) 

Having  computed  Aj«-  (...)  from  the  covariance  matrix  Kg  (...;  ...)  we  can  then  mode!  each  of 
the  processes  sn(t)  by  the  autoregressive  scheme  discussed  above.  Equivj'ent!).  we  can  mode! 
the  vector  process 

T.(i>  =  lsj(t)  s^O)  .  .  .  Sj/t))  (3-75) 

with 

E  IsJ,(t)  *  s^(w»*}  =  Matrix  {Ajj.(t.  w)J  (3-76) 


by  the  same  procedure. 

it  is  noted  that  the  received  signal  Y(t),  in  this  case,  admits  the  mean 
square  approximation 


N 


v  *  £y)  •  vv» 


(3-77) 


where  the  2  x  1  vectors  Vn(t)  are  given  by 


Vn(t)  =  /_*«,  X(t  -  t)  *  v>n(r)  •  dr.  n  =  1.  2.  3 . N. 


(3-7S) 
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3.2.2 


Optimum  Receiver  Design:  The  Case  of  a  Point  Target  in  a  Cloud  of  Chaff  Dipoles 
Let  Y(t)  he  the  received  signal  from  a  point  taigel  in  a  cloud  of  chaff  dipoles,  then 

Y«»)  =  Yjtti  +  Yc<ti.  <3.-9) 


where  y  and  ’l  r  are  the  received  signals  from  the  target  and  .’  .half,  respective);  .  and 
both  ^  -j-  and  'l  y  are  Gaussian  processes. 

The  binary  detection  problem  is  to  find  the  option  m  receiver  for  testing 


IV  Null  Hypothesis  *i.e..  Zero  Target  Signal): 

Yit)  =  Y^*i)  -  m).  for  .ili  t  in  jTj.  T>j;  i3-s0> 

Hi:  Alternate  Hypothesis: 

N 1 1 *  =  Y-jit)  *  Y<-U >  +  l-l t».  for  a'!  t  in  jTj.  Ts).  <3-Sl) 

In  the  above,  ii  t »  i>  an  external  additive,  statistically  independent  2  i  vector,  white  noise 
process.  T]  and  T s  an.  the  initial  and  final  observ.tion  times,  and  the  interval  | T j .  T^j  is 
assumed  to  he  long  enough  for  complete  observation  of  the  received  signal. 

3-2.2. 1  State  Space  Mode!  for  the  Point  Target  -  Without  lost  of  generality,  the  point 
target  is  modeled  bv  a  "2-step"  autoregressive  scheme.  Tlius  lit] nation  (3  54  >  now  becomes 


Sui  =  A.  *  Stt  -  it  +  As  •  St T 


2)  +  u<  J ) 


i 3-82 ) 


where  the  4  \  4  matnces  Aj  and  A-  have  been  computed  from  the  metr.ods  jiscusseu  above 
and  utt)  is  a  star  isti-caliy  independent  4  \  1  vector  white  noise  process.  Set 


Stt) 


Sit) 

Sit  -  lr 


■  v  O'*. 

» e-8.' ) 


Then  the  state  space  model  for  the  point  target  is 

S(t  -  1)  *  utt)  (3-34. 


y pi  t  >  = 

Xt:  -  r0)  *  (14  041  *  S-n 

— 

r  »<n  i 

where  urn 

=  ,-°4  J 

and  G_j  i 

is  a  4  v  1  zero  vector. 
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Stt  r  = 


A,  A: 


lU 


Vlt)  = 


fi'i  •  vlt) 


0 


4.1 


The  main  information  given  for  the  design  of  the  optimum  receiver  is  the  target 
received  signal  covariance  matrix  Ky  (t.  w),  the  chaff  dipoles  received  signal  covariance 

matrix  (t.  w),  and  the  white  noise  process  tit)  covariance  matrix  tt  -  w)  = 

K  *  5  (t  -  w). 


3.2.2-3  Estimation  of  Vj<{)  under  Hypothesis  Hj  -  In  the  state  space  target  .snu  chan 
formulations,  set 


AT 


.  and 


Vh  ' 


(3-93) 


Then  define 


(3-94) 


We  obtain  -  by  combining  Equations  (3-80),  <3-81).  (3-93).  and  (3-92)  -  the  following  state 
space  formulation  for  the  received  signal  Yft)  under  hvpothesis  H  | : 


Z(t) 


-  Z(t  -  1)  +  Wtti 


(3-95) 


Y(t)  =  Cj(t)  •  Z(t)  +  tit) 


where 


Cj(t)  =  (Xu  -  n.  0-,  pi  (0-,  §  /  X(t  -  r)  -  NcnI;dt.  4)  (3-96) 

=  r-p(t)  C^t)  (say) 

It  then  follows  that  the  problem  of  estimating  Y-j-(t)  now  becomes  that  of  estimating  the 
state  vector  Z(t).  This  is  discussed  below. 

3.1.2.4  The  Kalman  Filter  for  Estimating  Zft)  —  The  least  square  estimate  Z(t)  of  Z(t)  - 

basing  on  the  observations  Y(t).  Y(t  -  1) . Yi.O)  -  Is  obtained  by  the  discrete  Kalman 

filter.  The  derivation  of  the  Kalman  filter  is  standard  and  will  be  omitted  here.  Figure  3-7 
gives  a  block  diagram  of  the  filter. 
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Figure  3-7.  The  Kalman  Filter  for  Estimating  Z 


The  basic  equations  for  the  realization  of  the  Kalman  filter  are: 

t\ 

a)  The  estimate  at  time  t  of  the  state  of  Equation  (3-95)  is  Zttlt  -  1)  - 
conditional  mean  -  and  it  satisfies  the  recursive  equation: 

a  A  A 

Z(t|t  -  i)  =  A  *  Z(t  -  lit  -  2)  +  Gtt)  | Yin  -  Cj(t)  *  Zft  -  lit  - 

A 

with  the  initial  value  ZfliO)  (say)  given.  The  matrix  Gm  is  giver,  by 
GC)  =  A  *  P(t)  •  Cjfr)  |C,(\)  •  PG)  •  C,*(t)  +  K,,!"1. 

In  the  above 


A 


AT 

°4 


°4 

*Ch 


(Equation  3-95) 


and  1C,  is  the  covariance  matrix  of  the  white  noise  process  i>(t). 


b)  In  Equation  (3-98)  the  matrix  P(t)  is  the  covariance  of  the  estimation  error 
and  it  is  given  by 


P(t  +  1)  =  (A  -  G(t)  •  C,(!)J  Pit)  IA  -  G(t)  *  C  it)j  +  K  (3-HX 

1  ' 


+  G(t)  •  K  •  G*(t) 


^,th  the  initial  value  P(0)  taken  to  be  the  covariance  matrix  of  the  >nit-al 
state  ZiO). 


^•-•2-3  Estimation  of  VV-tt)  Imd?r  HVDO*hpc»c  _  Tu,  -l- ,  t-r  |  i  . 

_ -  v-  '  — :f_2_ ~  Q  1,,c  c,oud  received  signal  V^-ni 


can  of  course  be  estimated  by  a  Kalman  filter  as  in  the  case  of  esiimatinc  the  tarset  receiver 

signal  Y-j-tt).  However,  it  is  clear  that  once  the  shite  vector  7ms  has  neen  estimated,  the 
e»iin>ation  of  is  readiiy  ghcn. 


3. 2. 2.6  pie  Optimum  Receiver  -  The  optimum  receive!  is  described  in  Figure  3-8.  It 
compares  Ypt)  under  31,  with  Y(t)  -  Y^U  -  under  H0  against  the  decision  threshold. 


orcisio?. 

threshold 


Figure  M  The  Optimum  Receiver 


3.2.3  State  Space  Summary 

The  foregoing  approach  has  provided  a  state  space  model  for  the  design  of  an 
optimum  receiver  in  terms  of  detecting  a  point  target  in  ciutter.  The  receiver  implementation 
is  giver,  in  the  conventional  Kalman  filter  realization,  based  upon  knowledge  of  the  target/ 
clutter  scattering  properties  and  a  given  transmit  waveform.  Further  study  will  be  necessary 
to  also  utilize  the  state  space  model  in  the  design  of  optimal  transmit  waveforms. 

3.3  THE  FREDHOLM  INTEGRAL  EQUATION  APPROACH 

We  discuss  in  this  section  the  design  of  yet  another  form  of  -eceiver  for  the  bit  ary 
detection  problem  of  the  previous  section.  It  consists  of: 

a)  An  optimum  linear  filter  Eq  which  estimates  the  target  received  signal  Y-pt) 
from  an  observation  Y{t)  -  assuming  that  hypothesis  Hj  is  true,  and 

b)  A  prewhitening  filter  Qq  for  the  chaff  cloud  interference. 

This  receiver  a  shown  in  figure  3-9,  where  the  two  filters  Eq  and  Qq  are  solutions  of  tw'o 
matrix  Fredholm  integral  Equations. 

We  now  illustrate  this  approach  for  the  case  cf  a  point  target  in  a  cloud  of  chaff 
dipoles  considered  in  the  previous  section. 


at 


DECISION 

threshold 


Figure  3-9.  The  Correlator- Delector  Receh-er 


3-31 


3.3.1  integral  Equation  for  the  Filter  Qq 

The  dual  channel  fitter  has  the  weighting  matrix 

Qo(t,  u)  =  K'1  \\2  *  b  (t  -  u)  -  u)j,  (3-101) 

where  K  -  a  2  ~  2  constant  matrix  ~  is  tire  white  no»se  process  covariance  matrix,  and  the 
2x2  matrix  Dy  Q  satisfies  the  following  matrix  Fiedhoim  Integral  Equation: 

Dc(t,  a)  •  K  +  /2  D0(t.  s)  *  Kv.  (s  a)  *  ds  =  K.,  U,  u>.  (3-102) 

T.  '  ^ 


To  solve  fer  Dq  (...)  front  Equation  (3-102)  we  first  compute  the  chaff  dipoles 
reorired  signal  covariance  matrix  Ky  (l,  irt. 

We  have  -  from  Equation  (3-86)  of  the  previous  section  - 

(t  -  u.  r)  =  p  (l  -  u)  •  N (t)  '  P  (3-103) 

then 


OO 

«yc  ft.  u)  -  (t  -  «».  t)  *  X(t  -  r.  u  -  r)  *  dr 

•where  we  recall  that  ^  (*)  is  the  4  x  l  vector  with  entries  of  ^  (*).  and 
X{i  -  r,  u  -  r)  =  {x jit  -  r;  Xj(u  -  r)  ,  \j(t  -  r)  x2(u  -  r)  . 


(3-104) 


(3-105) 


x2(t  -  f)  Xj(u  -  r)  .  x2C*  ~  f)  x2(u  -  t)  }. 

Here  :»j  and  x2  are  entries  of  the  transmitted  signal.  We  have 

Ky^  (t,  u)  =  pft  -  u)  •  p  *  X(t  -  r.  t»  -  r)  •  N(-)  -  dr  (M06> 

=  p(t  -  u)  •  p  •  G*t.  u)  (say) 


Now.  let  gjj,  i.  J  =  1,  2,  be  the  entries  of  the  2  x  1  vcc**»r  G.  then 

f~ 3pjjU,  u)  +  g22(t.  u)  gj-^t,  u)  +  g-*j(t.  u)  1 
*yc  <*.  *0  =  p(t  -  «>  j  ~  "  \  p, 

Lgj-»(t,  u)  +  g?j<t.  u)  g  { *  ft,  u)  +  3g-*-»(t.  u)  J 

(3-107) 

=  0.9X2 


where  p 


3.3.2  The  System  of  Fredholm  Integral  Equations  for  Dq(i  u)  =  |Oq^I.  u>| 

Substituting  Equation  (3-108)  into  Equation  (3-104),  we  obtain  the  following 
system  of  integral  equations: 


“»  T 

pit  -  u)  Gj  j(t,  u)  p  =  D0n(t,  u)kj  j  +  Dq1-^,  u)k-»j  +  p  /  -  JDq^(i.  s)  Gjj(s.  u) 

T, 

+  D^-ft.  s)  Gp(s.  u)J  pis  -  u)  •  ds  <3-108i 

p(t  -  u)  G^d,  u)  p  =  Dq1  *(t,  u)k,  +  DQ,2(t.  u)k-,-,  +  u  |DqJ  *(t,  s)  Gp(s.  u) 

T, 

+  DQ12(t.  s)  G22(s-  Pts  -  u)  -  ds  J3-109) 

X 

p(t  -  u)  G21(t,  u)  <1  =  D02!(t,  a)kj  j  +  Dq"(».  uik-,,  *  -i  f  2  [D025(t.  s)  Gjjts.  u) 

*1 

+  Dq*2^,  s)  G2i(s,  u)J  pis  -  u)  •  ds  *3-1 10) 

^  X 

pit  -  u)  G2a(t,  u)  if  =  D0r‘(t,  u)kp  +  Bq— (t,  u)k-,2  +  p  /  2  |D02!(t.  s>  G,  ds.  u) 

Tj 

+  s)  G22(s.  u)l  pis  -  u)  -  ds  t3-l  n  > 


where  ky,  j  =  1,  2,  are  entries  of  the  constant  matrix  K.  and  Gy.  i.  j  -  !.  2.  ue 

entries  of  the  matrix  on  the  right  hand  side  of  Equation  (3-106). 

3.3.3  Integral  Equation  for  the  Filter  Eq 

The  optimum  filter  Eq  which  gives  the  minimum  variance  estimate  of  the  signal 
y-j-  -  under  Hj  -  is  given  by  the  matrix  integral  equation 

J 

EqO,  u)  *  K  +  / T2  Eq(i,  s)  JKy^  (s.  u)  +  Ky^  (s.  u)]  -  ds  -  (*.  u) 

1  (3-112) 


H-  re  Ky^.  (*)  is  the  target  received  signal  covariance  matrix  -  which  can  be  computed  by 

the  results  of  Section  3  2.  The  above  can  be  expanded  into  a  system  of  integral  equations, 
as  was  DqK,  u)  for  the  chaff. 
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•si 


J 


3.3.4  Approach  to  Solving  Equations 

The  design  of  the  Correlator-Detector-Receiver  is  complicated  by  the  fact  that  it 
involves  solutions  of  two  matrix  Fredholm  Integral  Equations.  It  is  clear  that  the  above 
integral  equations  can.  in  general,  be  quite  difficult  to  solve  “exactly”.  Their  solutions 
depend  on  the  functions  Gy  which,  in  turn,  are  dependent  on  the  transmitted  signal  x(t). 

In  most  cases,  these  solutions  can  be  approximated  by  appropriate  ortl.cnormal 
sequences  of  functions  (a  mean  square  approximation),  where  the  choice  of  these  functions 
depends  on  each  particular  case.  In  the  case  of  the  chaff  cloud  for  instance,  since  the 
interval  JTj ,  T-»j  is  taken  to  be  finite,  we  assume  that  the  functions  Dgy(t.  u)  can  be 
expanded  in  a  “double”  orthcnormal  series  of  Legendre  polynomials  P:(t)  Pj(u): 

Doij(t.  u)  =  £  £  nj  PK<t)  -  p,(u) 

0  0 


Substituting  in  the  system  of  integral  equations  we  then  oht.-;r>  a  system  of  algebraic 
equations  for  rhe  coefficients  jjk  j  y.  Sitnihrly.  one  can  approximate  the  integral  equation 
for  Egit.  u). 

It  is  also  noted  that  an  expedient  choice  of  the  transmit  signal  can  s.mpljfv  the 
problem  somewhat.  Referring  back  to  the  system  of  equations  for  Dq.  we  have 


gjdt)  =  /  N|t>  *  x;(t  -  7)  x^t  -  t)*  dr  i.  i  =  !.  2. 

y  — oo  ■  j 


(3*113) 


where  xKt).  i  =  1.2.  are  ihc  entries  of  the  transmitted  sisnai  X(t).  Since  Xm  can  be 

a  "*■ 

chosen,  and  due  to  the  specific  form  cf  the  function  X(r).  in  this  ease,  we  can  take 


V 


Xj(t)  =  Heg(t)  =  I.  and  x-»(t)  =  Hej(t)  =  {yfcioJX. 
the  first  two  Hcrmite  Polynomials.  Then  it  can  be  easily  verified  that 


<3-i  14) 


gjjit)  =  0  for  i  f  j 
P  0  for  i  =  j 


(3-115) 


.  i.  j  =  1.2. 


The  matrix  G(t)  therefore  becomes  a  diagonal  matrix,  which  will  simplify 
the  solution. 
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33.5  Fredholm  Approach  Summary 


The  above  approach  has  been  the  least  developed  during  the  course  of  this  study, 
ending  essentially  with  the  mathematical  formulation  of  the  problem.  The  difficulties 
envisioned  in  solving  the  sets  of  integral  equations,  would  tend  to  fav<  r  ihe  finite  dimensional 
and  state  space  approaches  as  the  more  preferred  methods  in  obtaining  solutions  to  waveform 
design  and  associated  coherent  receiver  processing.  Additions;  stud>  would  be  required  to 
verify  this  assumption. 

3.4  STOKES  FORMALISM 

The  final  method  reported  upon  for  utilizing  scattering  properties  of  target  aro 
chatter,  is  the  formulation  of  the  system  moo  el  in  terms  of  Stokes  parameters.  This  method 
is  summarized  briefly  below.  A  much  more  detailed  description  of  this  approach  can  be 
found  in  Rosien  et  al  1 1 9~9 1 . 


3.4.1  Slokes  Model 

The  previous  sections  have  described  the  general  concepts  of  scattering  theory  in 
terms  of  the  target  scattering  matrix,  as  given  by: 


H 


(3-116) 


Kennaugh  I1951J  has  shown  that  the  average  received  power  can  be  expressed  as: 


P  =  VAX 


(3-H7t 


where  X  and  Y  arc  four-dimensionai  vectors,  whose  elements  arc  the  Stokes  parameters, 
which  describe  the  transmitter  and  receiver  antenna  polarizations.  The  matrix  A  is  4  x  4 
real  symmetric,  and  is  known  as  the  average  Stokes  scattering  operator.  Tire  elements  of 
A  are  cor  posed  of  the  expectations  of  various  combinations  of  the  dements  of  the  scattering 
matrix  h  (see  Rosien  1 1 979 J  for  details;. 

If  one  is  given  the  Stokes  target  scattering  operator  A.  and  a  similar  operator  for 
clutter  C.  then  the  problem  becomes  one  of  maximizing 


YX\ 

=  s/c  «' 4-1181 

YCX 


The  solution  to  Equation  (3-118)  is  detailed  in  Rosien  |i979|.  and  consists  of  a 
constrained  maximization  of  the  ratio  of  two  bilinear  forms.  The  output  is  the  optimum 
pair  of  Stokes  vectors  for  X  and  Y.  which  completely  describes  the  transmit  and  receive 
antenna  polarizations. 

Results  of  this  approach  will  be  reported  upon  in  Section  5.2. 
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Section  4 


MODELING  AND  SIMULATION 

This  section  contains  a  brief  description  of  the  target  and  chaff  modeis  utilized  in 
the  computer  simulations  which  provided  the  performance  results  reported  in  Section  5. 

Both  the  target  and  chaff  models  have  been  developed  in  prior  work,  and  considerable 
additional  information  can  be  found  in  Rosien  { 1979J  _  The  major  items  of  interest  for  this 
study,  as  reported  below,  are  the  method  for  computing  the  target  scattering  covariance 
matrix  and  development  of  the  autoregressive  chaff  model  for  the  computer  simulation  runs. 
Also  included  are  brief  descriptions  of  the  various  software  packages  utilized  to  generate  the 
dual  and  single  channel  target  detection  performance  results. 

4.!  TARGET  MODELS 

4.1.1  Methodology- 

Since  exact  solutions  to  the  electromagnetic  response  of  complex  radar  targets,  such 
a<  sifCTaft.  are  almost  impossible  to  obtain,  most  studies  on  modeling  of  the  cross  section  for 
such  targets  employ  approximate  techniques,  which  are  valid  and  sufficient  for  the  frequency- 
range  of  interest  and  intended  applications.  We  are  primarily  interested  in  the  high  frequency- 
response  of  these  targets  (i.e..  wavelengths  smaller  than  the  characteristic  dimensions  of  most 
of  the  scatters  cn  the  targets).  Therefore,  we  have  assumed  that  the  target  response  can  be 
approximated  by  the  algebraic  sum  of  the  responses  of  a  collection  of  individual  scatters 
which  form  the  target.  Furthermore,  to  reduce  the  computation  labor  we  have  approximated 
most  of  the  scattering  components  of  aircraft  targets  by  simple  geometric  shapes  for  which 
there  are  either  exact  cross  section  expressions  or  approximate  expressions  derived  using 
Physical  Optics  (PO)  and, or  Geometric  Theory  of  Diffraction  (GTD)  techniques.  Both  of 
these  methods  are  high  frequency  approximations  and  arc  adequate  for  modeling  of  the 
scattering  characteristics  of  most  scatters  at  radar  frequencies. 

4.1.2  Target:  A  Collection  of  Scatters 

The  simple  geometric  shapes  which  we  use  in  modeling  aircraft  targets  arc  ellipsoids, 
elliptical  flat  plates,  cylinders  and  thin  wires.  In  addition  to  these  we  represent,  as  discrete 
scatters  surface  discontinuities  such  as  the  wing  fuse  age  joints  or  the  engine  fuselage  joints. 
Table  4-1  lists  the  various  types  of  scatterers  which  we  have  modeled  and  the  corresponding 
theory  (PO  or  GTDj  used  in  deriving  the  analytical  representation.  Figure  4-1  depicts  the 
BQM-34A  Drone  model  used  in  our  simulation  program.  Table  4-H  lists  the  individual 
scatterers  and  corresponding  scattering  models  used  in  the  BQM-34A  Drone. 

The  mathematical  model  of  each  scattcrer  is  embodied  in  ihe  j  -»tcring  matrix.  It 
is  represented  in  spherical  coordinates  in  the  scatterer  principle  axis  system,  which  is  selected 
such  that  a  diagonal  scattering  matrix  results.  In  mathematical  form,  wc  have 
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TABLE 4-1.  SCATTERER  TYPES  USED  FOR  THE  SIMULA TIOS 
OF  THE  RADAR  CROSS  SECTIONS 


Scatter «  Type 

Ellipsoid 

Elliptical  Flat  Plate 
Hollow  Elliptical  Cylinder 
Cylinder 
Thin  Wire 

Joint  Between  Two  Cylinders 
Exhaust  Duct 
Cone- Cylinder 

T  Doubly  Reflected  Ray  From  Two  Ellipsoids 

Cl 

g  Smooth  Junction  of  Curved  Suf races  (Slope  Discontinuity) 


Theory 

PO 

GTD 

GTD  &  PO 
GTD  &  PO 

From  Series  Expansion 
GTD 
GTD 
GTD 
PO 
GTD 


TABLE 4-U.  SCATTERER  TYPES  USED  IX  MODELING  THE  BQM-34A  DRONE 

Scatterer  Model 

Fuselage  Sections  1  and  3  Ellipsoid 

Fuselage  Section  2  Cylinder 


Wings.  Ta:'.  Stab  ’  -tcrs.  and  Vertical  Fin  Elliptical  Plate  ♦  Thm  Wire 

Wing;.  Tail.  Stabilizer,  and  Vertical  Fin  Joints  Joints 


Slope  Discontinuity  at  Engine  Interlace  and  Exhaust  Slope  Oiscon’inuity 

Engine  Open  Ended  Cylinder 


Engine  Blades 


Two  Thin  Wires 


Double  Reflections  Two  Ellipsoids 

Wing  -  Fuselage 
Stabilizers  -  Tail 
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0  0  0 


si  =  0  %i  0  * 


H^kpi 


where  Sj  is  the  scattering  matrix  of  the  i^!  scatterer  in  its  own  principle  axis  system.  and 
Spj  are  complex  functions  of  the  scatterer  size  3nd  aspect  angle  in  the  orthogonal  6  (vertical) 
and  p  (horizontal)  directions,  and  pj  is  the  one-wa>  distance  beiween  a  plane  wave  front 
passing  through  the  reflection  point  and  one  passing  through  the  center  of  the  target 
coordinate  svstem. 


To  compute  the  received  signal,  we  transform  the  scattering  matrix  from  the 
scatter  principle  axes  system  to  radar  coordinates.  This  includes  several  transforms  as  follows: 

Hj  =  P'A'B'Q'Sj  QBAP  (4-2) 

where  Hj  is  the  scattering  matrix  of  the  i^1  scatterer  in  the  antenna  coordinate  system.  P  is 
a  transformation  from  spherical  to  Cartesian  coordinates  in  me  antenna  frame.  A  is  a 
transformation  from  the  antenna  frame  to  the  target  frame.  B  is  a  transformation  from  the 
target  frame  to  the  scatterer  principle  axis  frame,  and  Q  is  a  transformation  from  Cartesian 
to  spherical  coordinates  in  the  scatterer  frame. 

In  the  case  of  a  single  pulse  much  longer  than  the  target,  the  scattering  matrix  of 
the  target  may  be  approximated  by  the  summation 

HT  =^CHi  ‘4-3> 

The  three  element  poian/awon  signal  vector  backscattcred  from  the  tareet  is  then  given  bv 

Hb  =  HTH,  (4-1) 

where  the  subscript  b  refers  to  the  backscattcred  target  return  in  the  antenna  frame,  and  the 
subscript  t  refers  to  the  transmit  antenna  polarization  sector. 

■*•1-3  Radar  Cross  Section  Validation 

The  target  model  of  the  BQM-34A  drone  was  validated  by  comparing  the  simulated 
RCS  plots  ror  variou*  transmit  and  receive  polarizations  with  RATSCAT  [  J977|  measurements. 
The  drone  model  consists  of  approximately  30  individual  scattcrers.  each  of  which  includes 
both  size  parameters  and  a  fcackscatter  coefficient.  The  size  parameters  were  selected  to 
resemble  the  target  configuration.  Then  from  a  knowledge  of  the  characteristics  of  each 
scat.erer  model  and  through  a  trial  and  error  process,  the  backscatter  coefficients  were 
adjusted  in  order  to  match  the  simulated  RCS  plots  to  those  obtained  from  measurements. 
Sample  results  of  the  comparisons  arc  shown  in  Figures  4-2.  4-3.  and  4-4  for  three  different 
transmit  and  receive  polarization  settings.  As  can  be  seen,  the  simulation  results  match  the 
measurements  rather  well  in  the  a,cas  of  general  curve  drape,  specular  type  returns,  lobins 
structure,  fluctuation,  and  relative  scale  (power  levels). 
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4. !  .4  Covariance  Matrix 

The  covariance  matrix  of  the  target  was  computed  using  300  samples  of  RCS 
obtained  from  contiguous  increments  in  aspect  angle  The  angle  stepping  was  performed  in 
nested  looos  with  yaw  being  stepped  in  the  inner  loop:  roll  in  the  intermediate  loop:  and 
pitch  in  the  outer  loop.  Yaw  was  varied  between  5  and  25  degrees  with  increments  of 
0.2  degree.  Roll  ranged  from  -30  to  -28  degrees  and  was  incremented  in  steps  of  1  degree. 
Pitch  was  varied  from  9  to  1 1  degrees  in  1-degree  steps.  These  3$pect  variations  actually 
yield  9C0  samples,  of  which  300  were  selected  for  covariance  matrix  computation.  When  the 
sample  time  was  assumed  to  be  1  msec.  300  contiguous  samples  were  used.  When  the  sample 
time  was  assumed  to  be  3  msec,  every  third  sample  was  used  to  compute  the  covariance 
matrix.  (Note  should  be  made  that  there  are  four  sets  of  RCS  samples  of  900  numbers. 

Each  set  represents  one  of  the  elements  of  the  scattering  matrix  <hjj.  bp.  h^ j -  ho-O). 

Plots  of  the  correlation  function  derived  from  the  simulated  target  modei  are  shown 
in  Figure  4-5  for  both  the  amplitude  and  phase.  Since  h->j  =  hp,  h-»j  is  not  explicitly 
shewn. 

4.2  CHAFF 

4.2. 1  The  Dipole  Chaff  Cloud  Model 

In  this  section  a  chaff  cloud  is  modeled  as  a  collection  of  rotating  dipole  scatterers. 
The  modei  provides  for  simulation  of  the  polarization  and  spectral  characteristics  of  radar 
signals  backscattered  from  a  cloud  of  rotating  dipoles  which  can  have  either  completely 
random  or  any  preferred  orientation.  A  vector  autoregressive  process  based  upon  the  above 
mode!,  is  then  used  to  obtain  a  Stochastic  model  for  the  simulation  of  the  scattering  matrix 
elements  of  the  dipole  cloud  and  signals  scattered  from  the  chaff. 

4.2.2  Theoretical  Covariance  Matrix 

A  description  of  the  theoretical  aspects  of  the  chaff  covariance  matrix  is  presented 
below  (more  detail  is  given  in  Appendix  A).  This  covariance  function  can  be  written  in 
terms  of  the  spectral  and  polarization  characteristics  of  the  chaff  cloud.  These  characteristics 
are  dependent  upon  the  geometrical  configurations  of  the  dipoles  which  constitute  the  cloud. 
As  a  consequence,  the  orientation  of  a  typical  dipole  (as  depicted  in  Figure  4-6)  will  be 
considered  first  in  the  development  of  the  theoretical  model. 

To  arrive  at  a  description  of  the  dipole  orientation,  we  carry  out  a  transformation 
consisting  of  three  successive  rotations  defined  by  the  Eulcrian  angles.  Starting  with  the 
unprimed  system,  we  rotate  the  double  primed  system  by  an  angle  q  counterclockwise  about 
the  z  axis.  The  primed  system  is  defined  through  a  rotation  by  an  angle  0  in  the 
counterclockwise  direction  around  the  x’  axis.  The  dipole  orientation  is  then  found  in  the 
primed  system,  where  the  axis  of  rotation  coincides  with  the  x'  axis  and  the  dipole  is  in  the 
(x\  y')  plane  as  defined  by  the  angle 
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In  order  to  be  consistent  with  our  definition  of  the  scattering  matrix,  we  orient 
the  -:np:.mcd  system  such  that  x  is  in  the  vertical  direction  and  v  is  m  the  horizontal 
direction.  The  wav.*  is  assumed  to  travel  along  the  z  axis.  The  reoriented  coordinate  system 
is  shown  in  Figure  4- 7  where  z  is  into  !he  page.  Notice  also  that  the  primed  sysicrn  is 
shown  with  the  dipole  rotating  in  the  (s'.  >')  plane.  The  primed  stnera.  as  shown  in 
Figure  4-7.  is  derived  from  Figure  4-o  bv  setting  p  =  2?(f  ~  Ao  =  -9*3'  *  ip.  and 
0  =  90s  +  it?.  The  angle  u  la  function  of  time)  is  g;vcn  by  y  =  ec.t  ■«  o.  where  cc.  is  the 
dipole  rotation  rate  and  a  S  the  initial  angle.  As  seen  .n  the  figure-,  she  perturbations  AO 
serve  to  orien*  the  axis  of  rotation  is')  and.  subsequent!),  the  plane  of  the  mating  dipole 
relative  to  the  x  axis.  Notice  that  with  AO  =  0  the  line  of  sight  iz  axis*  remains  in  the 
plane  of  rotation.  Only  when  A*?  0  i*  the  axis  not  in  the  plane  of  rotation.  The  limits 

on  the  perturbations  are  as  follows.  -90*  ^  it?  <  SO"'  and  —U>0'  -7  ip  <  ISO  .  if  t!> 
limits  are  reached,  the  dipole  may  have  any  orientation  whatsoever. 

A 

•f  the  dipole  direction  (denoted  by  the  uni:  vector  d»  ri  described  in  the  rx  .  y'j 
phne  bv  the  angle  p.  then  we  can  find  its  orientation  in  the  (x.  y.  z»  coordinate  system 
through  the  application  of  the  inverse  Fulcrum  transformation.  Tile  scattering  matrix  of  such 
a  dipoic  to  an  electromagnetic  wave  props  faring  clone  the  z  avis  N  given  oy  (Borison  (!96”i. 
Wong  et  ai  |  i  j  • 

f  hvv  fev  .  , _ 

H  =  1  '  =\0  9  >.dd*  14-5) 

lb*  hyy. 

A  t* 

where  >.  is  the  wavelength  of  the  incident  signal,  du'ts  defined  as  the  outer  product,  and 
h.vv  ~  i'svx-  *’roni  *s  clear  that  hyx.  hxv.  hvv  and  hvx  are  ill  functions  of  the 

variables  d.  q.  and  d-  and  X.  where  v  -  c;rt  +  a  is  a  function  of  the  dipole  rotation  mle  and 
time  The  orientation  of  the  dipole  rotation  plane  described  by  0  and  o  impact  the 
polarization  of  the  scattered  wave,  while  the  Doppler  is  dependent  upon  c  and  X 

before  proceeding  to  determine  the  temporal  and  pchrizatioxi  characteristics  of  a 
collection  of  dipole*  kt  us  define  a  new  vector 

h««;  ~  |hxxt»).  hxy(ts.  hyyfttj  (4-r>i 

Thcri  the  covariance  matrix  of  h(t>  is  given  by 

Kj|(r)  -  fcth(t)  h*tt  -  r)J  t4-“) 

•vhere  the  *  denotes  complex  conjugate  transpose. 


When  the  angular  rotation  rate  and  the  Doppler  due  to  the  drift  velocity  are  each 
described  by  distribution  functions  which  include  means  and  standard  deviations,  then  Kjjtr) 
wiil  be  a  function  of  to.  o..  tu^.  anti  o^.  where  the  subscript  r  denotes  rotation,  the 
subscript  d  refers  to  Doppler,  the  bar  means  average,  and  a  is  the  standard  deviation.  If  the 
dipoles  are  distributed  over  reel  rioted  .egion.*  of  both  $  and  <).  Kj«(r)  will  thus  be  a  function 


of  these  limits  6 j.  6 dj.  and  p-,.  As  a  result  the  covariance  matrix  includes  spectral 
content  due  to  the  dipole  activity  and  polarization  characteristics  due  to  the  dipole  distribution. 

The  covariance  matrix  can  be  written  such  that  the  spectral  and  polarization 
characteristics  are  separable  functions  tsee  Appendix  A): 

Kh(t)  =  Vd(T>  •  j  C  +  MjX2r)  -  D  +  M^’t)  •  B  J  (*-S) 

v. fiere  Md(r)  =  Nld( r.  u>^.  o^)  and  M,ir)  =  Mjir.  uT.  oT )  are  scalars,  and  where  C  and  D  are 
3X3  matrix  functions  of  0j.  0->.  Pj.  and  P-».  Notice  that  the  spectral  content  for  all  of  the 
elements  in  the  covariance  matrix  stems  from  the  same  two  sources,  namely  \!d  and  Mr 
However,  the  elements  may  exhibit  different  correlation  characteristics  as  a  function  of  the 
matrices  C  and  D.  In  effect,  the  correlation  functions  for  all  of  the  elements  in  Kh<t> 
result  from  the  summation  of  weighted  time  furctions.  where  the  weights  differ  for  the 
various  elements. 


4-2.3  Stochastic  Model  Derived  From  Autoregressive  Process 

Given  the  theoretically  derived  m'trix  covariance  function  Kjj(t)  of  the  chaff 
scattering  matrix  elements  ky.  the  stochastic  vector  hft)  mav  be  modeled  by  a  three  element 
complex  vector  autoregressive  process  excited  with  a  zero  mean  complex  vector  white 
Gaussian  noise  source  u(t)  described  as  follows. 

M 

hft)  =  2  Aj  hit  -  »At)  +  aft)  (4-9) 

i=J 

where  the  Aj's  are  3  x  3  complex  matrix  coefficients  to  be  determined.  The  matrix 
coefficients  A;  can  be  obtained  by  solving  the  matrix  Yule-Walker  equations 

M 

K}|fmAt)  =2  -Aj  Kh  f(m  -  i)At}  m  =  1 . M  (4-10) 

i=l 

As  the  order  of  the  autoregressive  (AR>  model  increases  by  I  from  M  to  M+l.  there  arc 
nine  new  matrix  elements  of  +  j  to  be  computed  in  addition  to  the  9  x  M  elements  of 
Aj  through  A^j  that  need  to  be  recomputed. 

!t  is  noted,  as  stated  previously,  that  the  Kjjfr)  in  Equation  (4-S)  can  be  decom¬ 
posed  into  three  components  as  foiiews. 

Kjjfr)  =  (Md(r)  •  C)  +  (Mdfr)  *  M/2T)  *  D)  +  (Md(r>  -  Mf<2r)  *  D) 

=  kjjjfr)  +  ^H2-t^  +  (4-11) 


i 


It  can  be  shown  that  the  covariance  function  of  the  sun;  of  the  outputs  of  three 
independent  linear  systems  each  being  excited  with  a  mutually  independent  zero  mean  while 
Gaussian  noise  source  is  equal  to  the  sum  of  the  autocovariance  functions  of  the  output  of 
the  three  linear  systems. 

Hence,  each  of  the  covariance  functions  KHrfT)  (i=i  .2.3)  cat  be  individually 
modeled  with  an  autoregressive  process. 


AHo.  we  note  that  Kj|j(r>  is  composed  of  ?  product  of  the  time  correlation 
function  M^u).  which  i=  a  scalar,  and  the  covariance  matrix  C  of  the  scattering  matrix 
elements  hy.  it  then  follows  that  the  complex  vector  autoregressive  process  can  foe  modeled 
with  scalar  coefficients  a- _ .a^j.  that  is 

A*  =  aji  ;=i.  ..  Ma 


!}(t*-^!l  Aj  hj(t  -  :At)  +  U] 


(4-1 2i 


a  a 

KHj(mAt)  =  2  ai  KH1{(m  -  iiAt!=  C  2  ^  •  Md((m  -  iiAtl 


The  same  argument  applies  to  K^-dr)  and 


Therefore. 


h?(D  =  2/  ^  1,Ml  ~iAt>  +  U7,,) 

i  =  1 


(4-13) 


Kn^roAt)  -  2  b-  Kjj-»I(m-  OAt] 


-  D  -  ^  foj  •  -Mdf(m  -  i)A|j  -  Mrl2(m  -i)At  j 


and 


He 

l»3<t)  =  2  Cj  h5(i  -  LAt>  +  113(1) 

i=l 


(4-14) 


Kj|3(mAt)=2  -  i)At} 

i=l 


M 

_  C 

^  D  *  2  Cj  Mdl(m  -  i)  At)  -  -  i)At) 

i=l 

Furthermore.  since 

D  -  M^)  +"d  -  M/2 t)  =  2  Rc)D  - 
it  follows  that  Cj  is  the  complex  conjugate  pair  of  bj  Le.„  Cj  -  bj 
and 


Mr)  =  hifT)  +  h2(r)  +  h3(r)  C4-15) 

The  scalar  coefficients  aj,  bj  and  Cj  can  be  obtained  by  solving  the  Yule-Walker 
equations  of  the  condition  functions  M^fr).  M^fr)  -  Mjfcr)  and  M^(r)  *  Mj(2t)  respectively. 

From  ja  closer  examination  into  the  properties  of  the  correlation  functions 
Md(r).  Mj(2r),  Mji  2r)  *n  Appendix  A  (Equations  10-14).  it  is  noted  that  these  functions  are 
expressed  in  terms  oc  exponential  functions  of  complex  exponents. 

Md(r)  =  exp  Hccjf  -  -7—  t2] 

Mj(2r)  =  exp  (-j2 ui^r  -  2or-r^J 
Mj(2t)  =  exp  IjZwjT  -  2or-r-J 

where  oj  denotes  the  average  value  of  w.  and  o~  represents  the  variance  of  «. 
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The  real  part  of  the  exponenis  are  ail  negative  and  thus  an  exponentially 
deco  Relating  function  of  time.  The  imaginary  part  un  is  associated  with  the  sinusoidal 
modulation  of  angular  rate  to  (i.e..  m(r)  =  exp(-jo;t).  a  line  in  the  power  spectral  domain 
since  this  imaginary  part  of  the  exponent  is  a  linear  function  of  t. 


Under  this  special  case,  the  real  and  imaginary  parts  of  the  exponent  can  be 
separated  by  associating  the  real  part  with  the  autoregressive  process  and  the  imaginary 
part  with  the  sinusoidal  modulation. 


Realization  of  the  autoregressive  process  for  the  general  case  of  Equations  (4-12). 
(4-13)  and  (4-14)  yields  complex  valued  coefficients.  For  the  special  case  the  coefficients 
ar*  real. 

A 

A  For  this  special  case  we  define  the  three  au*ore«ressive  <AR»  processes  as  «i.(t). 

h-,(t)  and  fi^(t)  and  the  associated  scalar  coefficients  aj.  Oj  and  Cj  which  are  real  values. 


Thcn  we  have 


Ma 

A  „  AA 

h j(t>  =  2J  3jhj<t  ~  +  Ujft) 

i=l 


(4-16) 


M„ 


M. 


A  „  AA  __  a  °d~  t 

K|||(mAt)  -  2^  a-K|^j{(m  -itAt)  =  C  ^  ai  *  CXP  l~ — -  i)~^t~) 


i=l 


i=  1 


A 

hMt) 


„  AA 

2^  bjh-«(t  -  iAt)  +  u->(t) 
i=l 


(4-17) 


A 

k|p(mAt)  — 


t=l 


A  A 

Oj  KjpKm  —  i)At] 


Mb  a 
»=  1 


Oj"  +  4<J  “ 

_ 1 _ L  <m  -  i)-At_l 


I) 


(4-18) 


A  A  A 

h3(t)  =  X  bi  h3(t  *  ,A,)  ^  “3<t) 
i=l 


Mb 

A  •r-^  A  A  .  . 

KH3<mAt)  =  X  bi  kH3l(m  '  ,)Al* 

i=l 


Mk 


m  (J  j  J  4o.* 

=  D  V  b-  exp  l-  — - (rn  - 

*■  i  •> 


-At-? 


;=? 


and 


~j“;r*yy 

h<t)  =  e  lhj(t)  +  c  hR(t)  +  e  h3(t>l  <4"19 

An  efficient  alsoriihm  developed  by  Durbin  i  DURBIN,  19601  :o  compute  the 
coefficients  of  the  AR  process,  is  used.  This  algorithm  computes  the  coefficients  in  a 
recursive  manner  increasing  the  order  of  the  AR  process  by  one  with  each  iteration  until 
the  desired  mean  squared  error  criterion  of  the  fitness  of  the  AR  model  to  the  theoretical 
compliance  function  is  met.  A  discussion  on  the  Durbin  Recursive  Algorithm  is  grren  -y 
Appendix  F. 

Realization  of  the  stochastic  modeling  of  chaff  is  accomplished  using  the  AR 
processor  in  the  form  of  a  digital  filter  with  an  appropriate  number  of  delays  T  and  an 
associated  weighted  feedback  k>rp  whose  weights  are  the  scalar  coefficients.  Two  functional 
diagrams  of  the  stochastic  chaff  model  are  given  in  Figures  4-8  and  4-9  for  the  general  and 
special  cases  respectively. 


4.2.4 


Covariance  Matrix  Validation 


Correlation  functions  of  the  elements  in  the  chaff  scattering  matrix  are  plotted  in 
Figures  4-10  and  4-11.  These  functions  were  derived  from  the  covariance  matrix,  which 
resuits  from  the  scattering  matrix.  Figure  4-10  lias  plots  of  the  theoretical  correlation  functions 
according  to  the  results  of  Appendix  A.  Simulated  correlation  curves  for  a  7-point 
autoregressive  fit  of  both  amplitude  and  phase  .*re  found  in  Figure  4-1 1.  As  would  be 
expected,  the  simulated  curses  differ  slightly  from  the  theoretical  curves.  However,  the 
simulated  curves  would  be  expected  to  approach  the  theoretical  curves  if  an  infinite  number 
of  samples  were  used  in  the  computation. 

4.3  SIMULATION 

4.3.1  Overview 

The  overall  simulation  is  carried  out  in  two  major  steps  as  indicated  in  Figure  4-12. 
The  first  step  consists  of  simulating  the  target  and  chaff.  Included  in  this  function  is  the 
generation  of  sampled  versions  of  scattering  matrices  for  both  target  and  chaff  as  would  be 
sensed  by  the  radar.  The  second  major  step  is  the  simulation  of  the  signal  processing 
performed  in  the  radar.  This  consists  of  computing  the  covariance  matrices  of  the  target  and 
chaff  from  the  sampled  dat2  and  performing  the  optimization  process  to  arrive  at  the  weights 
(and  waveform  in  the  matrix  weighting  case)  to  be  applied  to  the  pulse  train. 

This  simulation  approach  has  the  advantage  that  the  target  and  chaff  simulations 
need  to  be  run  only  once  for  a  given  set  of  conditions.  The  results  are  stored  on  tape  and 
can  then  be  used  as  needed  by  the  various  processing  routines. 

4.3.2  Target 

Simulation  of  the  target  consists  of  generating  samples  of  the  scattering  matrix  for  a 
series  of  aspect  angles  representative  of  a  target  maneuver.  Two  basic  modules  are  used  in 
the  simulation,  each  of  which  includes  several  subroutines.  One  is  the  geometry  module:  the 
other  is  the  scattering  matrix  module.  These  modules  and  routines  are  shown  in  the  target 
simulation  block  diagram  of  Figure  4-13. 

Two  loops  are  shown  in  the  diagram,  one  is  a  loop  on  all  of  the  scattcrer's 
composing  the  target.  The  second  ioop  controls  the  number  of  scattering  matrix  samples 
generated  by  the  program. 

The  input  to  the  program  includes  the  range  of  aspect  angles  and  angle  increments 
to  be  used  in  computing  the  scattering  matrix  samples.  The  output  consists  of  a  vector  sum 
of  the  individual  scattcrer  scattering  matrices.  This  sum  includes  both  amplitude  and  phase 
from  each  of  the  scatterers  for  eaco  of  the  elements  in  the  scattering  matrix.  This  is  what 
produces  the  deep  nulls  in  the  res  plots  as  a  function  cf  aspect  angle. 
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Figure  4-12.  Overall  Simulation  Block  Diagram 
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Figure  4-13.  Target  Simulation  Block  Diagram 


4.3.3  Chaff 


Chaff  simulation  consists  of  generating  random  samples  of  the  scattering  matrix 
from  representative  stochastic  processes.  A  block  diagram  of  the  chaff  simulation  is  shown 
in  Figure  4-14.  The  program  is  divided  into  two  basic  areas.  The  first  consists  of  a  one-time 
computation  of  parameters  and  values  needed  for  generation  of  the  random  samples.  The 
second  is  the  loop  composed  of  the  actual  steps  followed  in  the  recursive  generation  of  the 
stochastic  process  outputs. 

Inputs  to  the  program  include  parameters  to  the  theoretical  correlation  model,  anglt-s 
describing  the  geometrical  distribution  of  the  dipoles,  parameters  affecting  the  power  received 
by  the  radar,  and  the  number  (and  size)  of  tile  time  increments  required.  The  outputs  are 
samples  of  the  elements  in  the  scattering  matrix  as  a  function  of  the  time  increments. 

4.3.4  Processing  Algorithms 

4.3.4. 1  Weighting  Matrix  -  A  block  diagram  showing  the  main  steps  in  the  simulation  of 
the  weighting  matrir.  followed  by  a  maximum  likelihood  receiver  is  depicted  in 

Fisi  re  4-15.  Tire  inputs  to  this  processor  simulation  are  the  target  and  chaff  samples 
of  :he  scattering  matrices. 

The  first  section  of  the  program  is  devoted  to  computing  the  optimum  transmit 
waveform  and  receive  weighting/transform  matrix.  The  weighting 'transform  matrix  was 
computed  according  to  the  method  discussed  in  Section  3.  However,  the  transmit  waveform 
is  determined  through  a  more  direct  approach.  We  assumed  that  .an  antenna  in  a  practical 
system  would  be  able  to  transmit  polarization  states  defined  by  four  bits  of  phase  and  four 
bits  of  amplitude  information.  This  yields  a  total  of  256  states,  only  114  of  which  are 
nonredundant.  Thus,  the  algorithm  used  in  selecting  the  optimum  waveform  cycled  through 
these  1 1 4  states  for  each  puke  and  chose  the  state  which  minimized  the  Frobenius  norm  as 
discussed  in  Section  3. 

The  second  part  of  the  program  consists  of  a  loop  wherein  the  target  and  chaff 
samples  are  operated  on  by  the  optimum  radar.  This  includes  transmitting  the  optimum 
wavefom.  operating  on  the  return  with  the  weighting/transform  matrix,  and  comparing  the 
result  with  a  threshold.  The  number  of  threshold  crossings  is  then  counted  and  used  in 
determining  the  probability  of  detection. 

4.3.4.2  Weighting  Vector  -  A  block  diagram  of  the  vector  weighting  method  is  shown  in 
Figure  4-16.  The  simulation  consists  of  two  parts.  The  first  includes  a  loop  to  iterate  on 
the  optimum  weighting  vector  and  then  computes  the  signai-to-chaff  ratio  resulting  from  the 
application  of  this  optimum  vector.  The  second  part  computes  the  signal-to-chaff  ratio  for 
the  simple  matched  filter  case.  Plots  of  the  outputs  provide  comparison  of  the  optimum 
vector  approach  to  the  conventional  matched  filter  method. 
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Figure  4-N.  Chaff  Simulation  Mock  Diagram 


Diagram  of  Weighting  Matrix  Simulation 


I'llture  4- 16.  Mock  Diaxrmo  of  Vector  We  lull  non  Simulation 


4.3.4.3  Single  Chancel  -  A  simple  method  was  u  ed  in  performing  the  simulation  of  the 
single  channel  system.  The  approach  was  to  compute  the  target  and  chaff  covariance 
matrices,  introduce  a  waveform  with  equal  pulses,  account  for  matched  filter  weighting,  and 
then  to  calculate  the  theoretical  probability  of  detection  based  on  the  power  computations 
for  chaff  plus  noise  vs  target  and  chaff  plus  noise.  A  block  diagram  of  the  simulation  flow 
is  shown  in  figure  4-17. 
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Section  5 


COMPUTER  SIMULATION  RESULTS 

The  purpose  of  the  computer  simulations  was  to  evaluate  the  effectiveness  of  dual 
channel  systems  as  compared  to  single  channel  systems.  This  was  accomplished  by  exploiting 
the  polarization  properties  of  various  target  metiers  relative  to  those  of  chaff,  in  terms  of 
maximizing  the  receiver  output  sgnaTto-clutter  ratio  (s;'c).  The  major  fisure-of-mcrit  JFOMi 
considered  was  the  probability  of  targe?  detection  for  a  given  fake  alarm  probability  due  to 
dutter  plus  receiver  noise.  Associated  with  each  result  is  the  other  major  study  output  which 
yields  transmit  waveform  selection  coupled  with  optimum  receiver  weighting. 

The  results  reported  indude  those  of  the  matrix  and  vector  approaches,  as  well  as 
the  Stokes  formalism.  The  latter  was  noted  as  being  obtained  in  a  previous  study 
[Rooen,  1979] . 

5.1  SIMULATION  PARAMETERS 

This  section  describes  the  various  input/output  parameters  associated  w-th  obtaining 
results  for  the  aforementioned  three  approaches,  and  discusses  the  types  of  conclusions  one 
can  draw  from  each. 

5.1.1  Target  and  drrff  Models 

The  description  of  the  generation  and  simulation  of  The  target  and  chaff  models,  has 
been  covered  in  Section  4.  The  scattering  matrix  correlation  functions  of  the  models  utilized 
to  obtain  the  reported  upon  dud  and  single  channel  results  are  compared  for  reference  in 
Figure  5-1.  The  1  hsv  h^-**  1  parameter  is  shown. 

Three  target  models  are  noted:  the  standard  Swelling  i  and  two  others  designated 
as  RT1  and  R13.  The  time  base  for  the  correlation  functions  is  given  in  terms  of  an 
interpube  period  AT,  to  conform  with  the  use  of  transmit  signals  consisting  of  generalized 
puise  train  formats.  The  SWi  target  model  utilizes  the  same  scattering  matrix  as  the  BQM-34A 
model,  but  has  a  constant  correlation  function  (pube-to-pube  correlation!.  The  RT1  and  RT3 
target  models  are  also  those  of  the  BQM-34A  drone.  Tnis  shorthand  notation  refers  to 
interpube  sampling  intervab  of  1  msec  and  3  msec  respectively,  where  3  msec  represents  a 
greater  change  in  target  aspect  per  AT.  As  a  consequence,  the  RT1  model  b  more  closciy 
associated  with  the  Swcriing  I,  whereas  RT3  wouid  be  more  nearly  like  a  Swerling  II. 

The  chaff  model  dosely  matches  the  RT3  target  in  terms  of  its  hn  scattering  matrix 
correlation  function.  It  should  be  noted  that  the  same  chaff  correlation  property  is  used 
regardless  of  the  value  of  the  interpube  period.  In  this  way  one  obtains  a  variation  in  the 
spectrum  of  the  target  relative  to  the  duff.  For  example,  the  RT3  target  has  essentially  the 
same  spectrum  as  does  the  duff.  By  contrast,  the  chaff  has  a  much  broader  spectrum  than 
does  the  impulse  assodated  with  Swelling  I.  The  RT1  target  represents  a  choice  between 
these  two. 
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5.1.2  Definition  of  s/c  Rati 


Each  of  the  approaches  has  a  common  reference  for  s/c  ratio,  a  key  parameter 
associated  with  probability  of  target  detection.  The  ratio  utilized  is  defined  as  -he  total 
available  target  power  relative  to  the  total  available  chaff  power.  The  available  power  is 
aver,  as: 


?a  =  ih22?-  +  2!hI2«2-Hh:2!2  f5-l) 

where  the  H  parameters  are  relative  to  either  the  target  or  chaff  scattering  properties. 

Fquain  :i  (5-1 }  represents  all  of  the  energy  that  can  be  scattered  back  from  either  tarset 
or  chaff,  when  impinged  upon  hy  a  dual  channel  transmit  waveform.  Receiver  noise  is  taken 
into  account  by  scaling  Pa  ;o  a  normalized  noise  vector,  if.  for  example,  the  s-'n  ratio  is 
designated  as  20  dB.  then  Pa  for  the  target  would  be  100  times  greater  than  the  input  noise 
vector.  Similarly,  a  c.n  ratio  of  17  dB  results  ir.  a  Pa  for  chaff  which  is  50  times  the  noise. 
This  results  in  an  available  input  sv  ratio  of  -s3  dB.  which  is  then  the  value  plotted  on  the 
detection  probability  curves. 

5. 1 .3  Qualitative  Assessment 

There  are  various  d.fterences  noted  in  the  computer  simulations  for  the  three 
approaches  reported  upon.  The  Stckcs  formalism  compared  dual  and  single  hanne!  results  on 
the  bar,  of  single  pulse  probability  of  dele  tion.  As  a  consequence,  the  correlation  properties 
of  the  targe!  and  chaff  were  not  a  factor,  although  the  scattering  matrices  were  kev  parameters 
m  the  optimization  procedure.  The  Stokes  output  consisted  of  the  optimum  transmit  and 
receiver  antenna  poiariraiio&s  for  maximizing  the  s  c  ratio. 

The  vector  method  employed  a  transmit  sequence  of  five  puists  using  the  RT!  t arret 
model.  An  iterative  procedure  was  used  to  converge  to  the  maximum  s  c  ratio.  Comparisons 
were  nude  for  dual  and  single  channel  performance  based  upon  rebfiw  %  c.  but  not  to  targe, 
detection.  Other  significant  results  consisted  of  the  selected  five-puke  transmit  vector  and 
associated  dual  channel  receiver  weighting  function  which  maximized  the  s  c.  .As  previously 
noted,  these  reside  depend  upon  the  initial  choice  of  the  five-puke  waveform. 

The  matrix  resuits  were  bv  far  the  most  comprehensive,  comoanng  performance  for 
all  three  target  models  against  both  dual  and  single  channel  implementlionl  Results  weie 
derived  for  the  probability  ct  target  detection  Iigure-cf-merit.  Orhr.  outputs  of  interest  were 
the  optimum  .ransmit  vector  selected  from  among  a  library  of  1  waveforms  utilizh.g  a 
th.ee-puke  sequence,  and  it?  associated  6  \  6  receiver  weighting  matrix. 

Although  the  three  approaches  had  differences  associated  mo  choices  of 
transmit  vector  pulse  sequences.  the  results  were  cmsisicr.t  enough  to  be  agnifiowt  All 
three  exhibited  the  same  byric  result.  in  that  the  dual  channel  impkmn  ntat lou  does  provide 
enhanced  performance  over  conventional  single  oranne!  configurations. 


STOKES  APPROACH  RESULTS 


5.2 

The  above  was  reported  by  Rosien.  ei  ai  (1979).  Figure  5-2  is  excerpted  from  this 
report  and  is  representative  of  the  basic  results  obtained  using  the  Stakes  formalism. 

The  probability  of  target  detection  for  single-pulse  operation  is  shown  for  various 
sic  ratios  (as  defined  previously).  The  false  alarm  probability  was  set  to  iO~T.  the  chaff 
dipoles  were  uniformly  distributed  ±45  degrees  from  the  hori/.ental.  and  the  BQM-34A 
aspect  angles  are  distributed  as  shown  in  the  figure.  Optimum  dual  channel  performance  is 
plotted  relative  to  conventional  single  ctannel  operation.  Substantial  improvement  is  noted 
when  compared  against  single  channel  .lorizonta!  and  circular  polarized  s>  stems.  Approxi¬ 
mately  2  dB  of  improvement  is  seen  ve*»us  the  vertically  polarized  single  channel  system.  The 
lesser  improvement,  in  this  case,  is  -  result  of  the  chaff  model  dipole  distribution  which  favors 
the  horizontal.  The  vertically  polarized  signal,  therefore,  results  ir.  substantially  less  chaff 
bsckscattcr  than  does  either  the  horizontally  or  circularly  polarized  waves.  Since  the  target 
does  no!  have  a  strong  preferred  polarization  for  its  backscatter  retun;,  it  becomes  obvious 
that  a  vertically  polarized  system  would  be  the  choice  for  the  single  cnanne!  implementation, 
and  a  nearly  vertically  po’arized  system  would  also  be  implemented  for  dual  channel  operation. 


Figure  5-2  Sisj'le-Ful^r  Detection  Probabilities  for  various  ( trammit-fecrirr )  polarizations 
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The  major  fact  seen  by  these  results,  is  that  the  Stokes  method  does  yield  better 
performance  for  dual  channel  operation,  despite  the  obvious  high  level  of  performance 
exhibited  by  single  channel  vertical.  It  should  also  bi  noted  that  these  results  were  for 
single  pulse  transmission,  such  that  the  spectral  differences  between  the  target  and  chaff 
could  not  be  exploited  in  the  signal  processor. 

5.3  VECTOR  APPROACH  RESULTS 

The  vector  weighting  results  were  based  upon  transmission  of  a  five-pulse  sequence. 
The  figure-of-merit  utilized  was  the  relative  improvement  in  s-'c  for  various  dual  and  single 
channel  implementations,  as  a  function  of  the  chaff  dipole  distribution.  No  further  effort 
was  expended  in  this  area,  in  view  of  the  importance  attached  tc  die  development  cf  the 
more  encompassing  matrix  approach. 

Figure  5-3  depicts  typical  results  for  the  RT1  target  in  chaff.  The  s/c  is  shown 
relative  to  matched  right  circular  single  channel  operation.*  The  dipole  orientation  angles 
refer  to  the  distribution  within  the  chaff  cloud,  where  90  degrees  is  for  completely  random 
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Figure  >3.  Signil-lo-  Clutter  Impro  remen*  Ratio  for  Vector  Weighting  Method 


Manned  filter  is  standard  definition  -efened  to  nmiciirhjg  received  si n  ratio. 
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dipole  distribution,  and  15  degrees  refers  to  nearly  horizontal  distribution.  As  noted  in  the 
figure,  matched  filter  performance  for  single  channel  horizontal  and  circular  polarizations  is 
extremely  poor.  The  matched  vertical  performance  is  good  for  dipole  distributions  which  are 
nearly  horizontal  (less  than  30  degrees).  Beyond  this,  he  matched  vertical  performance 
degrades  rapidly,  and.  as  would  be  expected,  it  converges  to  the  circular  and  horizontal  levels 
for  the  completely  random  chaff  distribution  at  90  degrees.  This  convergence  at  90  degrees 
for  matched  filter  operation  was  also  evident  in  the  Stokes  resuits  for  single  pulse  transmission. 


Figure  5-3  also  indicates  four  variants  of  dual  channel  operation.  The  vertical 
transmit-optimum  receive  curve  is  for  a  dual  channel  receiver  optimally  weighted  for 
transmission  of  a  vertically  po'arized  five  pulse  sequence.  This  in  tact,  is  the  result  of  the 
first  iteration  associated  with  the  Rummler  waveform  convergence  technique.  Similar  curves 
are  shown  for  horizontal  and  circular  oolarized  five  pulse  sequences.  As  expected  inm  the 
predominantly  horizontal  chaff  dipole  distribution,  the  dual  channel  receiver  weighted  !*.r 
vertical  transmission  provides  the  best  overall  performance.  There  is  again  convergence  of 
results  at  90  degrees. 

The  optimum  curve  is  the  result  of  utilizing  the  Rummler  iterative  technique,  for 
which  both  the  transmit  vector  and  receiver  weight  are  interchanged  during  each  it-crabon. 

The  iterative  procedure  was  terminated  when  convergence  in  s  c  ratio  was  ach«ev*d.  Although 
the  optimum  curve  is  for  full  dual  channel  o;»eralion.  it  can  he  appreciated  from  prior  results 
that  the  optimum  performance  is  only  siightly  better  than  vertical  transmit-optimum  dual 
channel  receive.  This  is  a  consequence  of  the  optimum  transmit  wav  Torn  Kmg  nearly 
vertically  polarized.  A  significant  result  of  the  iterative  procedure  was  that  the  ronvergence 
in  transmit  wave-form,  was  indt  pendent  of  the  initial  choice.  The  iterations  were  started 
with  either  vertical,  circular,  or  horizontal  five  pulse  sequences.  In  all  cases,  the  optimum 
waveform  converged  to  the  identical  nearly  vertically  polarized  structure. 

5.4  MATRIX  APPROACH  RESULTS 

Computer  runs  were  made  for  probability  of  target  detection  Pr,  as  ;>  function  ot 
the  previously  defined  s  c  ratio,  at  a  faise  alarm  probability  of  10~T.  Performance  levels  for 
the  three  target  model  set  of  Figure  5-1  were  computed  for  chaff  dipole  distributions  corre¬ 
sponding  -o  15.  45.  and  90  degrees.  Comparisons  were  made  between  dual  channel 
operation  and  matched  filter  single  channel  operation.  The  optimum  transmit  vector  was 
chosen  from  a  set  of  1 14  three-pulse  sequences.  The  associated  dual  channel  receiver  weight 
corresponds  to  a  maximum  likelihood  detector,  which  is  a  function  ot  me  probability 
distributions  related  to  the  target,  chaff,  and  receiver  noise. 

5.4. !  lJud  Channel  Results 

Figure  5-4  is  a  comparison  of  dual  channel  performance  between  the  throe  target 
set  of  RTi.  RT3.  and  Sweriing  I.  The  limited  range  of  s<c  shown,  was  chosen  to  encompass 
the  P(j  region  of  interest  for  the  more  random  dipole  distributions:  namely  a  Pq  in  the  order 
of  0.5  for  0  =  45°snd  90?  At  dipole  distributions  near  the  horizontal  (0  -  15°).  performance 
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DIPOLE  ORIENTATION:  15°.  45°.  90° 


figure  5-4.  Comparison  of  Dual  Channel  Performance  for  RT3.  RTi.  SWt  Targets 


a«  these  s  -  ratios  is  v?r\  good  for  all  three  urget  representations.  This  would  he  expected  at 
this  dipole  distribution,  since  operation  v.’th  essential!)  *  rticalh  polar  zed  transmit  pulses 
yields  minima:  chaff  backscatter  (note  similar  good  results  at  6  =  15°  for  vector  and  Stokes 
approaches  to  using  vertical  transmission). 

As  !V  dipole  distribution  becomes  more  random,  performance  degrades  for  these 
s-'c  ratios  as  noted  in  the  tigure.  Performance  for  ail  target  models  at  a  Pp  ot  0  5.  for 
example,  is  about  5-6  dB  worse  at  0  =  90s  than  it  is  for  9  =  45°.  The  difference  in  detec¬ 
tion  performance  at  a  Pp  of  0.'  between  0  —  153  and  9  —  45°  is  no!  shown  in  Fisrure  5-4. 
The  basic  trend  of  the  curves  indicates  a  significantly  greater  difference  than  the  6  dB  noted 
in  the  prior  case.  A  rough  estimate  would  place  this  differential  in  the  order  of  10-12  dBu 
It  becomes  rather  obvious  from  these  and  prior  results,  that  the  structure  of  the  chaff  model 
emphasizing  horizontal  dipole  orientation,  strongly  influences  the  ultimate  choice  of  more 
nearly  vertical  transmit  polarization  vectors;  arid  ultimately  leads  to  poorer  performance  as 
the  distribution  becomes  more  random. 
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Regardless  of  the  chaff  distribution,  however,  target  detection  is  best  for  the  RT3 
model,  followed  by  RT1  and  SW1.  respectively.  This  trend  would  be  expected,  in  general,  in 
terms  of  target  detection  between  SW1  and  SWI1  models,  noting  that  the  RT3  model  tends 
toward  the  uncorrelated  puise-to-pulse  SW1I.  For  a  PD  of  0.5  and  a  PpA  of  10'4.  the  SWII 
requires  1  dB  less  signal  return  than  the  SWJ  in  the  case  of  a  three-pulse  sequence  operating 
onl>  against  receiver  noise.  As  noted  in  Figure  5-4.  the  RT3  target  requires  5  dB  Ilss  signal 
return  than  the  SW1  when  operating  in  the  dual  channel  chaff  plus  noise  environment  at 
0  =  45°.  This  added  enhancement  for  the  RT3  target  can  also  be  anticipated  when  utilizing 
the  polarization  discriminant.  An  analogy  can  be  made  if  one  examines  the  decorrelation  of 
echoes  from  a  complex  target  utilizing  a  polarization  agile  waveform  as  having  effects  similar 
to  frequency  agility  (Nathanson.  I969J  This  decorrelation  is  desirable  in  a  clutter  and/or 
noise-only  environment.  Averaging  of  independent  target  samples  will  then  essentially  reduce 
the  probability  of  a  null  in  the  target  cross  section  at  the  aspect  of  interest.  While  this 
averaging  process  does  not  reduce  the  total  clutter  power,  its  variance  will,  however,  be 
reduced.  Since  detectubilPv  is  proportional  to  the  ratio  of  target  power  to  the  variance  of 
the  tola*  clutter  plus  noise  power,  target  detection  v  11  be  improved.  This  decorrelation  of 
the  target  echo  suggests  a  movement  of  the  target  model  towards  the  more  desirable  SWII 
configuration. 

Although  complete  target  decorrelation  is  not  anticipated  for  the  three-pulse  sequence 
utilized,  it  is  noted  that  the  RT3  target  will  have  more  degrees  of  freedom  than  the  SWI  in 
terms  of  optimally  combining  the  target  echo  returns  relative  to  the  chaff.  It  should  be  noted 
that  the  SWT  target  will  also  become  par.:ally  decorrelatcd.  but  not  to  the  extent  of  either 
RT1  or  RT3.  This  then  results  in  the  enhancement  of  RT3  relative  to  both  RTI  and  SWJ. 
when  the  dual  channel  polarization  rcturtts  axe  optimally  combined. 

5.4.2  Matrix  Dual  is  Single  Channel  Results 

A  comparison  of  the  matrix  dual  channel  implementation  relative  to  conventional 
single  channel  matched  filter  systems  is  shown  for  the  three  target  models  in  Figures  5-5 
through  5-13  for  chaff  dipole  distributions  of  15.  45.  and  90  degrees. 

Figures  5-5.  5-6.  and  5-7  show  performance  comparisons  for  the  SWI.  RTI  and  RT3 
targets  respectively  for  9  =  15°.  The  horizontal  and  right  circular  polarized  systems  pert'erm 
very  poorly  in  the  matched  single  channel  configurations,  because  of  the  horizontal  distribu¬ 
tion  of  the  chaff  dipoles.  The  matched  vertical  system  performs  very  well  over  the  range  of 
s/c,  as  would  be  expected.  In  each  case  however,  the  dual  channel  implementation  prorides 
even  better  performance. 

Figures  5-8,  5-9,  and  5-10  show  similar  performance  comparisons  for  6  =  45s.  Single 
channel  horizontal  and  right  circular  systems  still  perform  poorly,  since  the  dipole  distribution 
remains  biased  toward  the  horizontal  plane.  All  three  dual  channel  systems  out  perfoim  the 
single  channel  vertical  systems.  At  a  PD  of  0.5.  for  example.  SWI  operation  is  2  dB  better, 

RTI  is  4  dB  better,  while  the  RT3  target  exhibits  9  dB  better  performance  in  the 
dual  channel  mode. 


Figure  5-9.  Dual  rs  Single  Channel  RT !  Target 
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Figure  5-11.  Dual  rs  Single  Channel 
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Figure  5-12.  Dual  rs  Single  Channel 


figure  5-/3.  Dual  vs  Single  Channel 


Figures  5-11.  5-12.  and  5-15  show  similar  results  for  0  =  90°.  With  the  dipoles 
now  randomly  distributed,  the  single  channel  matched  horizontal  and  vertical  systems  tend 
toward  the  same  performance  levels.  Here  again,  the  dual  channel  performance  remains 
superior  although  the  dB  differential  relative  to  single  channel  vertical  is  less  than  for 
6  -  45°.  This  is  to  be  expected  in  view  of  the  more  random  dipole  distribution,  which 
mitigates  against  the  exploitation  of  the  more  deterministic  chaff  polarization  at 
0  =  15°  and  45°. 


In  all  cases  of  dipole  distribution,  it  is  noted  that  more  effective  dual  channel  to 
conventional  single  channel  performance  is  passible  for  the  RT3  target  relative  to  either  the 
RT1  or  SWI  models.  This  again  illustrates  the  advantages  of  a  target  modci  which  is  not  as 
highly  correlated  on  a  pulse-to-pulse  basis,  parti-'ularly  when  exploiting  the  dual  polarization 
scattering  property. 

5.43  Transmit  Waveform  Selection 

The  transmit  waveform  dual  channel  vector  is  selected  by  a  minimization  procedure 
relative  to  the  Frobcnius  norm  (see  Section  3.1.2).  The  matrix  method  computation  relied 
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upon  a  library  of  prestoa-d  polarization  states  and  selected  a  3-pulse  sequence  which  most 
closely  matched  the  norm.  The  library  of  available  transmitter  polarization  states  is  illustrated 
in  Figure  5-14. 

This  polar  representation  maps  circular  polarization  at  the  origin.  linear  polarization 
along  the  circumference  of  a  circle  with  radius  =  1.  and  elliptical  polarizations  within  a  circle 
with  radius  <  1.  The  representation  is  such  that  the  direction  of  a  major  axis  of  the  ellipse 
for  each  polarization  state  is  along  the  radius  vector  which  goes  through  the  point  representing 
the  state  on  the  plot.  Similarly,  the  direction  of  linear  polarizations  are  along  the  radius 
vector.  Furthermore,  elliptical  polarizations,  whose  sense  of  rotation  is  right  handed,  arc 
trapped  on  the  right-half  plane  of  the  polar  plot,  and  left-handed  polarizations  are  mapped  on 
the  left-half  plane  of  the  plot.  The  cUiplicity  is  given  by  i-r.  where  r  is  the  radius  of  the 
circle.  In  this  way.  all  of  the  information  relating  to  the  polarization  state  becomes  readily 
available.  The  symbols  (A)  in  Figure  5-14  indicate  the  library'  of  prestored  polarization  states 
available  for  selection. 


Two  typical  waveform  choices  are  shown  in  Figure  5-14  for  the  SW1  and  RT3 
targets,  at  a  dipole  distribution  corresponding  to  0  =  45°.  The  RT3  sequence  consists  of 
three  pulses  all  linearly  polarized,  and  lying  from  11  to  22  degrees  from  the  vertical.  The 
first  pulse  is  at  22  degrees,  and  the  next  two  pulses  are  at  1 1  degrees  (note  small  double 
circle  in  figure).  The  SWI  set  consists  of  a  first  pulse  linearly  polarized  at  22  degrees  from 
the  vertical,  a  second  pulse  with  a  left-handed  ellipticity  of  0.56  and  the  major  axis  aligned 
21  degrees  from  the  vertical,  and  a  linearly  polarized  third  pulse  11  degrees  from  the  vertical. 
As  expected,  the  preferred  dual  channel  polarizations  do  lie  more  nearly  vertical.  The  SW! 
transmit  vce.jr  has  more  ellipticity  which  denotes  a  trend  towards  more  decorrelation  of  the 
SWI  model  through  the  polarization  discriminist. 


5.4.4  Receiver  Matrix 

The  receiver  weighing  matrix  is  computed  from  knowledge  of  the  transmit  waveform 
and  the  likelihood  function  related  to  the  target,  clutter,  and  noise  probability  distributions. 
Far  the  3-pulse  sequence,  the  optimum  receiver  matrix  GCpj  would  be  z  6  \  6.  The  SWI 
target  is  a  special  case  because  o*  the  invariance  of  its  pube-to-pulsc  scattering  matrix 
correlation  function,  which  yields  only  three  eigenvalues  in  the  maximization  procedure. 
Consequently.  G0pt  for  SWI  is  a  3  x  6  matrix.  A  typical  example  of  C.on,  for  SW!  Is  shown 
below  for  a  s/c  ratio  of  +3  dB.  and  a  dipole  0  of  45° 


Note  the  elements  of  G  are  written  in  polar  form 


Gopt 


(0.32.  -7*> 
(0.07,  67*) 
'0.003.  -72* 


(0.17.  -32’) 
(0.11.  -76°) 
(0.12.  -11°) 


(0.1.  -33* ) 
(0.05,  55’) 
(0.05.  -40*) 


(0.14.  -60*) 
(0.21.  -69°) 
(0.03.  -8.H*) 


(0.24.  -5*1 
(0.03.  -3o*i 
(0.02.  -31*) 


(0.1.  -45° > 
(0.04.  33*) 

(-14.  -5*). 
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Section  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  objectives  of  the  study  were  met  in  the  following  specific  areas: 

a)  An  analytical  and  systematic  appioach  vas  developed  for  dual  channel 
optimum  processing  applicable  to  the  radar  scattering  matrix  of  targets 
and  clutter. 

b)  The  matrix  approach  was  utilized  to  develop  an  optimal  transmit  waveform 
and  optimum  receiver  which  maximizes  the  signaTto-ciutier  ratio  for 
enhancement  of  target  detection  at  a  given  probability  of  false  alarm. 

c)  A  statistical  model  was  generated  for  the  polarization  scattering  matrix  of 
the  BQM-34A  target  drone,  and  for  a  dipole  chaff  cloud.  The  properties 
incorporated  in  both  models  focused  upon  the  temporal/range  behavior  of 
their  scattering  matrix  elements.  Two  forms  of  the  chaff  model  were 
developed;  a  theoretical  representation  composed  of  general  mathematical 
expressions,  and  a  stochastic  model  derived  from  an  autoregressive  process 
fit  to  the  theoretical  model. 

d)  The  above  models  were  utilized  in  computer  simulations  to  compare  target 
detection  performance  in  chaff  for  dual  channel  system  configurations 
against  conventional  single  channel  systems. 

Five  anaiytical  approaches  were  considered  for  polarization  processing  implementation; 


1) 

Matrix, 

2) 

Vector. 

3) 

State-space, 

4) 

Fredholm, 

5) 

Stokes. 

Of  the  above,  the  matrix  and  vector  approaches  fali  under  the  heading  of  fini't-  dimen¬ 
sional  space  techniques.  The  matrix  method  was  the  most  fully  developed  and  as  such  yielded 
the  most  promise  in  terms  of  meeting  the  study  objectives  for  dual  channel  waveform  and 
receiver  design.  It  utilizes  all  of  the  information  available  in  the  target/cluttcr  measurement 
space,  and  results  in  a  method  of  waveform  selection  based  upon  the  scattering  matrices  of 
the  target  and  clutter.  Once  the  waveform  has  been  chosen,  a  maximum  iikenhood  criterion 
is  utilized  for  target  detection  at  a  desired  false  alarm  rate. 

The  vector  method  can  be  considered  a  special  case  of  the  natrix  approach.  ;n  that 
only  a  part  of  the  available  information  is  processed  from  the  measurement  space.  Althou  h 
this  yields  better  results  than  conventional  single  channe’  systems  and  is  simpler  to  implement 
than  a  matrix  receiver,  its  performance  relative  to  the  matrix  receiver  is  obviously  not  ar  good. 

Considerable  effort  was  devoted  to  the  development  of  a  Kalman  filter  realization  of 
a  state-space  model  utilizing  the  target  and  clutter  scattering  properties.  This  approach  was 
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carried  to  the  point  of  an  optimum  receiver  design,  given  the  transmit  waveform  and  the 
relevant  scattering  matrices.  This  effort  exemplified  the  mest  consistently  unified  approach 
towards  utilizing  both  the  temporal  and  range  behavior  of  the  scattering  matrix  elements. 
Further  work  is  required,  however,  in  developing  this  approach  for  waveform  des*^. 

An  approach  '-‘ring  Fredholm  integral  equations  was  formulated  to  the  extent  of 
enumerating  ihe  system  of  equations  involved  with  the  target/clutter  scattering  elements.  A 
generalized  method  of  solution  was  presented.  No  further  effort  was  expanded  beyond  this 
point,  although  a  preliminary  conclusion  would  indicate  that  the  complexity  involved,  relative 
to  the  other  techniques,  might  not  warrant  aoditional  work. 

The  Stokes  approach  was  formulated  and  investigated  in  a  prior  study.  These 
results  were  reported  upon  for  completeness  and  indicated  the  capability  to  design  an  optimum 
transmit-receive  antenna  polarization  pair  for  maximizing  the  signal-to-clutter  ratio  for  use  with 
noncoherent  processing. 

Computer  derived  results  '-ere  presented  for  the  matrix,  vector,  and  Stokes 
approaches.  It  was  shown  that  for  either  single  pulse  or  coherent  train  processing, 
exploitation  of  the  polarization  discriminant  in  a  dual  channel  configuration  was  clearly 
superior  to  conventional  single  channel  matched  filter  performance.  The  matrix  results  also 
considered  the  differences  in  dual  channel  performance  for  three  variants  of  the  BQM-34A 
drone  point  target  model  in  chaff.  These  included  a  pulse-to-puisc  correlated  Swcrling  I 
model  and  two  less  correlateo  versions  tending  toward  Svverline  II.  Correlation  of  the  chaff 
scattering  matrix  was  set  identical  to  the  least  correlated  target  model.  RT3. 

Probability  of  target  detection  was  used  as  a  figurc-of-mcrit  for  a  false  alarm 
probability  of  10  ^.  Best  performance  was  obtained  for  the  least  correlated  target  model, 
with  worst  perform  ance  for  the  Swcrling  1  (a  typical  6  dB  differentia.  *br  a  Pp  of  0.5).  This 
is  consistent  with  tandard  results  for  Swcrling  11  relative  to  Swcrling  1  targets.  The  pulse-to- 
pulse  decorrelation,  inherent  in  .he  RT3  target  scattering  property,  provides  lor  more  degrees 
of  freedom  in  terms  of  optimizing  the  receiver  matrix  combination  associated  with  the 
individual  pulse  returns.  This  results  in  the  noted  better  performance  relative  jo  SWi.  It 
should  be  noted  however,  that  performance  of  the  Swcrling  !  target  was  also  enhanced  c-y 
dual  channel  polarization  processing  relative  to  single  channel  matched  filler  det**ction. 

The  major  conclusions  reached  from  the  foregoing  results  can  be  restated 
as  follows: 

a)  Exploitation  of  polarization  scattering  properties  will  enhance  the  detection 
of  targets  in  chaff  relative  to  conventional  single  channel  svstems. 

b)  Analytical  and  systematic  approaches  have  been  presented  to  develop 
optimal  transmit  waveforms  and  associated  receiver  designs  for  maximizing 
the  probability  of  target  detection  for  a  given  false  alarm  rate. 

c)  The  availability  of  previously  developed  polarization  sensitive  models  of 
targets  and  chaff,  provides  a  powerful  tool  in  the  systematic  evaluation 
of  the  above  techniques. 
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The  following  study  tasl  s  are  recommended  to  further  evaluate  the  practicality  of 
multichannel  polarization  processing: 

a)  Expansion  of  the  state-space  approach  to  include  waveform  design  as  a 
possible  adjunct  to  or  replacement  and  generalization  of  the  matrix  method. 

b)  Reconfigure  the  target  and  chaff  models  to  include  realistic  target 
trajectories  and  chaff  mean  velocities,  to  combine  MT1  processing’  in 
conjunction  with  the  polarization  discriminant  for  dual  channel  waveform 
and  receive:'  design  evaluation,  b  e  of  longer  pulse  trains  would  be 
investigated  for  this  implementatioi 

c)  Develop  and;or  utilize  other  ciuttc.  models  to  evaluate  performance  in 
different  interfering  environments.  EUipticitv  associated  with  rain  drops, 
for  example,  would  be  a  reasonable  contrast  to  the  dipole  structure  of 
the  chaff  model. 

d)  Develop  a  strategy  for  detection  of  targets  in  clutter,  without  prior 
knowledge  of  target  scattering  matrix.  Use  typical  surveillance  radar 
requirements.  Evaluate  performance  of  various  approaches:  for  example, 
use  waveform  to  minimize  clutter  backscatter  or  maximize  s'c  ratio  for 
‘average'  target. 

c)  Preliminary  design  of  dual  channel  system  utilizing  methods  developed 
during  study.  Provide  qualitative  cost  assessment  against  single 
channel  system. 

f>  Develop  an  approach  using  the  range  spread  properties  of  the  scattering 
matrix  in  conjunction  with  its  polarization  and  temporal  behavior,  to 
design  the  waveform  and  receiver  for  optimizing  target  detection  in  clutter. 
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Abstract 

.•?  ret  h‘r  mitt -regressive  prut  css  is  used  lu  obtain  a 
Slot  hast  it  model  lor  the  simulation  of  both  the  v pcttral 
and  fmlarizaiion  characteristics  of  a  chaff  cloud  consisting 
of  a  collection  of  rotating  dtfmles.  whit  I:  mav  hai  t  either 
t  uni  f 'letch  random  or  preterre!  orientations. 


INTRODUCTION 


Recently  there  has  been  considerable  interest  in  the  use  of  polarization  diversity  to  improve 
signal  detection  in  the  presence  of  background  clutter/  In  the  design  of  such  polarization 

agile  radar  systems,  one  must  consider  both  the  autocorrelation  and  cross  correlation  functions  of  the 
baokscattercd  signals  for  two  orthogonal  polarizations/*^ 

What  follows  is  a  model  for  a  random  collection  of  rotating  dipole  stutterers.  This  model 
is  a  generalization  of  a  model  by  Wong  e!  al/4)  and  allows  the  simulation  of  both  the  polarization 
and  temporal  characteristics  of  radar  signals  backseattcred  from  a  collection  ol  rotating  dipoles,  which 
have  either  completely  random  or  preferred  orientations.  This  model  makes  use  ol  a  vector  auto¬ 
regressive  process  for  the  simulation  of  the  scattering  matrix  elements  of  a  dipole  cloud  and  can  be 
used  for  the  statistical  simulation  of  signals  scattered  from  chaff. 


THE  SCATTERING  MATRIX  OF  A  ROTATING  DIPOLE 

When  a  dipole  is  rotating  about  the  s-axis  of  an  orthogonal  coordinate  system,  (s,  u,  v)  as 
shown  in  Figure  A- 1,  the  direction  of  its  dipole  moment  is  given  by 


A  A  ,  A  . 

d  -  u  cos  ^  +  v  sin  v 


(1) 


where  i//  is  the  angle  ol  the  dipole  with  the  ti-axis  given  by 


=  Wj-t  +  Ol. 


(2) 


wr  is  the  dipole  rotation  rate  in  rad  .'sec 


and  a  is  a  constant  phone  angle. 

If  the  pol.u  angles  of  the  above  s-axis.  relative  to  a  Cartesian  frame  (x,  y,  z),  are  0  and  0, 
as  shown  in  Figure  A-l.  and  if  the  u-uxis  is  chosen  to  lay  in  the  x-y  plane,  then  by  a  coordinate 
rotation  from  the  (s.  u,  vi  to  the  (x,  \.  /i  system  we  obtain  that 

d  ••*  -x  tsin  <}>  cos  y <  r  cos  0  cos  <>  sin  vl't 


4  y  (ios  <£  cos  i p  -  cos  0  sin  <*>  sin  \fr) 

+  /  sin  0  sin  ^ 

When  an  electro  maguetii  wave  is  propagating  along  the  z.  axis,  the  scattering  matrix  ol 
the  above  dipole  is  given  by  (Borisoi/  *'*.  Wong  el  al<4)) 
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_  /hxx  ^xy\ 
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where  (X  is  the  radar  wavelength). 


From  equations  <3)  and  14)  we  find  tlut  the  elements  of  the  scattenng  matrix  H  can  be 
expressed  as 


and 


hxx  =  X  IC,  +  D|  c*2*  +  D,*  c~}2V i 


hyx  =  hxy  =  vr°l)  X  lC2  +  *  t>-** 


hvy  =  v7  0.9  X  1C,  +  I),  J2*  *  D,*  e"*2*I 
where  the  superscript  *  de  soles  complex  conjugate. 

]  ‘’I  “)  **  . 

C,  =  —  lsin~  <>  ^  cos-0  cos-<*l 

1  S  s  s  s 

I) j  =  j  (sin-  0  -  cos-  0  cos-  Oi  -  j  sin-  (■  cos  0)  . 

i  s 

(i  =  -v  -an  0  vO>  o  *.*n~  0  . 

5  s 

! )  •»  :  I  sin  0  cos  0  *  sin  0  c*»s  0  cos-  '}  -  j  j 'o>  cos  ®1 

|  '  s  s 

(  ;  :  ,  !l'm"  o  *  Mil”  0  cos”  til  . 

I  "%  ^  X 

1),  -  —  Iscift-  o  -  sin*  o  cos-  fit  t  J  sin  i20>  uss  0\ 


(5) 


(6i 


(7) 


STATISTICS  OF  hx%.  hyy  FOR  A  RANDOM  COLLECTION  OF  ROTATING  DIPOLES 

A  study  of  the  temporal  and  polarization  response  ot  a  random  collection  ot  rotating 
dipoles  requires  estimates  of  autocorrcbtioas  and  cross  correbtion  functions  for  Jl:-e  elements  of  the 
scattering  matrix  H  of  Equation  «4).  such  as 

h!!ixxtt»  hxx*tt-ri|  .  Ejhxxtt»  i»vy*«-T»|  . 

Eliixx(t>  hxy**t-r»!  .  FJh  tti  hxy*4t-T)!  .etc. 

where  the  expectation  is  over  time  and  all  possible  orientations  for  the  axes  of  rotation  of 
the  dipoles. 

These  .ebtions  can  he  expressed  in  a  compact  torm  by  use  of  the  covariance  matrix 
K||tn  =  l-.lhtn  h*  (t-m  « 8> 

of  a  vector  hf tl.  defined  as. 


hit)  = 


if  the  distribution  function  for  the  an-ubr  rotation  rate  is  pw  twfi  art'  for  tin-  Doppler 

shut  due  to  the  drift  velocity  is  p„  Iwji.  t:xn.  unne  tlie  results  of  Hone,  et  al.*'**  and 

"J  3 

substituting  iron!  Equations  t2).  <5>.  tb)  and  (~»  into  liquations  (8)  and  <">).  »c  obtain  tlut  Kjjtr) 
can  be  expressed  as 

K„!t>  =  M.CtI  |C  ♦  2 Re  «M„  t2r»  *  D»i  <10> 


where  the  matrix  <  is  erven  b\ 


F-C jl~  IfCjCs  i  HC,C? 
H<2rj*»  l*si-  ItC,Cs*|  | 

Etc/}*)  ll<3r2*;  !;</-  j 
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*he  matrix  D  is  given  by 
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and  hit.  (( .  i  .  I-IDjD:  |  arc  sj*herkal  random  averages  of  the  orientation  of  the  rotation  aces 
(figure  A- 1 1  given  by 


K/'i 


M 


do 

-'0,  .-B 


<?r 


sin  1/  tiO 


(13) 


i°:  |fl: 

I  do  }  'ii!  0  dr*  =  (o->  -  0|»  uos  ttj  -  ci*s  git 
Oj  -•*, 


ixz  0}  «  O'  <  2x  .  fl,  ^  e  ^  r 


When  the  dipcdcs  have  all  possihk  orientations  and  the  variables  ttj  and  wf  arc  no-mal. 
with  densities  N(ij.  °,}"t  and  Ni»;  °r~>-  respectively.  »v  obtain  that  the  covariance  matrix 
ts  psvs  t»> 


Kj|(t)  -  exp 


r  2’!  /!  0  i\  r  ,  .  A 

[*«dT  - -4—  0  £  c  \  ♦  exp  J-IWjT  -  —  j  «*  (2WT) 


STOCHASTIC  MODELING  OF  THE  SCATTERING  MATRIX  ELEMENTS  OF  CHAFF 

If  we  assume  that  the  vector  h  (f >  of  Equation  (8)  is  a  sample  from  a  Gaussian  vector 
process,  then  use  of  a  vector  autoregressive  process  yields  a  rvlativelv  simple  linear  model  that  has 
h  (t)  as  its  output.  In  general,  a  vector  autoregressive  process  requires  the  use  of  -iatn\ 
coefficients'6*.  However,  because  of  the  stricture  of  the  covariance  function  of  h  <ti.  given  in 
Equation  *  10;.  the  stocliastie  vector  h  (tl  can  be  modeled  as  the  sum  of  three  ve-.Iui  autoregressive 
processcs  h  jiU  li  sit)  and  h  jit)  with  scalar  coefficients  as 


htnAt.  =  ^  h  |(nAt)  ♦  hdndt)  «p  f-j^iyitn) 

♦  h  jin  At)  exp  t-»j2wrAln)  ^  exp  i-jw-jAtn) 


where 


M. 


hjinAt)  =  2  2-hj  (in-i»At)  +  X„ 
Vs 

-  t*  h>  (in-iiAt)  •»  Y 

«-I  *  ' 


hst.iAn  = 


Mi 


h  nAi '  i 


b,  *\;  *  *n 


1 151 


ji:.I  the  jj's  and  are  scalar  constants  obtained  front  the  Yule-Walker  1  cputntns  if-. I 


i 


L  a,R |  [il-:lAl] 


R;ikAt«  .  L 


1  M,  . 
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’O'  , 


iW/Ma1  illtJXus&tfix 


and 

Mi 

E  ‘  bjR->  [(k-i)Atj  =  R-dkAt)  .  k  =  1  .  ...  Ms 
i=l 

which  can  be  solved  by  a  fast  recursive  procedure  developed  by  Durbin  '  . 

The  orders  of  autoregression  Mj  and  M>  depend  on  the  sampling  interval  At  and  tlic 
desireu  accuracy  in  approximating  the  correlation  functions  RjtkAt)  and  RstkAt).  where 

RiikAt)  *  Mj(kAt)  evaluated  for  wj  =  0 


RsfkAtt  =  MjlkAt)  *  Mw  (2k At)  evaluated  for  wj  “  0  and  wf  =  0 
The  vectors  Xn  .  Yn  and  Zn  arc  independent  noise  vectors  with  sutocovariancc  matrices 

fc«X|Xj*l  «  if 
HY.Y/i  =  4,j0 


and 


HZ/-,  *  “  *,!*> 

where  t!u  i>ar  denotes  «.nsopk-\  conjugate. 


In  numeruai  simulations  ’isesc  .an  be  generated  Iroin 


where 


^0  c0  =  C 


D0  O0  *  D 


The  asterisk  denotes  complex  conjugate  liansposc  and  the  vecloi  L(  V  and  ft,  are  Ihroe- 
•iimenaonal  complex  white  noise  sectors  with  Gaussian  uncorrelated  components  In  mans  ca*es 

-v 

tlie  covariance  matrix  l)  is  real  or  wf  =  0.  If  either  of  these  conditions  is  true  then  one  can 

snow  that  htoAtt  can  lie  expressed  as  the  sum  of  only  two  autoregressions.  which  reduces  the 
computations  needed  for  simulation  of  the  elements  of  the  scattering  matrix 

NUMERICAL  RESULTS 

As  an  example  of  the  results  obtained  by  use  of  the  above  model,  we  show  tsi  figure  A-2. 
a  pobr  plot  of  the  polan/c‘!ion  states  of  waves  hackscattered  from  a  random  dipole  cloud  for  three 
different  transmitter  antemu  pobn/ations  vertical,  horizontal  and  circular  I>ata  tor  these  plots 
were  obtained  from  the  above  model  for  utj»o!cs  oriented  primarily  within  j»0  degrees  irom  the 
horizontal.  The  pobr  representation  (r.  <>i  of  the  polarization  states  in  these  ligurcs  is  such  that 
the  length  of  the  radius  xcctoi  is  r  =  1  -  c.  where  e  is  the  elhpticity.  i.c...  circubr  polarization  is 
mapped  at  the  origin,  linear  polarization  along  the  circumference  of  a  circle  with  radio-  =  1.  and 
elliptical  pobrizations  within  a  circle  of  radius  =  1.  The  polar  a->glc-  <>  is  such,  that  direction  of 
the  major  axis  of  the  ellipse  for  each  pobn/aiion  state  is  along  the  radius  vector  which  goes  through 
the  point  representing  the  state  of  the  plot.  Ssmibrlv.  the  direction  of  linear  polarizations  are  alone 
the  radius  vector,  furthermore,  elliptical  pobri/atsor.s  whose  sense  of  rotation  rv  ngbt-handcd.  arc- 
mapped  on  the  right-half  pbne  of  the  pobi  plot.  Ihu>.  tlie  above  graphic  representation  is  capable 
u*  representing  all  polarization  states  and.  in  eflect.  maps  the  -urfacc  of  tfe  Poincare'  pobn/alion 
sphere  onto  a  circle.  Tlie  only  poun/atior.  state  which  is  ambiguous  in  this  representation  is 
circubr  pob'-./atron.  Both  right-handed  and  lelt-fiandcd  circular  pobrizations  are  n» »pped  on  the 
•ame  point  at  tlie  origin  of  the  coordinate  system  at  the  center  of  the  polarization  circle.  Tor 
elliptical  polarizations  whose  major  axis  is  along  She  vertical,  we  have-  aruitrarilv  chesses  to  map 
left-handed  pobn/ations  onto  tin:  lower  halt  ot  the  vertical  axis,  and  nght-handed  pobn/alions  onto 
the  apj>c:  half  of  the  vertical  axis,  f urthermotc.  tlie  size  of  the  letter  used  Sot  indicating  tfie 
pobn/ation  of  the  received  vv  ivc  is  varied  m  accordance  with  the  total  signal  power  in  the  w^ve. 

IX..  maximum  rower  ;hal  can  be  extracted  fV«»m  the  wave  by  a  matched  antenna. 

FORMULAS  FOR  THE  EXPECTATIONS  EiCjC^t  AND  ElOjDj*} 
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/•'foim*  A-2.  Ilavkscattvml  Wave  Polarizations  from  dipoles  with  orientation 
r  ,io°  from  liorl:atlonal 
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Appendix  B 

LIKELIHOOD  DETECTION 
BINARY  HYPOTHESIS  TESTING 

In  order  to  make  a  decision  whether  a  target  is  present  in  chaff,  binary  hypothesis 
testing  is  performed  to  determine  its  decision  region  [see  Van  Trees  1968].  Referring  to 
Figure  B-l,  the  two  hypothesis  Hq  and  H  j  with  respect  to  the  received  signal  Y  (with 
additive  narrowband  white  noise)  are: 


Hypothesis  Hq: 
Hypothesis  H|: 


y  =  Y.  +  n 


Y  =  YT  +  Yr  n. 


The  input  signal  Y  is  then  processed  through  the  optimally  weighted  matrix  G  whose  output 
Z  is  given  by: 


Z  =  G  •  y 

Hence, 

Z|Hq  =  G  -  (Yc  +  N0) 

Z1H,  =  G  •  (Yx  +  Yc  +  N0) 

and  the  respective  covariance  matrices  are 

CovCZIHq)  =  G  ♦  E*vcYc*)G*  +  G  -  E[Nq  Nc*l  •  G*  =  Rq* 

Coy(Z1Hj)  =  G  -  E{YTYT*JG*  +  Rq^  *  Rx  + 

In  genera],  the  covariance  matrices  Rj  and  R^  are  not  diagonal,  and  it  is  desirable 
to  perform  a  coordinate  transformation  such  that  the  resultant  covariance  matrices  are 
simultaneously  diagonalized.  The  new  set  of  transformed  random  variables  are  now  statistically 
indepenoenL  In  order  to  find  the  desired  coordinate  transformation  we  solve  the  generalized 
eigenvalue,  eigenvector  problem:  Find  the  eigenvalues  /ij  and  eigenvectors  ^  such  that 

rT  $i  =  #*i  RCN  i  =  1,  .  .  .  k 


«s  satisfied. 


%  v-w'  sw 


Further,  we  define  another  transformation  P  such  that 
r  =  P-Z  =  P  -Q  *  Z  =  P  *  Q  *  G  •  Y 


P  =  [I-(ji+lHl^  =  k/l+jij 


A(r)  =  f*  r  +  const 


and  the  decision  criteria  is 


u  =  r*  r 


=  2  i  r- !~  >  T  H,:  Target 


present 


Consequently,  the  probability  of  false  alarm  Pj.-^  is  siren  by 


PFA  =  Jj  P<uiH0>d- 


and  probability  of  detection  b  given  by 


*D  ~  f°  F(u|H|  >du 


and  the  optimal  receiver  weighting  matrix  b  given  by 


Gopt  =  P*<?*G 


■  _  ‘  \ 


WM4U 


THRESHOLD  DETERMINATION 


The  general  form  of  the  probability  density  function  p(u)  is  given  by 
k 


p<u>  =£  * 


i 


-  u/Oj2 


k  j  or~  \  of 

-sr*) 

j^i  \  / 


:  u 


where 


u  = 


V 

Z-r 

S—l 


«T  =  f'irji- 


p(rj.  ....  rk)  =  II  P<ri> 


,*  k  ^  -  T/o-~ 

P(T)  =  JT.  p(  aj  =2rf 


i=l 


n  (J-Oi'/oj2) 
j=l  -  1 


For  a  given  PpA  the  threshold  T  can  be  computed  from  the  above  equation  using  the 
Newton-Raphson  iterative  technique  with  the  appropriate  op  for  condition  and 

P(T)  =  PFA- 


Appendix  C 


CONSTRAINED  OPTIMIZATION  OF  S/C 

The  constrained  optimization  problem  is  as  follows.  Given  the  target  and  clutter 
convariance  matrices  Kyj  and  Kyc  fhc  optimum  W  that  maximizes  the  signal  to 
clutter  ratio  defined  by 


S/C  A 


W^KyjW 

w*kycw 


(1) 


where  W  is  a  N  xl  complex  column  vector, 

*  denotes  complex  conjugate  transpose 
and  Kyp,  Kyp  are  both  positive  definite  hermitian  matrices. 

This  problem  is  equivalent  to  maximizatior.  of 

W*KYTW  (2) 

subject  to  the  constraint 

WKYCW  =  1  (3) 

Let  X  be  the  Lagrange  multiplier  and  define  F(w)  as 

Hw)  =  WKyjW  -  X  (W*KycW  -  n  (4) 

In  order  tc  use  the  Lagrange  multiplier  for  maximizing  F(w).  both  real  scalar 
functions  W*KyjW  and  W*KypW  must  be  difTerentiable  with  respect  to  W  which  is  complex. 
Applying  the  Cauchy-Riemann  Equations  to  test  for  differentiability,  the  test  fails  and  it  is 
eenduded  that  the  scalar  functions  are  nowhere  difTerentiable  in  CN  (complex  N  dimensional 
space).  An  alternate  approach  is  to  solve  the  differentiation  problems  in  real  space 


Let  W  =  X  +  JY 

(5a) 

Kyr  -  Re  iRyyJ  +  j  Im  |Kyy] 

(5b) 

Kyc  =  Re  {Kyc1  +  j  lm  |Kyc| 

(5c) 

C-l 


ISSSMB 


S**'****’  ♦’-ifcil  -T  *Z'*Sk  Jscvwi 


For  any  hcrmltian  matrix  M  we  have 


M  =  Uel-Mj  +  j  I_IM1 


MM]  =  (Rel.M])T.  lmlM|  =  — <Imf MJ ) s 


where  superscript  T  denotes  tranpose. 


Therefore.  W*MW  can  be  rewritten  as 


W* MV  =  XTRe|M]X  +  YTRe{M]Y  +  YTImtM]X  -XTIm|M]Y 

+j{ XTRe] M 1 Y  -  YTRe]MlX  +  XTlmlM]X  +  YTIm]M]Y) 
=  XTRefM;X  +  YT*Re[M}Y  +  YTImlM]X  -  XTIm]M]Y 


=  RejW’MW] 


Note  that 


XTnetM]Y  =  (XTRetM]Y)T  =  <Re]MJY)T(XT)T  -  YTRe(Ml  X 


-XJIU 

XT,m]M]  X  =  «XTImlMlX)T  =  «!miM|X)T  (XT)T  =  XTt-ImJM] >X  =  0 

W*MW  can  now  be  expressed  as  follows 

W*MW  =  <XT.>T»  TrcIM]  -  1_{M]"I  /X  \  =  fXT.YT)  M,  /x\  «7| 

\V )  \Y  / 

where  Mj  is  a  positive  definite  symmetric  real  matrix. 

New  we  introduce  the  gradient  operator  7^-  defined  as 
„  .  .1  3  3  3  3  3 


3  3 


V«eA  i —  -  - - - 

3Xj  3X-,  '  dXN*  3Yj  3YV 


man 

ItuA  *******  W* V  f  |  ${] 


m 


'»  f  t  ...  A# 


Wa* 


That  is,  from  Equations  (7)  and  (!0) 


VWF(W)  =  ^^I^KyyW  -  X  (WKYCW  -  1)1 


iKYci||  M  -  0 


•  -  UKYtS|  pel^vd  •  -  W^vdl  A 

•mSKyrl  fe(K„l  j  "  [yKYC!.  feJKYClj  \Y 


Equation  (12)  is  equivalent  to  two  simultaneous  equations  set  to  zero 


(Rc(KYTj  -  XRe[KYCI)X  -  (l^Kyy!  -  XSm|KycJ)  Y  =  0 
Ojj-JKyjI  -  Mn-IKycl)  X  ♦  (RefKyjl  -  XRe}Kyc])  Y  -  0 
Therefore,  we  can  express  the  two  simultaneous  equations  into  one  as  follows 
(RelK^j  -  XfeiKYCJ)X  ♦  Kl^fKyyl  -  >^1^}  K)Y) 


-  XlmiXYrl»  X  *  (RelKyrl  -  XRe-[Kyc}  )Yj  -  0 


a  equivalently 


KRefK^j  -  ARelKYrl>  ♦  j  (yiCyjl  -  X  y^M  •  JX  +  iY) 
”  ^YT  *  ^  Kycl  W  =  0 


(13a) 


<  13b) 


VWF(w)  -  2  (Ky-J  -  X  Ky(-)  *  W  -  0 

The  solution  to  (Kyy  -  X  Ky^)  *  W  =  0  can  be  obtained  using  the  simultaneous 
diagonaiization  of  two  hermitian  matrices  discussed  in  Appendix  D. 


It  then  follows  that  there  are  N  such  extreme  values  Xj  and  transformation 
vectors  W-  (equals  bj  in  Appendix  D)  that  satisfies  Equation  (16)  as  given  below 

Wi*  KYT  Wj  =  Xj 


*’i*  Kyc  Wj  =  I 


Substituting  Equation  (17)  into  Equation  (1)  yields 


(S/C);  = 


Wi**YTWi 

Wi*KYCWi 


and 


^ax(S/Q  =  X^ 

Hence,  by  selecting  the  largest  eigenvalue  X^=  Xma_)  and  the  corresponding 
transformation  vector  W-{=  WQpt).  then  (S/C)  is  maximized  over  all  W  satisfying  the 
constraint. 


Appendix  D 

SIMULTANEOUS  DIAGONAUZATION  CF  TWO  COVARIANCE  MATRICES 


Let  A  and  B  be  the  two  complex  NXN  covariance  (positive  definite  and  hermitian) 

matrices. 

We  show  that  the  two  covariance  matrices  A  and  B  can  be  simultaneously  reduced 
to  a  diagonal  form  by  a  NXN  matrix  d  such  that, 

d*Bo  =  A 

d*Ad  =  I  (1) 

where  A  is  a  NXN  diagonal  matrix  whose  elements  are  real  and  I  is  a  NXN  identity  matrix. 
The  columns  of  d  are  the  eigenvectors  dj  satisfying 

(B  -  AjA)dj  =  0  (2) 

where  the  Aj  are  the  eigenvalues  of  the  matrix  equation 

jB  -  XA|  =  0  (3) 

where  1  *  !  denotes  determinant  of  a  matrix. 

Simultaneous  dia canalization  of  the  matrices  are  performed  in  two  steps.  The  first 


aiep  is  to  diagonalize  matrix  A,  that  is. 

W*AW  =  SI  (4) 

Tnc  columns  of  W'  are  the  eigenvectors  wj  satisfying 

A*Wj  =  UjWj  i=l,...,N  (5) 

where  ouj  are  elements  of  the  diagonal  matrix  SI  and  the  eigenvalues  of  the  matrix  equation 
!  A  -  cul  |  =  0  (6: 

Note  that  A  is  hennitian  and  thus  SI  is  real  and  W  is  a  unitary  matrix,  i.c.. 

W*W  =  vHw  =1  (?) 

From  Equation  (4). 

A  =  WS2W*  (S) 

D-l 


Since  A  is  hermitian,  the  eigenvalues  cjj  are  real  and  thus, 


£2  =  £2'*  *  £2’-  (9) 

Now  A  can  be  rewritten  as 

A  =  W  -  £2*-  *  £2--  -  W*  (10) 


and 


£2-K-  -  W*  -  A  *  W  •  £2~-  =  !  (11) 

Substituting  Hqu.»t:on  (10)  into  liquation  (3)  yields 
lB-XAi  =  iB  -  XfW£21-*£2,%»)i 

=  i  W£2  •*  (£2'’%*B  -  W£2_!-  -  X!)  £2%*  j 

=  i  W£2  - 1  -  i  frWB  *  Wn_!-  -  XI  |  •  i  £2*%*  !  (12) 

=  0 

Both  !W£2l-i  and  i  |  are  non -zero  hence 

j  £2-!-'\V*BW£2-1-  -  XI  1  =  0  (13) 


Let 


K  =  £2~’I\V*  ♦  B  *  W£2_,i  (14) 

Note  that  K  is  abo  hermitian. 

The  second  step  is  to  diagonalize  K.  The  same  procedure  3S  the  dia  socialization  of 


matnx  A  >s  used. 

V*KY  =  A  (15) 

where  the  columns  of  the  unitary  matrix  V  are  the  Msrnvectors  v;  satisfvine 

Kvj  *  XjVj  :  =  i . N  (;6) 

and  Xj  are  the  eigenvalues  of  the  matrix  equation 

|K  -  XI  |  =  0  (17) 


Substituting  K  from  liquation  114)  into  liquation  (15).  we  have, 

V*S2"!HV*BW  -  n_!-V  =  A  (IS) 


Premult. ply  and  postmultiply  Equation  (11)  by  V*  and  V.  respectively. 

v*sr,2\v*A\vir!-v  =  v*iv  =  \*v  =  i  (19) 

The  last  equality  in  liquation  t!9)  is  due  to  the  property  of  a  unitar.  matrix. 
Defining  th-*  matrix  operator  p  as  follows 

o  =  WST*=V  (20) 

Then  liquations  (IS)  and  (19)  yield 
o*Bp  =  A 

0*Ao  =1  (21) 


Appendix  E 


BHATTACHARYYA  DISTANCE  B  FOR  GAUSSIAN  DENSITIES 

It  is  assumed  that  the  conditional  probability  densities  under  the  two  hypotheses 
HQ  and  HI  are  zero  mean  Gaussian.  It  is  noted  that  the  random  variable  Z  is  complex  and 
the  real  and  imaginary  components  are  independently  zero  mean  Gaussian.  Hence,  the 
waveform  is  described  in  terms  of  Rayleigh  distributed  amplitude  (envelope)  and  uniformly 
distributed  phase. 

The  Gaussian  statistical  representation  is  a  reasonable  one  for  chaff  as  well  as  for 
targets  of  the  Sweriing  1  and  II  class  where  the  target  amplitude  iluctuation  is  described  by 
a  Rayleigh  distribution. 

Let  the  two  conditional  probability  densities  pfZjHQ)  and  p(Z|Hl)  be  zero  mean 
Gaussian  densities  under  toe  two  hypotheses  HO  and  HI.  The  zero  mean  Gaussian  densities 
densities  arc  { Helstrom  1968] 


where  K^q,  ^Z1  arc  ***  covariance  matrices  of  HO  and  HI. 

Note  that  |*i  denotes  the  absolute  value  of  the  determinant  of  a  matrix. 


Substituting  Equations  (1)  and  (2)  into  the  integrand  of  Equation  (3-34)  in 
Section  3,  the  integrand  becomes 


E-i 


r-  <W{^  „' 


and  evaluation  of  J  using  Equation  (3)  yields 


Jf 


*vi 


,,  -1  .  .  -!j  -1 

kZ0  +  *Z1  I 


=  /l(Z)dZ  = 


IKzali4-|Kzll'4 


r _ i 

^  szo-1  + 

a  (2»)k - 


r  *1  a  v  “1  —1 

kZ0  *  KZ1 


/Hzo"1  "kzi_1\ 

-  exp  -ViZ*E - :; - jZ  *  < 


czo_1  +  Kzi  1 1  1 


Kzo1*  •  ikzi{1A 


The  term  in  the  numerator  of  Equation  <4)  (KZq  *  +KZj  *)  030  ^  rewritten  as 


<KZ0',+KZr‘)  -  Kze*‘ (I  +  KZ0KZ,-1)  =  Xzo  1  (KZ1  +  KZ0)KZ1 
Then,  the  determinant  of  (KZg~!  +  KZj-1>  becomes 


ir _ -1 


IKZO  1  +KZ1  *•  =  »ZXTl  (K7!  +  ‘:zo)Kzi  ~  iKzor1  *  i^Zl  +  KZ0!  ’  *KZl' 


=  !k% 


I 


Hence,  substituting  .quation  (5)  into  Equation  (4)  yields 


Kzo  +  Kzir' 


IKZ0IR  •  IKZII% 


The  Bhattacharyya  distance  B  is  given  by 


B  =  -fci  J  =  -Sn 


m".  •  IKz.i*] 


[~  IK/.0  +  KZ 1 1  1 

P  IKzol’5  •  IKZ1IW 


as  shown  in  Equation  (3-35). 


E-2 


Equation  (7)  can  be  further  rewritten  as 


lKZ0  +  KZli“ 

B  =  rifin  — r— tv — r  +  k/2  fin  V* 


iKZ0!  *  IKZjl 


In  Equation  (3-36)  KZj  and  KZq  are  given  by 
KZ1  =  w  ^-yT  W*  +  WK-C  w* 


KZ0  =  WKyc  W* 


(8) 


(9) 

(10) 


Substituting  KZq  =  W  KyC  W*  for  A  and  W  Kyj  W*  for  B  in  Appendix  D,  it  is 
shown  that  there  exists  a  coordinate  transformation  matrix  4  such  that 


4*64  =  o* 

(WKyTW*)p  =  A 

(ID 

4**4  =  d»(WKyCW«)d  =  i 

(12) 

Let  Z  =  0z 

With  this  coordinate  transformation,  we  get 

Z*KZ1Z  = 

Z*W  K>t  W*Z  =  z*  d*W  KyT  W*  dz  =  z*  -Vz 

(13) 

Z*KZ0Z  = 

Z*W  Kyc  W*Z  =  z*  6*W  Kyc  W*  pz  -  z*  Iz 

(14) 

Substituting  Equations  (13)  and  '14)  into  Equations  .1)  and  (2), 

Kzi  =  A  +  i 

OS) 

*z0  =  I 

(16) 

Substituting  Equations  (IS)  and  (16)  into  Equation  (18)  yields 

B  =  V&n 

r  ->] 

KA  +  I)  +  ir 

A  +  !  j 

=  *£n 

IA  +  2»J“ 

|A  +  I[  +  k/2  in  V* 

(t?) 

E-3 


A  is  a  diagonal  matrix  whose  elements  are  the  eigenvalues  Xi,  i  =  l,  .  ,  k. 
Therefore  the  determinants  of  (A  +  2?)  and  (A  -i-  !)  can  be  expressed  in  the  produce  torrn 
2s  follows 


•A  +  211  =•  n  (Xj  +  2) 
i=  1 


1A  +  1!  =  n  (Xj  +  1) 

i=i 


Finally  we  obtain  the  desired  result 


B  =  ViSnl 


k  >1 
n  (\  +  2»- 

— - +  k/2  bi  '<* 


n  (X,  +  1> 

.*“»  ‘  J 


k 

=  £  -  Cn 
i=l 


—  — 
t.Xj  +  2>~ 

p-TT,  +  €n  Vi 


k  r  i  i 

=  Vh  s  Cn  Oj  *  1)  *  jjpj-j  +  2  >  k/2  ^  ra 


equation  (20)  is  the  result  shown  in  Equation  (3-27)  in  Section  ?. 


Ref:  Helstrom,  C.  «V_  '‘Statistical  Theory  oi  Signal  Detection,”  Appendix  D.  Perea  nun  Press, 
1968 


Appendix  F 

AUTOREGRESSIVE  MODELING 


The  Autoregressive  Process 

Without  lack  of  generality  wc  consider  a  scalar  random  process  u{  with  zero  mean 
and  variance  ou~.  A  second  process  {  Xt }  is  defined  to  be  an  autoregressive  process  (AR)  of 
order  m  if 

Xj  =  “j  Xt_j  +  -  -  -  +  “in  *t-m  +  “t  (?) 

For  the  first  order  case,  Le..  m  -  .1,  we  have 

*t  =  al  Xt-1  +  ut  & 


Therefore  by  successive  substitution 


Xt  ~  ar  iai  ^-3 +  ^-2)  +  ai  “t-i  +  «t 

Then  Xt  can  be  expressed  as 

xt  =  ut  +  ji  “t-i  +«r  ut-2  + . 

That  is.  X|  is  expressible  as  an  infinite  order  moving  avenge  process. 

For  the  general  order  case  we  define  the  operator 

Bk  Xf  =  Xt-k  for  all  k 
Then  Equation  (1)  becomes 

fl-o,  3-...-om  Bml  X,  =  u, 

Set 

fTB)  =  M-s.  B  -  .  .  -  a,,,  B®]”1  (3) 

=  l+p,B  +  ...+0mBn,  +  ...  (say)  (4) 


F-l 


4 


We  find 


X,  =  f(B)  ut 


It  then  follows  that  E  X(  =0  and  the  variance  is  finite  as  soon  as  £  f%~  is  convergent. 
Vale-Walker  Equations 

Let  us  assume  the  process  to  be  stationary.  Then  multiplying  Hquation  (1)  by 
Xj_j.  and  taking  the  expectation,  we  find  for  each  k>0 

f(k)  =  aj  f(k-D  + . +  «„,  ?(k-m)  (6) 

where  j'  ts  the  autocorrelation  function  of  the  process.  The  set  of  Equation  (6i  is  the  well 
known  Yule-H’aiker  equations. 


The  Yule- Walker  equations  are  particularly  useful  in  estimating  parameters  of  an  AR 

process,  that  is,  to  estimate  the  paiameters  aj . a,,,  ,  given  the  observations 

Xj  .  X->  .  .  .  .  ,  Xj^r  ,  (say)  of  the  process  f  Xj|_  It  is  evident  from  Equation  (6)  that  we 
can  solve  for  the  ctj's  by  substituting  the  sample  autocorrelation  functions  into  the  firs:  m 
Yule-W  lker  equations.  Let  us  next  consider  two  recursive  methods  of  solving  for  the 
coefficients  of  the  AR  process;  one  for  the  scalar  (Durbin’s  Method)  case  and  the  other  for 
the  vector  case. 


1 

at 

a 


Durbin’s  Scalar  Recursive  Algorithm 

In  fitting  a  stochastic  (stationary  )  process  with  an  autoregressive  process,  lire  order 
of  the  process  N  must  be  decided.  The  criteria  may  be  the  least  squares  fitness  to  the  mvdcl. 
where  lire  mean  squared  error  terminates  the  process  when  it  becomes  smaller  than  seme 
predetermined  threshold. 


The  conventional  method  of  solving  for  the  coefficients  is  to  form  N  simultaneous 
Yulc-Walkcr  equations  and  solve  fo:  the  Aj’s  by  matrix  inversion.  However,  with  this 
method,  each  time  the  order  of  the  AR  process  is  increased,  an  NXN  matrix  inversion  is 
required.  A  more  efficient  method  is  to  use  the  Durbin’s  Recursive  Algorithm 
described  below. 


The  Akorithm 


Given  the  N-l  coefficients  :  . 


kjs»_|  jsf_j)  for  the  N-l  order  AR 


process,  the  N  coefficients  1  A\,  j . A>»  >^)  of  the  order  AR  process  is  obtained 

from  tin.  two  recursive  equations. 


or  in  brief  notation 


rN-1  ::  PN-1  *  aN  -  1  °r  R-(N-n  "  PN-1  A  “(N-l) 
where  the  coefficients  of  Aj^_j  aie  known.  Note  that  ?  matrix  is  positive  definite. 

The  ma'rix  form  of  the  Nth  order  Process  is  thus  given  by 


(10) 


Rn  =  PN:  -AN  01) 


This  can  be  rewritten  as 


An  A  <aNl . aN.  N-?) 


Then,  we  have 


rN  -1  =  PN  -i  *  An  +  R-<N-1>  *  ~N.N 
PS  =  r-4N-1)An  +  aN.Nfl0 


From  (13) 

aN.N^0  =  PN  "  ?  ~»N-1)  *  An 
From  Equations  (12  and  10) 

An  =  (RN_j  -  R^n'-U  £<n.n  1 

=  aNM  '  °XN  **^1  RhN_j, 


=  Av^,  -  aN  N 


(12) 

(13) 


(14) 


(15) 


Substituting  Equation  (15)  into  Equation  (14)  yields 


aN.N^O  ~  PN  -  Vk-1)  *  ANM  +  ®N,N  *  R-(N-1)  *  A-(N-I) 


and 

PS  ~  R-tN-I)  *  AN-i  PN  "  aN-1  *  R  \-] 

aN.N - - - "  = - - - 

•°0  -  R-(N-1)  *  \N-1)  Pq  ~  an_i  *  Rn_s 


(16) 


From  Equation  (15)  and  we  form 

A\  =  aN-1  -  “N.N  *  A-<N-1) 

aN,N  =  0  +  aN,N 
ana  i!  fellows  that 


Equations  (16)  end  (17)  are  the  recursive  equations  of  Equations  (7  and  8).  The  mean 
squared  error  is  obtained  from 


N  ^  N 

Oyf-Q  R  ~  2  QjX;t-i)l-  =  ElOPt)  -  2  a;X(t-i))  *  X(t)l 

ir  r  * « 


i=i 


i=I 


which  yields 


Hence 


°S~  =  PO  ~  %  ai  Pi 


N 

I 

5”  i 


°N~  -  PQ  -  an  *  rn 

Substituting  Equation  (17)  for  Ajj  above  yields 

7  ~  -w 

°N~  =  PQ  ^AN-I  RN-i  -  “N.N  A  -<N-j)  R\_|  +  “N.N'PN* 
=  *0  -  AN-j  r-N-l  '  °N.\5  lpN  '  A  -{X-D  *  RX-l> 

=  p0  *  AN-1  *  RN-1  -  °NLN  (PN  -  AN-1  *  R-<N-1)) 


ilS) 


( 19i 


F-5 


and  from  Equation  (16) 


°N~  =  (>0  ~  aN-1  RN-i  ~  ‘  0>o  ~  AN  -i  RN-1> 

=  “  AN-I  *  RN-1)  (i  -  aX,N~' 

=  »N-I~n  -  ttN  K")  <20> 

which  is  the  result  of  Equation  (9) 

Matrix  Recursive  Algorithm 

For  the  general  case,  the  multivariate  stochastic  process  can  be  modeled  by  a  Vector 
autoregressive  process  Z(t)% 

Consider  the  following  Nth  order  vector  autoregressive  process  Z(t) 

N 

Z(t)=  £  AXjZ(t-i)  +  0(t)  (21) 

i=l 

where 

Z(t)=  fzj(t)  zi(t) . Zp(t) 

;  pxp  matrix  coefficients  of  the  Nth  order  AR  process. 

U(t)pxi  vector  of  independently  distributed  zero  mean  white  noise  process, 
and 

~  denotes  transpose 

The  matrix  coefficients  Aj  can  be  computed  from  the  matrix  Yule-Walker  equations 
N 

Km  =  Z  AN4  m  =  1 . X  (22) 

i»l 


F-6 


where  Km  is  a  covariance  matrix  (hermitian)  given  by 


I  L-c 

—  £  Z(n+k)  *  Z*(a)  ■  =  0.  i.  2  .  .  . 

k=l 


K 


IE 


Superscript  *  denotes  complex  conjugate  transpose 


The  mean  squared  error 
from  Equation  <21 ) 


of  the  AR  process  may  be  obtained 


Sn  =  K°" 


N 

2 

i=j 


ANj  Ki 


(25) 


Real  "rocess 

The  matrix  recursive  algorithm  solving  for  the  Aj^  j  is  derived  below  for  a  real 
valued  multivariate  stochastic  process.  For  complex  case,  see  the  next  subsection  or.  complex 
processes. 


The  mean  squared  error  2.^  of  Equation  ? 23)  combined  with  the  \  Yule-Walker 
Equations  of  Equation  (22i  can  be  expressed  in  matrix  form  as  follows 


ll-  _AN.l.  _AN.2. - _AN.NJ 


K0 
K-l  K, 


1  kn 


0 


K 


N-l 


LK-N  K-<N-1)  K0  J 


*[SV 

N 


or  in  brief  notation 

H.  ~an1  *  *n  =  >2 2  0 . oi 

N 

For  N+lth  order,  we  pic1:  an  such  that 


"■  -4n-  01  -  I  £a.  0. 

N 


(24) 


(25) 


F-? 


0.  o^l 


which  is  the  representation  of  Equation  (24)  and 


Kn+i  -AS  ,  KS  -AN  2  KN_!  -----  AN.NKJ  '  0  -  K0  -  aN 


or 


N 

aN  =  ^N+l  "  2  AN.i  *  KN+1 
i=l 


(26? 


Wc  note  that  if  aN  -=  0.  then  AN+1  =  fA,\-  0)  and  the  desired  solution  is  obtained. 
If  this  does  not  happen.  aN  must  be  forced  to  zero.  In  order  to  do  this,  we  introduce  the 
backward  (adjoint)  equation. 

N 


Z(t-N)  =  2  BNj  Z(t-S  i) 
t=! 

b 

10-  -3N.N.  -BN.N-1.  -  -  -  -Bn.1.  *  KN+1  =  ^N.G.  -  •  •  .0.  I  ^ 


or  in  brief  notation 


!G.-BN  II  kn+1  *  l.?N.o.  -  - .  0.  Ij 


which  is  a  representation  of 


and 


l-BN  11  *  kn  =  10 . 0.  £ni 


N 

=  KhN+1)  *  2  ®N.i  *S-N-1 

i=l 


We  now  form  a  weighted  combinat:on  of  Equations  (2?)  and  <2r). 

([I.  -An.  01  +  QNa  (0.  -Bn  Ilj  ’  kn+i 


(27) 


(28) 


1  £3  +  0 . 0.  x#  +  0^  £n1 

N 


(29) 


or 


By  setting  the  last  element  of  the  vectoi  on  the  RHS 
<*N  +  "  0 

Qsa  =  -o'n  i  S*r*  =  0 


(30) 


then  we  have  the  solution  for  AN+1  in  the  form  of  Equation  (24)  for  N  provided  we  hav 
Bn.  and 

II.  -Ay+jl  =  U-  -Ay.  ol  +  Qfc  10.  I]  13 


and 


X  =  2*  +°N 

N+l  N 


(32) 


In  order  to  obtain  the  By  and  £  .  we  perform  a  similai  operation  as  in 
Equation  (29)  as  follows  ^ 

HO.  -By.  I!  +  QNb  II.  -Ay.  01}  •'  Kn+1 


=  I%  +  QNbSN°  • 


0. 


2'n 


<*fe 


i  QN 


.) 


133) 


from  which  we  set 

0N  +  Q>I  Zy-=  0 

or 

<&=-eN  iSy  r* 

and  it  follows  that 

l-gN+1'  11  =  l°-  -In-  +  °N  1K  “=N-  °’ 

vb  =  yb  +  Ofcory 
^N+l  WN  N  N 


F-9 


J 


■-  *  'I 


a 


jt 


r4 

1 

J 

4 


3. 

The  matrix  recursive  algorithm  is  stated  as  follows:  Given  <AN;  X^  ftN>  30(1 

2  ■»  b 

(3N:  Xb:^j)  for  the  Nth  order  AR  process.  (AN+1:  X  :0tN+l>  “«*<&*1:  2  : 

_!  ^  *  N+t  .  ; 


?N+1)  for  the  N+lth  iteration  is  obtained  from  the  following  matrix  recursive  equations 


a  b  -1 

la)  Q^:-  -Ojx  l  Xn  1 


2a)  U-  -^j]  =  (I.  -An.  0]  +  10.  -Bn.  I) 


3a)  X 


'n+i  =  £N  +qn 


IVTj 

a)  aN+j  =  Kn+2  -  Xj  aN+1J  *  KN+2- 


b  2 

!b)  QN=-^n  I  INi 


2b)  (-|N>5.  li  =  |0.  -BN.  II  +  Qjs-  U-  ~^N-  °1 


3i”  Xn+i  '  Zn  +  qn  ‘  "n 


%+i  =  K  -(N+2)  “  §  BN+12  *  Ki-N-2 


The  initial  conditions  are 


a  °  IT 

XQ  =  K0  »  XQ  :  og  =  Kj  and  aQ  =  K,  and  £&  -  K_j  -  K, 


wc  note  that  «N  =  0y;  •  *--c*  ar  av'iioint  relationship 


jg 


1 


I 


Complex  Process 

In  dealing  with  a  complex  random  variable  /  -  x  +  jy.  we  note  that  z  is  in  fact 
comprised  of  bivariate  random  variables  x.  y  .  and  hence  a  complex  scalar  stochastic 
process  must  be  treated  as  a  rea!  vector  of  dimension  two.  For  a  complex  N  dimensional 
process,  the  multivariate  dimension  is  2N. 

As  described  earlier,  in  the  process  of  generating  an  AR  model.  if,,  matrix 
coefficients  ol  A-  must  be  determined  when  the  .tochastii  process  is  complex,  special  care 
must  be  taken  ;p  formulating  tne  matrix  equations. 

lints  t!;c  following  rule  must  iv  observed 

i)  Coven  3  NxN  tomolex  matrix  M.  decompose  this  main v  into  the  real 
imaginary  parts. 

M  =  R..1M)  +  i  US  Ml 


2i  Define  a  2N\2\  real  mairx  M  wh 


here 


M  = 


Substitute  M  for  M 

The  complex  form  of  the  matrix  VuleAV-.lker  F.sjuation  of  l  q:tat:r>r  1 22 >  .ar 
be  rewritten  as 

\ 

— Nj  '  br.-.-i 


K-  *  I,  -'Xa  • 

can  be  rewritten  as 

\ 

Re‘Km»  ♦  i  WC,!  *  ♦H„IAXJil-  <RciK:1._:P  j 

?RciAx.ii  -  ReiKM-ji  -!i.-,AN.i»  *  Vu^m-i** 

•  •  RclKm-ii  *  UAyl  • 


F-l  1 


MiSSION 

of 

Rome  Air  Development  Center 


HAVC  plans  and  executes  xeseardi,  development,  test  and 
selected  acquisition  programs  In  support  cl  Command,  Contact 
Communications  and  Intelligence  (C* z)  activities.  Technical 
and  engineering  support  icithln  areas  technical  competence 
Is  provide d  to  BSD  Program  Q$ilces  ( PCs }  and  other  BSD 
elements.  The  principal  technical  mission  areas  are 
comvjtications,  el ?ctroKog netic  guidance  ted  control,'  sur¬ 
veillance  oi  ground  and  aerospace  objects ,  inteiilgzrcz  data 
collection  and  handling,  information  system  technology,, ^ 
ionospheric  propagation,  solid  state  sciences,  micrcuave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


/ 


